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PREFACE 


A sound mathematical training for the modern engineer 
is a sine qua non for attaining new heights in all aspects 
of engineering practice. One of the areas of mathematics 
that plays a big role in the mathematical education of 
the engineer is vector analysis, which is now invariably 
included in the curriculum of higher mathematics in 
engineering colleges. 

The present collection of problems in vector analysis 
contains the required minimum of problems and exercises 
for the course of vector analysis of engineering colleges. 

Each section starts with a brief review of theory and 
detailed solutions of a sufficient number of typical prob- 
lems. The text contains 100 worked problems and there 
are 314 problems left to the student. There are also a 
certain number of problems of an applied nature that 
have been chosen so that their analysis does not require 
supplementary information in specialized fields. The 
material of the sixth chapter is devoted to curvilinear 
coordinates and the hasic operations of vector analysis 
in curvilinear coordinates. Its purpose is to give the reader 
at least a few problems to develop the necessary skills. 

The exposition in this text follows closely the lines 
currently employed at the chair of higher mathematics 
of the Moscow Power Institute. 

The present text may be regarded as a short course in 
vector analysis in which the basic facts are given without 
proof but with illustrative examples of a practical nature. 
Hence this problem book may be used in a recapitulation 
of the essentials of vector analysis or as a text for readers 
who wish merely to master the techniques of vector 
analysis, while dispensing with the proofs of propositions 
and theorems. 

In compiling this problem book, the authors made 
extensive use of material in published courses of vector 
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calculus and collections of problems. Many problems 
were made up by the authors themselves. 

This collection of problems is designed for students of 
day and evening departments at engineering colleges and 
also for correspondence students with a background of 
vector algebra and calculus as given in the first two years 
of college study. 

We would like to express our sincere gratitude to 
Professor V. P. Gromov (the Moscow Krupskaya Pedagog- 
ical Institute), Professor A. V. Efimov and Associate 
Professors I. M. Petrov, B. I. Fridlender, and V. N. Zem- 
skov (the Institute of Electronics) for their thorough 
scrutiny of the manuscript of the book and for valuable 
tuggestions and remarks that were made full use of in 
she final editing. 


M. L. Krasnov 


A. I. Kiselev 
G. I. Makarenko 


Moscow-Dubna, 1977. 


CHAPTER | 


THE VECTOR FUNCTION 
OF A SCALAR ARGUMENT 


Sec. 1. The hodograph of a vector 
function 


Definition 1. A vector r is said to be the vector function 
of a scalar argument ¢ if each value of the scalar taken 
from the domain of admissible values is associated with 
a definite value of the vector r. This can be written as 
follows: 

r=r (t). 


If the vector r is a function of the scalar argument ft, 
r=r (t), 


then the coordinates z, y, z of the vector r are also func- 
tions of ¢: 


z=z(t), y=y(t), 2=2 (2). 


Conversely, if the coordinates of the vector r are func- 
tions of t, then the vector r itself is also a function of t: 


r=z(tjity(t)jt+2()k. 


Thus, specifying a vector function r (¢) is the same as 
specifying three scalar functions z (t), y (t), z (t). 

Definition 2. The hodograph of the vector function 
r (t) of a scalar argument is the locus described by the 
terminus of the vector r (t), as the scalar t varies, when 
origin of the vector r (t) is fixed at a point O in space 
Fig. 1). 

The hodograph of a radius vector r = r (t) of a moving 
point is the trajectory L of that point. Some other line L, 
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(Fig. 2) is the hodograph of the velocity v = v (t) of that 
point. Thus, if a material point (particle) is in motion 
around a circle with constant ve- 
A locity, | v | = constant, then its 
hodograph of velocities is like- 
tl wise a circle with centre at O, 

4 and with radius equal to | v |. 
Example 1. Construct the ho- 


0 B dograph of the vector r = ti+ 
ty + ?k. 
Fig. 1 Solution. 1°. This construc- 


tion may be carried out by using 
points and setting up a table: 


ci i ol 4 


t 0 
o| i+j+k 31+ 3j + 9k | 4i + 4j + 16k 


r 2i-+ 2j+-4k 


2°. Alternative solution. Denote by 2, y, z the coordi- 
nates of vector r; we have 


z=, y=ut, z= ff, 


Eliminating the parameter ¢ from these equations, we 
get equations of the surfaces y = xz, z = 2’, the line L 


Fig. 2 


of intersection of which is what defines the hodograph of 
the vector r (t) (Fig. 3). 


Sec. 2] The limit and continuity of a vector function 441 


1. Construct the hodographs of the following vectors: 
(a) r= 2i-+ ¢2j — ?k. 


_ Met. 2, 
(b) r= Gyape i+ Gea. 


(c) r=cost-i+sint-j+k. 


Fig. 3 


(d) r= t+ 2 Pj4+— bk. 


_ 2ti+2tj+ (¢2—2)k 
(e) r= ar oo . 


Sec. 2. The limit and continuity 
of a vector function of a scalar argument 


Suppose a vector function r = r (¢) of a scalar argument 
t is defined in some neighbourhood of the value ¢, of the 
argument t, except perhaps for the value ¢, itself. 

Definition 1. A constant vector A is said to be the 
limit of the vector r (t), as t—> to, if for any e >O there 
is a § >O such that for all t ~t, that satisfy the con- 
dition | ¢ — ty | < 6 the following inequality holds true: 


Ir(t) —Al<e, 
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As in the case of ordinary calculus, we write lim r (¢) = 


t+to 

= A. 

Geometrically (Fig. 4), this means that the vector 
r (t) tends, as ¢ > fy, to the vector A both in length and 
in direction. 

Definition 2. A vector 
a (t) is said to be infini- 
tesimal, as t— to, if 
a(t) has alimit, as t >t, 
and that limit is zero: 


lim @ (t) =0, 
(+t 
or, what is the same 
thing, if for any e >0 
Fig. 4 there exists a 6 >0 such 
that for all ¢ ~ t, that 
satisfy the condition |t—t, |< 6, the inequality 
| a(t) |<e holds true. 
Example 1. Show that the vector @ (¢) = ti + sin ¢j is 
infinitesimal when t > 0. 
Solution. We have 


ja (t)| = |ti+ sin tj] <|t| + |sint| <2 |#]. 


From this it is evident that if for every ¢ >0O we take 
& = e/2, then for | — 0 | < 6 = e/2 we have |a (t) |< 
< e. By the definition, this means that @ (¢) is an infini- 
tesimal vector when t— 0. 

2. Show that the limit of the modulus of a vector is 
equal to the modulus of its limit (if the limit exists). 
3. Demonstrate that for a vector function r (t) to have 
a limit A, as tt, it is necessary and sufficient that 
r(t) be representable in the form 


r(t)=A+a(t), 
where a@ (¢) is an infinitesimal vector when t — fp. 


4. Show that if the vector functions a (t) and b (t) have 
limits as t—> ft, 


lima(t)=A, limb(t)=B, 
tt 


t{+to 
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then their sum a (t) + b(t) and their difference a (t) — 
— b (t) also have limits as t—> ?¢,, and 

lim [a2 (t) +b (t)) =A+B. 

t+to 


5. Let 
lima(#)=A, limb(t)=B, 
toto 


t+to 
Prove that 
lim (a (¢), b (t)) =(A, B), 
t~to 
where (a (t), b (¢)) is a scalar product of the vector func- 


tions a(t) and b (t). 
6. Let 


r@=z@ity@jtz@k, A=ait+ajt+ 
+ a,k. 
Show that if lim r(t)=A, then 
to 


lim z(t)=a,, limy(t)=a,, limz(t)=as. 
t+to toto ttg 


Find the following limits: 
7, lim (= 14 Set j+ek). 


: eat i t 


e sint 
9. lim (= i+cost- it+az): 


10. lim ( sine jy Shest yt), 


11. i (£55 i+p> j+2k). 


Definition 3. A vector function r =r (t) defined in 
some neighbourhood of the value t = ft, is said to be 
continuous when t = t, if 


lim r (¢) =r (5). 
tte 
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In other words, r = r (¢) is continuous for ¢ = ¢, if for 
every € >O there is a 6 >O such that for all ¢ that 
satisfy the condition |[t—t,|<(5 the inequality 
|r (t) — r (to) |< e holds true. 

The hodograph of a continuous vector function of a 
scalar argument is a continuous curve. 
12. Start with the familiar inequality |a—b]|> 
=> ||a]— |b || and demonstrate that the continuity 
of a vector function implies the continuity of its modulus. 
Is the converse true? 
13. Show that if a (¢) and b (t) are continuous for t = fo, 
then the vector function a (¢) + b(t) is also continuous 
for t = ty. 
14. A vector function a(t) + b(¢) is continuous for 
t = ty.. Does it follow from this that the vectors a (t) 
and b (t) are also continuous when ¢t = ty? 
15. Prove that if a(t) and b (é) are continuous vector 
functions, then their scalar product (a(t), b (¢)) and 
vector product [a (t), b (¢)] are also continuous. 


Sec. 3. The derivative of a vector 
function with respect to 
a scalar argument 


Suppose a vector function r = r (t) is defined for all t 
on the interval (to, ¢,). Take some value ¢ € (éo, ¢;), then 
give ¢ an increment A¢ such thet t + At € (fo, ¢,) and find 
the corresponding increment Ar = r (t + At) —r (t) in 
the vector function r (t). Now consider the ratio Ar/At. 

Definition. If, as At — 0, the ratio Ar/At has a limit, 
then that limit is called the derivative of the vector 
function r = r (t) with respect to the scalar argument ¢ 
for a given value ¢ of the argument and is denoted as 


dr (t)/dt or r’ (t) or r (t). Thus, 


a =lim a lim 
t At—+0 At—+0 


r(t+At)—r(d) 
At . 


In this case the vector function r = r (t) is said to be 
differentiable. 
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16. Show that if the vector function r =r (z) has a 
derivative for some value ¢ of the argument, then it is 
continuous for that value ¢. 

The derivative of a vector function r (t) of a scalar 
argument ¢ is a vector directed along the tangent to the 


ar 
Mat 


Fig. 5 


hodograph of the original vector at the point under con- 
sideration (Fig. 5). The vector dr/dt is in the direction 
of the terminus of the vector r (¢) as it moves along the 
hodograph when the parameter ¢ increases. 

Suppose r =r (é) is the radius vector of a moving 
point. Then the vector v = dr/dt is the velocity vector 
of that point. 

Suppose 


ri) =z(t)ity@®itz(@)k, 


where the functions z (¢), y (¢), 2 (t) are differentiable at 
the point ¢. Then there exists dr/dt for that value of ¢ and 


dr __ dz dy ., dz 
ae ar ae ae k. (1) 


Example 1. Find dr/dt if r = ia cost + jbsint (the 
point is moving in an ellipse). 
Solution. From formula (4), 
dr 


a= —iasint+ jbcost. 
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By analogy with the differential of a scalar function, 
the differential of a vector function r = r (t) is a vector dr 
defined by the equality 


where df = At is the increment in the scalar argument ¢. 
As in the case of scalar functions, 


Ar = dr + a-At, 
where a = a (t, At) 0 as At—>0. 


Basic rules for differentiating a vector function 


Assume that all functions being considered (both scalar 
and vector) are continuous and differentiable. 
1°. If ¢ is a constant vector, then de/dt = 0. 
2°. The derivative of a sum of vector functions is equal 
to the sum of the derivatives a the summands: 
d(a(t)+b(t)) 4 ab 
dt =F +r 
3°. Suppose a vector function a(t) is multiplied by 
a scalar function m (t) of the same scalar argument. Then 
d(ma) _ py da, dm 
de ae a 


GP (0) + (Se) 


5°. srl — (= ’ b]+[a. al 


(In this formula, the order of the factors a and b in the 
right-hand member must be the same as that in the left- 
hand member.) 

Let us prove formula 4°. We set 9 (t) = (a(t), b (t)). 
Give ¢ an increment At; then, by the distributive prop- 
erty, we have for the scalar product 


Ag = 9 (¢ + At) — 9 (t) = (a + Aa, b + Ab) — (a, b) 
= (Aa, b) + (a, Ab) + (Aa, Ab), 


$= (42, b)+ (att) +(42, a8). @ 
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It is given that the functions a (¢) and b (¢) have deriva- 
tives for the value ¢ of the argument and, hence, are 
continuous for that value of ¢. Therefore 


fe_ da tim Sh = and lim Ab=0. 


areo Of dt? agg At de * At—+0 

Passing to the limit in (2) as At > 0, we obtain 
d(a,b) da ab 
at ar >)+(a a): 
17. Given r=r(t). Find the derivatives: 
d d d ad d 

@TP), bae(n gz). Os[rz]- 
18. Prove that if the modulus |r | of the vector function 
r = r(t) remains constant for all values of t, then dr/dt_Lr. 
What is the geometrical meaning of this fact? 
38: Prove that if e is a unit vector in the direction of 
the vector E, then 


[E, dE} 


[e, de] = TB 


20. Suppose 
u = uy (z, y, 2, t) i + Us (z, y, 2, t) } + Us (z, y, 2, t)k, 
where U;, Us, Us are continuously differentiable functions 
of their arguments, and z, y,z are continuously differentia- 
ble functions of ¢. Show that 

du _ gu, du dz ou dy | mu ae 
a atts at ay att as de 
21. Find the trajectory of motion for which the radius 
vector r(¢) of a moving point satisfies the condition 
dr/dt = {a, r], where a is a constant vector. 

The derivative dr/dt of the vector function r (¢) of a 
scalar argument is a vector function of the same argument. 
If there exists a derivative of dr/dt, then it is called the 
second derivative and is denoted d*r/dt*. Generally, 

cy Cs ae 
dim de \ den-l }» Ee eats, 
22. Given the radius vector of a point moving in space: 
r{asint, —acost, bt?} 


(t is time, and a and b are constants). Find the hodographs 
of velocity and acceleration. 


2-910 
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23. Given: r = acos wi + bsin wi, where , a, b are 
constants. Prove that 


(1) [r, 3] =(wa, by, 
(2) r+ o% =0. 


24. Show that if r = ae®t + be-°!, where a and b are 
constant vectors, then d*r/dt — w*r = 0, 

25. Show that the modulus of the differential of the 
radius vector of a point is equal to the differential of 
the length of the arc described by the point. 

26. Suppose a = a (u) is a vector function of a scalar u, 
where wu in turn is a certain scalar function of the basic 
scalar ¢. Assuming a (uw) and u = u (t) to be differentiable 
the necessary number of times, find an expression for the 
derivatives of the composite function da/dt, d*a/dt*. 


Sec. 4. Integrating a vector function 
of a scalar argument 


Definition 1. We will say that the vector function A (t) 
is the primitive of the vector function a (t) whent, << t< 
< t, if A (t) is differentiable and 


x =a(t), t€ (to, ty). 


Definition 2. The collection of all primitive functions 
of a(t) is termed the indefinite integral of the vector 
function of a scalar argument a = a(t). As in integral 
calculus, the indefinite integral of a vector function is 


denoted by the symbol j. and we have 
j a(t)dt =A(t)-+C, 


where A (¢) is one of the primitive functions of a (t), and 
C is an arbitrary constant vector. 

The following properties hold true for integrals of 
vector functions: 


lage faa (t) dt =a a(?) dt (ais a numerical constant). 
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2°. { (a(t)+b(t)) dt= jac) dt + j b(t) dt. 


27. Show that if ¢ is a constant vector, and a (t) is a 
variable vector, then 


j (c, a(t)) dt= (e, j a(t) at) F 


j le,a()]at=[e, Jao at). 
If 
a (t) = a, (t) i + a, (t)j + ay (¢)k, 
then 


face dt =i j a, (t)dt+j j a, (t) dt +k j a,(t)dt. (1) 


That is, the integration of a vector function reduces to 
three ordinary integrations. 

Example 1. Find the indefinite integral for the vector 
function a(t) = icost + je + k. 

Solution. According to formula (1), 


fawdt=i j costat+j | etat+k { dt = 
=isint—je*+kt+e, 


where ¢ is an arbitrary constant vector. 
Find the integrals of the following vector functions: 
i k 
28. a (t)=teli+ sin? t-j|— aay. 
ti a: 
rar + tej+cost-k. 


30. a (t) =cos tel"! .§ —t cos t?-j+k. 


29. a(t)= 


34. a(t)=— —tsin t-j+2'k. 


Let a vector function a (¢) be defined and continuous 
over a certain interval [t,, 7], which is the range of the 
argument ¢. 

Definition 3. We define the definite integral of a vector 
function a (t) on the interval [t), 7) as the limit of the 


2° 
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vector integral sums 
n-1 


= P2 a (Tr) At,, TR c [t, tn+1] 


as the length At of the largest of the subintervals [t,, ty ‘al 
(k =0, 1, ..., m—1) into which the interval [t,, 7 
is partitioned tends to zero: 
T n—-1i 
faw dt=lim > a(t,) Ath. 
At+0 
to h=a0 
The following formula holds true: 
T 
J a() dt=A(T)—A (6), (2) 
to 
where A (t) is some primitive for the function a (¢) on the 
interval [t,, T 


a (t) = a, (t)i + ay (t)j + a, (¢)k, 
then 
T 


T T T 
fatdr=i j a, (t) dt-+j { a, (t) dt +k j as(t) dt. (3) 
to to 


to to 


m/2 
Example 2. Compute j a(t) dt, where a(t) micost— 
0 


—jsin?¢. 
Solution. By virtue of formula (3), 
n/2 0/2 n/2 
\ a(t)dt=i| cost at—i sin? t dt 
0 0 0 
- [e/2 ft — sin 2e \ [2/2 
=isine|" —j (5-3) =1-+ij. 


Compute the following integrals: 


x 
32. J ater, where a=sin?¢ cos ¢-i+cos?¢sint-j+k. 
0 
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iets jet/2 


33. \ a(t)dt, where a=—,—+-5 + ke’. 


34. ttt 


35. 


Ct ay Che Qe Cte ae 


a(t) dt, where a= 3m cos nt-i——_ 4 2tk. 


a(t) dt, where a=(2¢-+)i+¢sint-j+ ak. 
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Example 3. An electric current J flows upwards along 
an infinite wire that coincides with the z-axis. Find the 


vector H of the magnetic Coen ear te 
field intensity set up by 4 
this current at an arbitrary 
point M (2, y, z) of space 
(Fig. 6). 

Solution. We consider a 
sufficiently small element 
PP, = dt ofthez-axis. By 
the Biot-Savart law, the 
intensity dH of the mag- 
netic field set up at point 
M by the current flowing 
through element df of the 
wire coincides in direction 
with the vector product 
(dg, r,], where d{ = PP,, 
|d§ |=dg, r1=PM (see 
Fig. 6). By this same law, the 
modulus of the vector dH is 


I YN 
|@H| =a sin (dG, r,) dC, 


Fig. 6 


PM(Z,Y,2) 


/~ 
where (d§, r,) is the angle formed by the vectors d§ and rj. 


Since 
: /N 
[[4§, ry]| =rydE sin (df, v4), 
we can write 
I 
dH =—~ (de, r,]. 


(4) 
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In order to obtain the desired vector H at the point M 
we have to sum all vectors dH pertaining to distinct 
elements PP, of the wire, that is, we have to integrate 
the expression (4) over the whole z-axis: 


+c 


H= j = (a, r). (5) 


We have 
r, = OM — OP. 
But 
OM = zi + yj + 2k, OP = Ck, 
and therefore 
r, = zi + yj + (2 — O)k 
so that 


ry = [ry] =V 22+ y? + (2— 0)? = V p?+ (20), 


where p = Vz + y? is the distance of point M from the 
axis of the wire. 

For the vector product [df, r,] we have 
ijk 
00 a 
zy 2z—C 


jd§, ry] = = —iydC+ jrdg, 


and formula (5) takes the form [the point M (z, y, 2) is 
fixed, J = constant] 


+o 


To compute the integral on the right-hand member of (6), 
make the substitution 
4 Shep 79%" 1 


~ _ pdt 
C—z=ptant, d= - 
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We then have 
+00 nf2 


a = j pat 
| ero 9 STOR OF tant 
co 1 
4 n/2 
=p \ cost dt= —> 
1/2 


Thus, the intensity vector H of the magnetic field is 
in our case given by 


T , : 
B= (yi eat 
or 
2 
H=—y Il, rj, 


where I = J-k is the current vector, r is the radius vector 
of point M (z, y, 2) of the field, and p is the distance of 
M to the axis of the wire. 

Example 4. The motion of an electron in a homogeneous 
magnetic field. 

1°. Suppose a magnetic field H is set up in some region 
of space; let it be constant in magnitude and direction 
(a homogeneous field). Suppose at time ¢t = fy, an elec- 
tron enters the field with an initial velocity v,. Determine 
the path the electron wil! take. 

Solution. First suppose the vector vo is perpendicular 
to H and that the initial position of the electron is at 
point M,. Choose the origin O at an arbitrary point of the 
plane P passing through M, at right angles to the vector H 
(Fig. 7). Let the initial radius vector OM, be ro, let r be 
the radius vector of the electron at the current instant of 
time ¢, and let v be the instantaneous velocity at that 
instant. The basic differential equation of motion is 

ar 
It will be recalled that the force F acting at time ¢ on the 
electron by the magnetic field is 


F = —e,[H, v], 
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where é@, is the absolute value of the electron charge. 
Thus, 


dtr 
m r= eolv, H). (7) 


At every instant ¢t, the force F is perpendicular to the 
direction of velocity and to the direction of the field H; 


Fig. 7 


at every instant, it will force the electron to deviate from 
a rectilinear path and to describe a certain curvilinear 
trajectory. 

Let us rewrite (7) as 


dv dr 
mar =eo|[ ar. 1] 


and integrate from t, to ¢ with respect to t. This yields 


my — mv, = é,lr, H) — eglry, H) 
or 
mv = eglr, H] + (mvo — eo [ro, H)). (8) 


Now choose a coordinate origin O’ such that the term in 
parentheses in the right-hand member of (8) vanishes, that 
is, so that 

€9 (ro, H] = mvp. (9) 
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From (9) it follows that the initial vector r, must be 

perpendicular to the vector v, and must lie on the straight 

line M,K, which is perpendicular to the plane of the 

vectors v, and H. By virtue of (9), the modulus of the 

vector rp must satisfy the relation 
éoltol-|}H|]=m|vol, 

whence 


_ mlVvol 
\Tol = eo [HI ° (10) 


This determines the position of the new origin O’. Relative 
to this origin, equation (8) is rewritten thus: 


mv = ey [r, HJ (11) 
or 
m & —e,{r, H). (12) 


From equation (41) it follows that the trajectory of the 
electron is a plane curve lying in the plane P because at 
every instant the vector v is perpendicular to H. Now take 
the scalar product of both sides of (12) by r: 


m (r. Gp) =e0(r. Ir. HI). (13) 
The mixed product in the right-hand member of (413) 


is zero, so that 
dr 
(r. =) =0, 


£ (r2) =0 or + (r2)=0, that is, r?=constant. 
This is the equation of a circle lying in the plane P with 
centre at the chosen point O’. The radius of the circle is 
found from formula (10) since the initial point M@, must 
also lie on that eircle. Thus, we finally have 
r= T= (10) 
Thus, if an electron enters a homogeneous magnetic 
field H with an initial velocity v, at right angles to H, 
then it will describe, in that field, a circular trajectory 
lying in the plane P perpendicular to H and passing 
through the initial point. The radius of the circle is 


whence 
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given by formula (10) and its centre O’ lies on the straight 
line perpendicular to the plane of the vectors v, and H; 
note that a rotation from v, to H must be seen from point 
O' as a counterclockwise rotation. 

From (10) it is evident that the radius r, of the circle 
is inversely proportional to | H |. Thus, the greater the 
intensity of the magnetic field, the greater the curvature 
of the trajectory. 

From formula (11), 


mv = ey lr, Hi, 


it is clear that if r is constant in modulus and is all the 
itme perpendicular to H, then also the velocity v of the 


- 


Fig. .8 


point will be constant in magnitude, 
|Jv|/ == constant, 


so that the electron is in uniform motion in the orbit. 
The period of revolution 7 is 


T = 270 On (14) 
Vo 


af Mw 
eo |H] * 
This formula does not involve the initial velocity vo. 
Thus, irrespective of the initial velocity v, which is 
perpendicular to H with which the electron enters the 
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homogeneous magnetic field H, it will perform a single 
orbital revolution and always in the same time 7. 

2°. Now suppose an electron enters a homogeneous 
magnetic field H with some initial velocity V that is 
not perpendicular to the vector H. This velocity may 
then be resolved into two components: the vector v, 
at right angles to the field, and a vector v, parallel to the 
magnetic field. 

From the formula 


F =e[V, H] = ey Ivo, HI 


it is evident that the “twisting” force F is given only by 
the perpendicular component v, and that it imparts to 
the electron a rotational motion about the circle (centred 
at O’) discussed above. As for the other component v,, 
the electron will retain it by inertia and, besides having 
a uniform circular motion, it will have a rectilinear and 
uniform motion in the direction of H with a velocity 
v, = |V|cosa. The combination of these motions 
yields a helical curve with axis parallel to the vector H 
and passing through the point O’ (Fig. 8). 


Sec. 5. The first and second derivatives 
of a vector with respect to the arc length 
of a curve. The curvature of a curve. 
The principal normal 


Consider a curve L in space. On it, choose a point M, 
as the origin and also choose a direction along L that will 
be regarded as positive. For a parameter, take the arc 
length s reckoned from M, of the curve (Fig. 9). Then 
the radius vector of a point M of the curve is 


r=r (s). 


With that choice of parameter, 


where t° is a unit vector directed along the tangent to 
the curve L in the direction of increasing values of the 
parameter s. 
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If the vector r is given by the coordinates 
r= zi+ yj + zk, 
then 
dz ., dy,, d 
Va Fit ei+gk 
and 
dz \2 dy \2 dz\2_ 
V (E+ (#)'4 (E) =" 
Since | 1° | = 1, the vector dr%/ds is orthogonal to the 


vector +°. 
The modulus of the vector dr°/ds is 


dx? 
| |=x. 
Here, K is the curvature of the curve L at the point M. 


The straight line having the direction of the vector 
dv°/ds and passing through the point M of the curve is 


Fig. 9 


termed the principal normal of the curve at the point M. 
Denoting the unit vector of that direction by n°, we have 
ee Rik (1) 


ds 
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The inverse of the curvature of a curve at a given point 
is called the radius of curvature of the curve at that point 
and is denoted by R: 


1 
R=7. 


Thus, formula (1) may be rewritten as 


dr dy in? 
ds3 ds R° 
From this, 
1 d*r 
K= R | as 
or 


K=4=) ($)'+(4)+(S). 


Using (2), we can compute the curvature of a curve at any 
point if the curve is specified by parametric equations in 
which the parameter is the arc length s. 

In the particular case of a plane curve, a circle of 
radius a, 


$ 

z=acos—, 

@ 

—_ . $ 

y=asine, 

we have 

d®*z 1 s dy ; 
aw et at ae oe NS. 


and formula (2) yields 


ae 2k rr ee 
K=—>= ar cos ata sin ao ae 


This means that the curvature of a circle of radius a is 
constant and is equal to the inverse of the radius of the 
circle. 
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If the curve Z is given by the vector-parametric equa- 
tion r = r (¢), where the parameter ¢ is arbitrary, then 


lar ae | 

4 at’ dt? 

Rr dr |8 (3) 
| 


Formula (3) permits computing the curvature of the 
curve at any point provided we have an arbitrary para- 


metric specification of that curve. 
Example 1. Compute the curvature of the helical curve 


r=acost-i+ asin t-j + htk. 
Solution. Since 


K= 


ae —asint-i+acost-j+hk, 
a 
ore —acost-i—asint-j, 
the vector product 
+ i j «kk 
(+. oe = —asint acosth 
—acost —asint0 
=sahsint-i—ahcost-j+ atk. 
Consequently, 
d d* area d = Le 
[S$ Grllaevere, |t|=Vaere. 
By virtue of (3), 
1 a 
K=3-aym 
or 
3 
Ratt = constant. 


Thus, a helical curve has a constant radius of curvature. 
Find the radius of curvature of each of the given curves: 

36. r = In cos t-i + In sin ¢-j + Y 2t-k. 

37. r = fi + 284. 
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38. r = 325i + (38¢ — ¢*) j + 2k for ¢ = 1. 

39. r= a (cost + t sin ¢)i + @ (sin t — ¢ cos 2) j 

for t = n/2. 

40. r = acosh t-i + a sinh ¢-j + atk at any point ¢. 


Sec. 6. Osculating plane. Binormal. 
Torsion. The Frenet formulas 


The plane passing through the tangent line and principal 
normal to a given curve J at a point M is termed the 
osculating plane at the point M. 

For a plane curve, the osculating plane coincides with 
the plane of the curve. 

If the vector r =r (t) has a continuous derivative 
dr/dt in the neighbourhood of a point ¢, and, besides, a 
Second derivative d’r (t,)/dt? such that 


dr (to) d*r (to) 
[ ae ae | #9, 


then at the point ¢ = ¢t) there is an osculating plane ‘o 
the curve r = r (t) whose vector equation is 


(erty, [ ul, S40 }) 0, 


where p = ¢ (é) is the radius vector of the current point 
of the plane. 

The normal to the curve at the point M, which normal 
is perpendicular to the osculating plane of the curve at 
that point, is called the binormal of the curve at the given 
point M. 

Denote by b® the unit vector of the binormal oriented 
so that the vectors t°, n°, b® form a right-handed trihedral 
(Fig. 10). Then 


bo? —1, b°=[t°, n°]. 
For the derivative db°/ds we get 


The vector db/ds is perpendicular both to the vector 1° 
and to the vector b?, that is, it is collinear with the vec- 
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tor n°. Set 
dbo | 
ds | T° 
We then have 
ab? _ to 
= 


The quantity 1/7 is termed the torsion of the given 
curve, and 7 is called the radius of torsion of the curve. 

The torsion of a curve 
is given by the formula 


1 dr dr dr 
pu (agra): 
where the symbol (a, b, c) 
denotes a mixed product of 
the vectors a, b, c, that is, 
(a, b, c) = (a, [b, cl). 

For the case where the 
curve is given by the vec- 
tor-parametric equation r = 
=r (t), we have 


(4,5 a 
dt’ dt?’ at) 


see a (") 


aia 
r= 


“dt? at? 


Example 1. Find the torsion of the helical curve 
r=acost-i+ asint-j + htk. 


Solution. We find the derivatives of the given vector: 


f= —asint-i+acost-j+hk, 
d*r : , F 
v7; ie —acost-i—asint-}, 


=, =asint.i—acost-j. 


d*r 
at® 
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The mixed product of these vectors is 
—asint acost h 
)= —acost —asint 0}=ah. 
asint —acost 0 


(+ da dr 
dt’ data’ dt 


In example 1, Sec. 5, we found that 
d dap 4/2 
[a> ar] f=e@+e. 


Using (1), we obtain for the torsion 


Ak 
T at+ha° 
Thus, the torsion of a helical curve is the same at all 


its points. 
Example 2. Write the equation of the osculating plane 
at the point ¢ = 0 of the helical curve 
r=acost-i+asint-j + htk. 


Solution. We find the values of the derivatives of the 
given vector and its derivatives dr/dt and d*r/dt? at the 


point t = 0: 
a 
r(0)=ai, HO ait hk, S50) 
Consequently (see example 1, Sec. 5), 


dr (0 d3r (0 : 
(SP. Set |= —ani+or. 


The vector equation of the osculating plane is 


dr (0 d*r (0 
(p—r(0), au ’ Ol) =0 


= —ai. 


or 
(pe — ai, —ahj + a’k) = 0. 

Since the radius vector of the current point of the osculat- 
ing plane p = zi + yj + zk, it follows that by passing 
to coordinate notation we obtain an equation of the 
desired plane in the form hy — az = 0. 

Formulas expressing the derivatives of the vectors 
1°, b°, n° are called Frenet formulas: 


ax o @bo 1 4 dn? 15 14,5 
—_—_ = n, —_— — n — Rv 7h. 


4 
ds R ds T° ds 


3-910 
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41. Write down the equation of the osculating plane at 
the point t = 2 of the curve 

r=ti—tj+ > ek. 


42. Write down the equation of the osculating plane at 
the point t = 0 of the curve 


r=e4+e"j+ VY 2tk. 
43. Find the torsion at the point ¢ = 0 of the curve 
r=e'cost-i + e' sin t-j + e'k. 
44. Find the torsion at any point ¢ of the curve 
r = acosht-i + asinh ¢-j + atk. 


CHAPTER II 
SCALAR FIELDS 


Sec. 7. Examples of scalar fields. 
Level surfaces and level lines 


Definition. If a value of a certain quantity is defined 
at every point of space or a portion of space, then we say 
that the field of the given quantity has been specified.v 

The field is termed a scalar field if the quantity in 
question is a scalar quantity, that is, if it is fully de- 
scribed by its numerical value.-: 

Examples of scalar fields are: a temperature field, an 
electrostatic field. 

Specifying a scalar field is accomplished by specifying 
the scalar function of a point M: 


u =f (M). 
If a Cartesian coordinate system zyz is introduced in 
space, we have 
u =f (x, y, 2). 
Geometrically, a scalar field is characterized by a 
level surface; this is a locus of points at which the scalar 


function of the field assumes the same value. The level 
surface of a given field is defined by the equation 


f(z, y, z) =C, where C = constant. 


In the case of a temperature field set up in a homo- 
geneous and isotropic medium by a point source of heat, 
the level surfaces are spheres centred at the source (this 
is a central-symmetric field). 


3e 
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In the case of an infinite uniformly heated wire, the 
level surfaces (isothermic surfaces) are circular cylinders 
whose axes coincide with that of the wire. 

Example 1. Construct the level surfaces of the scalar 


field 
u=2x+ 2y + 3z. 
Solution. The level surfaces are given by the equation 


xz+2y+3z=C, where C = constant. 
This is a one-parameter family of parallel planes, 
Example 2. Find the level surfaces of the scalar field 
u=2?+ y? — 2, 
Solution. The level surfaces are given by the equation 
2+ y?—2z?=C, where C = constant. 
For C = 0, we obtain a circular cone. For any C >0, we 
obtain a hyperboloid of revolution of one sheet with the 
axis coincident with the z-axis. For C <0, we obtain 


a hyperboloid of revolution of two sheets. 
Example 3. Find the level surfaces of the scalar field 


Zz 
u=aresin Vr : 


Solution. The domain of definition of the given scalar 
field is found from the inequality 


|<1, that is, O<—2— <1, 


worst, 


| vate ' 


whence 0 < 2? < 2? + y’. This double inequality shows 
that the field is defined outside a circular cone 2? = 
=z? + y’ and on it, with the exception of its vertex 
O (0, 0, 0). 

The level surfaces are found from the equation 


arcsin =C, where —F<C<F. 


4 
VzI+y 
That is, 2/Yz?+ y2=sinC or 2 = (x? + y*) sin? C. 
This is a family of circular cones located outside the cone 
2? = g* + y* with a common axis of symmetry, Oz, and 
a common vertex, O (0, 0, 0), at which the given field 
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is not defined; note that the cone itself, 22 = 2? + y?, 
is also included in the family. 
Example 4. Find the level surfaces of the scalar field 


u= ea ‘), 
where a is a constant vector and r is the radius vector 


of a point. 
Solution. Here, 


r= {z, y, 2} =zi+ yj + 2k 
and let 
a = {a;, a2, as} = a,i + a,j + ask 

Then the scalar product 

(a, r) = a,z + agy + agz. 
The equation of the level surfaces is 

e@N=C, C>0, 
whence 
(a, r)=InC 

or 

a,x + ayy + agz = In C. 
This is a family of parallel planes. 

Find the level surfaces of the following scalar fields: 


_ a yt | et 
45. u=Ttot+e: 
46. u=2?4 y2—z. 
47. ponte. 


z 


48, u=2y? + 922. 
49, u= 3*t2u-*. 
(a, 1) 

50. u= (ber) (a, b are constant vectors). 
Si.u=In|r!. 
52. u = ea.) (a, b are constant vectors). 

A scalar field is said to be plane if there is a plane 
such that in all planes parallel to the given plane the 
scalar field is the same. 
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If we take this plane as the zy-plane, then the scalar 
field is given by the scalar function 


u = f(z, y), 


which is to say it is not dependent on z. 
An example of a plane scalar field is the temperature 
field of an infinite, uniformly heated wire. 


y.-z y py 
c<0 
6>0 0>0 
zr 
c<0 
Fig. 11 


Geometrically, plane scalar fields are characterized 
by level lines; these are loci in which the scalar function 
has one and the same value. 

Example 5. Find the level lines of the scalar field 


use — 7, 
Solution. The level lines are given by the equations 
xz? — y? =C, where C = constant. 
When C = 0, we obtain a pair of straight lines: 
y=2, = —2z. 


For C 0, we obtain a family of hyperbolas (Fig. 11). 
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Find the level lines of the following plane fields: 
53. w= 2r—y. 


b4.u=in 2. 
2 


55.u=L. 
z 

56. u= ev'-v. 

57. Find the level lines of the scalar field u given im- 
plicitly by the equation 


utarlinu+y=0. 


Sec. 8. Directional derivative 


Suppose we have a scalar field defined by a scalar 
function 
u =f (M). 


In the field, take a point M, and choose a direction 
indicated by the vector 1. Then in the field take another 
point M so that the vector M,M is parallel to 1. Denote 
by Au the difference 


Au = f (M) — f (Mo) 


and by Al the length of the vector MM. The ratio Au/Al 
defines the average rate of change of the scalar field per. 
unit of length’ in the given direction. Allow the point M 
to move towards the point M, so that the vector MyM 
is always collinear with the vector I]. Then Al— 0. 

Definition. If, as Al > 0, there is a limit to the ratio 
Au/Al, then it is called the derivative of the function 
u =f (M) at the given point M, in the direction of | 
and is denoted by the symbol du/dl so that, by definition, 
we have 


elim 34 = tim LORE) MM IIL. 


at~o 44 aro 

This definition of a directional derivative is invariant, 
that is, it is not connected with any choice of coordinate 
system. 
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Suppose a Cartesian coordinate system has been intro- 
duced in space and suppose the function f (M) = 
= f (z, y, z) is differentiable at the point M, (Zo, Yo, 25). 
Then 


Ou = 


du x ou Ou 
© le =e lo, O88 + Ge cos B+ o> ial (1) 


M M, 


where cosa, cos, cosy, the direction cosines of the 
vector 


1 = a,i + a,j + ak, 
are found from the formulas 


aad a NB a! 
cosa=Ty cos B my ©oSv=TP> 
[1] = V a? +a? + a3. 
Ou Ou Ou Ae 
ee symbols Be hag? ay ai,9 Oe lad: signify that the par- 
tial derivatives are taken at the point Mo. 
For a plane field u = f (x, y), the directional derivative 
I at the point My (xo, Yo) is 


ou Ou 


“@ Im, oz 


Ou 
cosa+—|_ sina 2 
Mo = oy la , (2) 


where @ is the angle formed by the vector ! and the 
z-axis. 

Remark. The partial derivatives du/dz, du/dy, du/dz 
themselves are derivatives of the function wu in the direc- 
tion of the coordinate axes Oz, Oy, Oz respectively. 

Formula (1)—used to compute the directional derivative 
at a given point—holds true even when the point M@ 
tends to M, along a curve for which the vector } is the 
tangent line at the point M,. 

Example 1. Find the derivative of the scalar field 


u = zyz 
at the point M, (1, —1, 1) in the direction from M, to 
M, (2, 3, 1). 

Solution. We find the direction cosines of the vector 
MoM, = (1, 4, 0}, the length of which is |MoM,| = 17, 
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and we have 


cosa = cosB= Ta cosy =0. 


1 
Viv’ 
The values of the partial derivatives of the function 
u = zyz at the point M, (1, —1, 1) are 
ul 4 ml iy me 
“Oz |M, * @y Im, ”s@z 
Using formula (1), we get 


Ou 1 


@ lm, OU 47 


4 3 

a ya a 
Ou 

The fact that aE lng 


M, increases in the given direction. 
Example 2. Compute the derivative of the scalar field 


=—1. 


Mo 


> 0 means that the scalar field at 


u = arctan ry 


at the point M, (1,1), which belongs to the parabola 
y = 2’, in the direction of the curve (in the direction of 
increasing abscissas). 

Solution. The direction of 1 of the parabola y = 2? at 
the point M, (1, 1) is the direction of the tangent to the 
parabola at that point (Fig. 12). 

Suppose the tangent 1 to the curve at M, forms with 
the z-axis an angle a. We then have 


y’=2r, tana=y’ |x, = 2, 


whence the direction cosines of the tangent line are 
1 1 


cos 80 = ——————— = — 
y t+tan?a 5° 


cosB=sina 


_ tana _ 2 
Vittantta y5- 


The values of the partial derivatives of the given function 
u (z, y) at the point M, (1, 1) are 


1 1 du 1 1 


du 
mM, 1+z%yilm, 2° oy lm, 1+z%y?im, 2° 


“Or 
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Substituting them into (2), we obtain 
ou 1 1 1 2 


cae ee 235. 
a” 2 V5 2 V5 2yY5° 
Example 3. Find the derivative of the scalar field 


Fig. 12 


u = zz" + 2yz at the point M, (1, 0, 2) along the circle 
zr=14 = 


y=sint—1, 
z=2. 


Solution. The vector equation of the circle is of the 
torm 


r (t) = (4 + cost) i + (sin t — 1) j + 2. 
We find the vector + tangent to it at any point M to be 
os = —sint-.i+cost-j. 
The given point M, (1, 0, 2) is found in the 2zz-plane in 


the first octant and is associated with the value of the 
parameter t = n/2. At this point we have 


eae nm, : 
t|m,= — sin > -i+ cos zi=ai. 
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From this we obtain that the direction cosines of the 
tangent to the circle are equal to cos a = —1, cos fp = 0, 
cos y = 0. The values of the partial derivatives of the 
given scalar field at the point M, (1, 0, 2) are: 


Ou = 2 a =2s| =a 
oz |M, M, éy IM, M, 
Ou 
“Oz am, = (228 + 2y) My 


Hence the desired derivative is 


Ou 


Ou _ Ou 
al 


IM, ov |m 


=4.(—1)+4.044.0=—4. 


In the following problems, it is required to find, for 
the given functions, the derivative at the point 
M, (Zo, Yo: Zo) in the direction of the point My (ty, Yr» %)- 


PB u=Vrt+y+a2t, My(1,1,1), My (3, 2, 4). 
58 u=s2y+ 222-2, My(i,41, —1), Mi (2, —4, 3). 
iO. u = xeY + ye* — 22, M, (3, 0, 2), M,(4, 4, 3). 


61. u=——t,M o(4, 1), My (4, 5). 
62. Find the derivative of the scalar field 
u = In (2? + y’) 


at the point M, (1, 2) of the parabola y? = 4z in the 
direction of the curve. 

63. Find the derivative of the scalar field u = arctan y/z 
at the point M, (2, —2) of the circle x? + y? — 4z = 0 
along an arc of the circle. 

64. Find the derivative of the scalar field u = x? + y? 
at the point M, (zo, yo) of the circle z? + y? = R? in the 
direction of the circle. 

65. Find the derivative of the scalar field u = 2zry + y’ 
at the point (VY 2, 1) of the ellipse 22/4 + y?/2 = 1 in the 
direction of the outer normal to the ellipse at that point. 
66. Find the derivative of the scalar field u = 2? — y? 
at the point (5, 4) of the hyperbola z? — y? = 9 in the 
direction of thé curve. 


~~? 
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67. Find the derivative of the scalar field u = 
= In (zy + yz + 22) at the point M,(0,1,1) in the 
direction of the circle x = cost, y = sint, z = 1. 

68. Find the derivative of the scalar field u = 2? + y? + 
+ 2? at the point M, that corresponds to the value of 
the parameter t = x/2 in the direction of the helical 
curve z = Roost, y = Rsint, z = at. 


Sec. 9. The gradient of a scalar field 


Suppose we have a scalar field defined by a scalar 
function 
u=f (z, y, Z) 
where the function f is assumed to be differentiable. 
Definition. The gradient of a scalar field u at a given 
point M is a vector denoted by the symbol grad wu and 
defined by the een 


grad u=it+S = +k. (f) 


Using formula (1) of Sec. 8 for ‘he directional derivative, 
we have 


= (grad u, 1°), (2) 
where 1° is a unit vector in the direction of I, that is, 


9 = =icosa + joos B-+k cos y. 


Properties of a gradient 


- 1. The gradient is in the direction of the normal to the 
level surface (or to the level line if the field is a plane 
field). 

2. The gradient is in the direction of increasing values 
of the function of the field. 

3. The modulus of the gradient is equal to the largest 
directional derivative at a given point of the field: 


max St :\grad u| = Y (2 a) +(= )’ +(3 my. 


* The maximum is taken over all directions of | at the given 
point of the field. 


* 
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These properties yield an invariant characteristic of 
the gradient. They state that the vector grad u indicates 
the direction and magnitude of maximum change of 
a scalar field at a given point. 

Example 1. Find the gradient of the scalar field 


u=a2 — 2Qy + 3z. 
Solution. By (1) we have 
grad u = 1i — 2j + 3k. 


The level surfaces of the given scalar field are the planes 
xz — 2y + 3z = C; the vector grad u = {1, —2, 3} is the 
normal vector of the planes of this family. 

Example 2. Find the greatest steepness (rate) of rise 
of the surface u = x” at the point , 2,4 

Solution. We have 
grad u = yr¥-4ji + 2¥ ln zj, gradu |y = 4i + 4 In Qj, 

(FF) ns = |grad u| = 4 VY 1+ (In 2)?. 

Example 3. Find the unit vector of the normal to the 

level surface of the scalar field 


u=2t yt 2. 
Solution. The level surfaces of the given scalar field 
are the spheres 


Y+y+2=C (C>0). 


The gradient is directed along the normal to the level 
surface so that grad u = 2z-i + 2y-j + 2z-k defines the 
vector of the normal to the level surface at the point 
M (z, y, 2). For the unit vector of the normal, we obtain 
the expression 

gradu tityjtzko or 
[ered ul Vatpyipa 


jr] ° 
Example 4. Find the gradient of the field u = (a, b,r), 
where a and b are constant vectors and r is the radius 
vector of the point. 
Solution. Let 


n? = 


a = {4,, dg, ay}, b = ({b,, bg, bs}, r= {z, y, 2}. 
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Then 
@, ay ag 
u= by be b3 ‘ 
zy 2 
By the rule for differentiating a determinant" we have 
@, dz ay Q, a, ag Gy a, a3 
a= by by 63], = by be bs |, #2 = by by bs}. 
10 0 010 001 
Hence 
G2 a3/, | 44 as a; a, 
cai Ae leas) ve fas 
ijk 
=]@, a, d3|=[a, bj. 
by by by 


* Given a determinant D (t) whose elements a,; are differenti- 
able functions of t: 


1, (t) aig (t) «-- Gin (t) 
D(t)= y(t) gg (t) ..- Gen (t) 


Gn (t) Gng(t) »-. ann (t) 


Then the derivative of the determinant, D’ (t), is found from the 
formula 


ai, (t) Gig (t) «++ Gin (t) 
D’(t)=| 221) Gan (4) «++ Gan (t) 


Gn; (t) @ng(t) --- Gnn (t) 
Gy, (t) yg (t) ..- yn (t) 
4 Go, (t) Gg (t) --- aon (t) 


ec) 


Gn1 (ft) ng (t) --- Snn (2) 
Gy; (t) yg (t) --- ain (2) 
ee fy] 2 aaa) «++ Gan 


Gny(t) ng (t) --+ Gan (¢) 
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Example 5. Find the gradient of the distance 


where P (2, y, z) is the point of the field being studied 
and P, (29, Yo, 20) is some fixed point. 
Solution. We have 


Or ., Or or 
gradr=>o its its k 


— = 20) i+—yo) I+G—)k _ Yo 
V (220)? + (y— Yo)? + (2— 20)? 
which is the unit vector of the direction P,P. 


’ 


Fig. 13 


Example 6. Let us consider the scalar function 
= ry + To, 


where r,, r, are the distances of some point P (z, y) of the 
plane from two fixed points, F, and F,, of the plane. 
Solution. The level lines of this function are ellipses. 
We have (see example 5) 
grad (r, + r,.) =r? + ro. 
This shows that the gradient is equal to the diagonal of a 
rhombus constructed on the unit vectors of the radius 


vectors drawn to point P from the foci F, and F, (Fig. 13). 
Consequently, the normal to the ellipse at some point 
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bisects the angle between the radius vectors drawn to that 
point. 

Physical interpretation: a light ray coming from one 
focus enters the other focus. 

Example 7. Find the angle 8 between the gradients of 
the functions 


u=Var+y? and v=r+y+2V ay 


at the point M, (4, 4). 
Solution. We find the gradients of the given functions 
at the point M, (A, 1): 


__avityi | _ iti 
gradu |e = PT la, VE 


gradvimy=[(1+V/ 2) it (1+ 2) i||,,,= 2+24- 


The angle ® between grad u and grad v at the point M, 
is found from 


2 2 
Kod (gradu, gradv) ly, vat V2 _1 
[graduly -leraduly, = =64-2YR 
From this we have 
6 = 0. 
Example 8. Find the directional derivative of the 
radius vector r for the function u = sin r, wherer = |r |. 


Solution. By (2), the directional derivative of the 
given function of the radius vector r is 


du 7 
3r = (grad sinr, r°). (3) 
We find the gradient of the function: 

i — 


Pada Or. d(sinr) or 
Dae a ser a 


Q Ge 9 (sin r) k 


grad sinr = aoa ay Ce ae j 


_ d(sinr) Or. 
~~ dr oz | 


=($1+5 itg 5~k) cos r = 19 cos r. (4) 
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Substituting (4) into (8), we get 
aad =(r°cosr, r°) = (r°, r°) cosr =cosr. 

Example 9. Find the derivative of the scalar field 
u = f (z, y, z) at the point My (Zo, Yo, Zo) of the curve J 
specified by the system of equations 

f(z, y, 2) =a, 
@ (z, y, 2) 20 
in the direction of the curve. 

Solution. The direction of the curve / is given by the 
direction of its tangent vector t, which, by definition, 
is a vector that is tangent to the surface f (x, y, z) = a. 
The surface f (z, y, z) = a is a level surface of the given 
scalar field u =f (z, y, 2). Since 


} » @=constant, 


ae (grad u, 1°) = (gradu, t° 
al ) 


and the vector grad wu is perpendicular to the level surface 
f (z, y, 2) = a, it follows that grad u is perpendicular to 
the unit vector 1°, and therefore 

ou 

ol |M, 

Example 10. At the point M, (1, 1, 1), find the direc- 
tion of the greatest change in the scalar field u = zy + 
+ yz + 22 and the magnitude of that change. 

Solution. The direction of the greatest change of the 
field is indicated by the vector grad u (M). We find it 
thus: 

gradu (M)= (y+z)it@+2ji+yt+a2)k 
and hence grad u(M,) = 2 (i+j+k). This vector 
determines the direction of the greatest increase in the 


given field at the point M, (1,1, 1). The magnitude of 
maximum change of the field at this point is 


max 3¢ = | gradu(M,) | =2)3. 


69. Find the gradient of the scalar field u = ln (z? + 
+y? + 2) at the point M, (1, 1, —1). 

70. Find the gradient of the scalar field u = ze**+¥*+2? at 
the point O (0, 0, 0). 


4-910 


= (gradu, T°) |v, =0. 
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71. Find the angle m between the gradients of the func- 
tion w= arctan z/y at the points M, (1, 1) and 
M, (—1, —1). 
72, Find the angle @ between the gradients of the func- 
tion u = (x + y)e**” at the points M, (0, 0) and 
M, (4, 1). 
73. Find the angle @ between the gradients of the func- 
tions u = V2? + y* + 2 and v = In (z* + y? + 2) at 
the point M, (0, 0, 1). 
74. Find the points at which the gradient of the scalar 
field u = sin (x + y) is equal to i + j. 
75. Find the points at which the modulus of the gradient 
of the scalar fieldu = In VY 2? + y? + 2? is equal to unity. 
76. Let-u = u (z, y, 2) and v = v(z, y, z) be functions 
differentiable at the point M (z, y, z). Show that 

(a) grad (Aw) = A grad u, A = constant; 

(b) grad (wu + v) = grad u + grad v; 

(c) grad (uv) = v grad u + u grad v; 


(d) grad (=) = VeretuS oes v0. 
77. Show that 
grad u(p) = grad 9, 


where p = 9 (z, y, 2) is a differentiable function and 
u = u (qg) has a derivative with respect to 9. 
Find the gradients of the following scalar fields if 


r=zit+yjtzk, r=|r|=)V22+y?4 22, 


and a and b are constant vectors. 


78. u = Inr. 

79. wu = (a, r). 

80. u = (a, r)-(b, r). 
81. u = | [a, rl [?. 


82. Show that 
(grad u (r), r) = wu’ (r)-r. 
83. Show that 
{grad u (r), r] = 0. 


84. Let w = f (u,v), where u = u (z, y, 2), v = v (z, y, 2). 
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Prove that 

gradw= a gradu+ ot grad v 
if f, u, v are differentiable functions. 
85. Suppose G is a convex region in space (that is, a 
region such that if two points M and N belong to G, 
then the whole line segment MN lies in G). Let there be 
given in G a scalar field u (M) which at all points has a 
gradient that is continuous and bounded in G: 


}gradu(M)|< A, MEG, A = constant. 


Prove that for any points M and N of G we have the 
inequality 
Ju(V) —u(M)|<A|MN|. 


86. Find the derivative of the function u = 27/a? + 
+ y?/b? + 27/c? at an arbitrary point M (z, y, z) in the 
direction of the radius vector r of that point. 
87. Find the derivative of the function u = 1/r, where 
r= |r| in the direction of the vector 1 = cosa-i + 
+ cos B-j + cosy-k. Under what condition is the 
derivative equal to zero? 
88. Find the derivative of the function u = 4/r, where 
r= |r|, in the direction of its gradient. 
89. Find the derivative of the function u = yze* at the 
point M, (0, 0, 1) in the direction of its gradient. 
90. Find the derivative of the scalar field 
u =u (zg, y, 2) 
in the direction of the gradient of the scalar field 
v =v (zg, y, 2). 

Under what condition is it equal to zero? 
91. For the following scalar fields, find the direction and 
magnitude of greatest change at the nee points M,: 

(a) u(M) = Ji irda ie (1, 0, 0). 

(b) u (M) = zyz; M, (2, 1 


4° 


CHAPTER III 
VECTOR FIELDS 


Sec. 10. Vector lines. Differential equations 
of vector lines 


Definition 1. We say that a vector field is given if a 
vector quantity a = a (M) is specified at each point M 
of space or of a portion of space. 

If a Cartesian coordinate system is introduced in the 
space, then specifying the vector field a= a (M) is equiv- 
alent to specifying three scalar functions of the point 
P (M), Q@(M), R(M) so that 


a(M)=P(z,y,z)it+Q(z,y, zi + RG, y, 2k. 


Definition 2. A vector line of a vector field a is a 
curve at each point M of which the vector a is directed 
along the tangent to the curve. - 

Let a vector field be defined by the vector 

a= Pi+ Qj + Ak, 
where 
P=P(z, y, 2); Q@=Qcz, y, Z), R= RQ, y, 2) 
are continuous functions of z, y, z that have bounded par- 
tial derivatives of the first order. 

Then the differential equations of the vector lines are 
of the form 

ae _ wv _ as (t) 

Integrating this system of two differential equations 

(1) yields a system of two finite equations: 


Gi (zt, y, 2) = Cy, 93, y, 2) = Cy 
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which, taken together, define a two-parameter family of 


vector lines: 
’ ’ =C , 
1 (z,y 8 1 } (2) 
P2 (Zz, y, 2) = Co. 
If the conditions of the theorem of existence and 


uniqueness of a solution are fulfilled in a certain region 
G for system (1), then a unique vector line 


Ps (ZY, 2) =P (Lo, Yoo Zo)s 
P2 (2, Y, 2) =P2 (Zo, Yor 20) 
passes through every point M, (Zo, Yo: 20) € G. 
Example 1. Find the vector lines of the vector field 
a = [e, rJ, 
where c is a constant vector. 
Solution. We have 
e=ci+cj+tek, r= zi+ yj + zk 
so that 
ij k 
a=([c, r) =] cy Cy C3 
ry z| 
== (C22 — Cgy) ¥ + (Cyr — C42) § + (Cry — Cyr) ke. 
The differential equations of the vector lines are 
aa = ais = joa Z (3) 


Multiply the numerator and denominator of the first 
fraction by z, the second by y, the third by z and add 
termwise. Using a property of proportions, we have 

dz dy dz zdz+ydy+zdz 

C—O tye C“‘i<~i ;‘CS 
whence 
zdz+ydy+2dz=0 

and this means that 

r+y+2=A,, A, = constant >0. 
Now, multiplying the numerator and denominator of the 
first fraction of (3) by c,, the second by c¢,, the third by c, 


54 Vector Fields (Ch. IT] 


and adding termwise, we obtain 
dz dy dz __¢) dz+-¢, dy +e, dz 
0 


bas ._—4 
Cg2— Cgy CyZ— C12 €);¥— Cz 
whence 


ec, dz + ¢,dy +es,dz=0 
and, consequently, 
Cz + cgy + C32 = Az, Ay = Constant. 
The desired equations of the vector lines are: 
P+y+2=A;, 
ye Cy + Coy + C32 = Ap. 

These equations show that the vector lines are obtained 
via the intersection of spheres (having a common centre ut 
the origin of coordinates) 
with planes perpen- 
dicular to the vector ce = 
=¢i + cj + ck. From 
this it follows that the 
vector lines are circles 
whose centres lie on a 
straight line passing 
through the coordinate 
origin in the direction of 
the vector c. The planes 
of the circles are per- 
pendicular to the indicat- 
ed straight line (Fig. 
14). 
Fig. 14 Example 2. Find the 

vector line of the field 


a= —yi+ zj + bk, 


which line passes through the point (1, 0, 0). 
Solution. The differential equations of the vector 
lines are 


whence we find 
py =— cy Cy >0 
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or, introducing the parameter ¢, 
z=VC,cost, y=VCsint. 


In this case, the equation 


takes the form 
VC,costdt dz + 
Wiese eee 


whence we obtain 


Thus, the parametric equations of the vector lines are 
r= V C, cos t, 
y=VC,sint, (4) 
z= bt + C2. 
If we require that the vector line pass through the point 
(4, 0, 0), we will have 
1=VC,cost, 0=VCysint, 0=bt+Cp. 
The first two equations of this system are satisfied for 
t = 2kn, k = 0, +1,...and for C, = 1. Taking k = 0, 
we get t = 0 and the last equation of the system yields 
C, = 0. The desired vector line passing through the 


point (1, 0, 0) is 
zr=cost, 
y=sint, 


2= bt. 

This is a helical curve. 

Find the vector lines of the following vector fields: 
92. r = zi + yj + zk. 
93. a = a,i + a,j + ask, where a,, a,, a, are constants. 
94. a= (z—y)it+ («—2z)jt+ yy —2z)k. 
95. Find the vector line of the field 

a= ri — yj + 2k, 


which line passes through the point (?/,, —1/,, 1). 
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A vector field is said to be plane if all the vectors of a 
are located in parallel planes and the field is the same 
in each of the planes. 

If a Cartesian coordinate system rOy is introduced in 
one of the planes, the vectors of the field will not contain 
any components along the z-axis and the coordinates of 
a vector will be independent of z, that is, 


a= P(z, yit+Q, y)i. 
The differential equations of the vector lines of a 


plane field are of the form 


dz dy az 
P(z,y) Q(z, y) OO 


or 
dy _ Q(z, y) 
dz” P(z,y)* 
z= constant. 


From this it is evident that the vector lines of a plane 
field are plane curves lying in planes parallel to the 
zy-plane. 

Example 3. Find the vector lines of a magnetic field 
of an infinite current conductor. 

Solution. We will assume the conductor is in the direc- 
tion of the z-axis and that the current J flows in that 
direction. The intensity vector H of the magnetic field 
set up by the current is 


H=< 11, rl, (5) 


where I = /-k is the current vector, r is the radius vector 
of the point M (z, y, z), and p is the distance from the 
axis of the wire to the point M. Expanding the vector 
product (5), we obtain 


H=—2 alz 


or ita). 


The differential equations of the vector lines are 


Sec. 10] Vector lines 57 


whence 
r+ y= R, 
z=C. 
That is, the vector lines are circles with centres on the 
z-axis (Fig. 15). 
Find the vector lines of the following plane vector 
fields: 


96. a = zi + 2yj. 
97. a = zi + 2k. 
98. a = zi — yj. 
99. a = 2273 + 4yk. 
100. a = xi + yj. 
101. a = 2j — yk. 


The differential equations of 
the vector lines 


may be written as 


dz dy az 
awa? =O GaR 


or, in vector form, as 
$= a(M). (6) Fig. 15 


This form of the equations of vector lines turns out to be 
convenient in the solution of a number of problems. 
Example 4. Find the vector lines of the field a = [e, r], 
where ¢c is a constant vector. 
Solution. Applying (6), we get 
dr - 
at = [e, tl]. (7) 
Forming the scalar product of both sides of (7) by ¢ and 
using the properties of a mixed product, we find 


(e, +) = +(e, r)=0. (8) 
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Similarly, forming the scalar product of both sides of (7) 
by r, we obtain 
(r, F)=F(r, =0. (9) 
From equation (8) it follows that 
(c, r) = constant 
and from equation (9) it follows that 
(r, r) = constant. 


The vector lines are lines of intersection of the planes 
(c, r) = constant with the spheres r? = constant. 
Find the vector lines of the following vector fields: 
102. a = f (r)-r. 
103. a = (a, r) bo, where ao, by are constant vectors. 


Sec. 11. The flux of a vector field. 
Methods of calculating flux 


J. The flux of a vector field. Suppose we have a vector 

field 

a(M) = P(z, y, z)it+Q(z,y, 25+ Re, y, 2)k, 
where the coordinates P (z, y, z), Q (2, y, 2), R (z, y, 2) 
of the vector a (M) are continuous [the field a(M) is 
continuous] in some region G. Let S be a smooth or piece- 
wise smooth two-sided surface in which a definite side 
has been chosen (an oriented surface). 

Definition. The flux IT of a vector field a (M) through 
an oriented surface S is defined as the surface integral 
of the first kind, over the surface S, of the projection 
of the vector a (M) by the normal n (MM) to that surface: 


Ii = jJ prna dS = J (a, n°) dS, 


where n° is the unit vector of the normal n to the chosen 
side of the surface S; dS is the area element of the sur- 
face- S. 

In the case of a closed surface, we will always choose 
the outer normal n that is directed outwards from the 
region bounded by the surface S. 
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If a, B, y are the angles that the normal n forms with 
the coordinate axes Ox, Oy, Oz to the surface S, then 


the flux may be expressed in terms of a surface integral 
of the second kind: 


I = J (a, n°) dS = \J [P (z, y, 2) cosa 


+ Q(z, y, 2)cosB+ R(z, y, 2) cos y} dS 


or 
N= (a, n°) dS = P (zx, y, 2) dy.dz 
FFeoanas {I 
+Q(z, y, 2)drdz+ R(z, y, 2) dzdy, 
where 


cosa dS = dy dz, cos fh dS = dr dz, cos y dS = dx dy. 


Basic properties of the flux of a vector field 


(a) The flux reverses sign when the orientation of the 
surface is changed (that is, when the orientation of the 
normal n to the surface S is changed): 


ff @nas=— Jf (a, nyas, 


S+ s~ 


where S+ is the side of.the surface S on which the normal 
n is chosen, and S~ is the side of S on which the normal 
—n is taken (see [7)). 

(b) Linearity: 


\ (Aa+ pb, n°) dS =A \J (a, n°) dS +p jf (b, n°) dS, 


where A and p are constant numbers. 

(c) Additivity: if the surface S consists of several 
smooth parts S,, S,, ..., Sm, then the flux of the vector 
field a (M) through S is equal to the sum of the fluxes 
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of the vector a (M) through the surfaces S,, S,,..., Sm: 
m 
n= > ff (nas. 
het “S, 


This property permits extending the notion of flux 
to piecewise smooth surfaces. 


Fig. 16 


Example 1. Find the flux of the vector a = i through 
an area perpendicular to the z-axis and having the shape 
of a rectangle with sides 7 and 2 (Fig. 16) in the positive 
direction of the z-axis. 

Solution. According to the definition of the flux of 
a vector through a surface S, we have 


l= \5 (a, n°) dS. 


In our case, a = i, n° = i so that (a, n°) = (i, i) = 1. 
Taking into account that the area of the rectangle is 
equal to 2, we obtain 


A a 


Remark. If we had chosen the unit vector of the normal 
to the area S so that n° = —i, we would have got IT = —2. 
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Example 2. Compute the flux of the vector field 
a =r, where r is the radius vector, through a right cir- 
cular cylinder of altitude hk, base 
tadius R and the z-axis. 

Solution. The surface S consists 
of a lateral surface o,, an upper 
base o, and a lower base gs of the 
cylinder. By the additivity proper- 
ty, the desired flux II is equal to 
ll = 1, + MM, + My, where M,, I,, 
II, are the fluxes of the given field 
through 0,, 02, Os respectively. 

On the lateral surface 0, of the cyl- 
inder the outer normal n° is paral- 
lel to the zy-plane and therefore 


(a, n°) = (r, n°) = praor = R 


Fig. 17). Hence 
Nl, = j J (ann) aS = 2 j j dS = R.2nRh=2nRh. 
0 


On the upper base o, the normal n° is parallel to the 
z-axis and therefore we can put n° =k (see Fig. 17). 
Then 

(a, n°) = (r, k) = pro, r =h 
and so 
11, = { { (a, n°) dS =h | ( dS =h-nR2=nRh. 
[frome 

On the lower base o, the vector a = r is perpendicular 
to the normal n° = —k. Therefore (a, n°) = (r, —k) = 0 
and 


N= | ft, n°) dS =0. 


The desired flux is then 


n= i) (a, n°) dS = 3n Rh, 
s 
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Example 3. Find the flux of the vector field 
Tr 
2= TP 

through a sphere of radius R with centre at the origin 
of coordinates. 

Solution. Since the normal n to the sphere is collinear 
with the radius vector r, we can take n° = r® = r/|r |. 
Therefore 


(et ee 
(a, n°) = (Tar rt) Stee = Tee = Te 
But on the sphere S we have |r | = R, and so (a, n°) = 


= 1/R?. 
The desired flux II is 


n= (a, n') dS = ar dS = 4n 
s s 


since the area of the whole sphere S is equal to (as = 


s 
= 4x R?. 
104. Compute the flux of the vector a = 3j through an 
area having the shape of a triangle with vertices at the 
points M, (1, 2, 0), M, (0, 2, 0), M; (0, 2, 2) in the direc- 
tion of the coordinate origin. 
105. Find the flux of the vector 


a = ai + 6j + yk, 


where a, B, y are constants, through an area perpendicular 
to the z-axis and having the shape of a circle of radius R, 
in the positive direction of the z-axis. 
106. Find the flux of the vector a = r through the outer 
side of a circular cone whose vertex lies at the origin of 
coordinates; the base radius is equal to A and the altitude 
is h (the axis of the cone is along the z-axis). 
107. Find the flux of the vector a = f (|r |) r through 
a sphere of radius R with centre at the coordinate origin. 
II, Methods of computing the flux of a vector. 
1°. The method of projection onto one of the coordinate 
planes. Let an open surface S be projected one-to-one 
onto the zy-plane into a region D,,. In this case, the 
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surface S may be given by the equation z = f (x, y) and 
since the area element dS of the surface is 


dzdy 
|cosy|’ 


dS = 


it follows that computing the flux [I through the chosen 
side of the surface $ reduces to computing a double integral 
via the formula 


n= Jo aas= { f (a, 2”) dxrdy. (1) 
Ss ry 


|cosy | |z=s(x, y) 


Here the unit vector n° of the normal to the chosen side 
of the surface S is found from the formula 


d {z—/(z, y)) dx dy 
n° = + gra = + ° 
| grad [z—f(z, y)] | a\2, (a\z.,’ 
V (zy +(H) 
(2) 
and cos y is equal to the coefficient of the unit vector k 
in formula (2): 


1 
cos py = tf (3) 
Of \2 of \2 
V (4) +(4,) +4 
If the angle y between the z-axis and the normal n° is 


acute, then in (2) and (3) the plus sign is taken, and if 
the angle y is obtuse, the minus sign is taken. The symbol 


(a, n°) 
| cosy | |z=f(x, y) 


signifies that f (z, y) must be substituted for z in the 
integrand. 

If it appears to be convenient to project the surface S 
onto the yz- and zz-planes, use is made of the following 
formulas to compute the flux IT: 


(a,_ 0°) 


| | cosa | 
i h4 


dy dz (4) 


z= Ply, z) 
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or 


1a (a, 0°) 
Il == \ RT ene, » 27 (5) 
x2 

Formula (4) is used when the surface S is projected 
one-to-one into the region D ,, of the yz-plane, which means 
that it may be given by the equation z = 9 (y, 2); cosa 
is found as the coefficient of the unit vector i in the for- 
mula 


d([z—@{y, 2)] ~ “éy az 
ee en noel eh 2 LS) a en Se. ay 
| grad [z—(y, 2) | VY ap \2 ap \2 
1+ (3F) +(3r) 
y Oz 
That is, 


4 
cosa=+ 7/7 op \2, | Opie . 

V 1+ (32) +(2) 

The plus sign is taken if the angle a between the z-axis 
and the normal n° is acute, and the minus sign if the 
angle a is obtuse. 

Formula (5) is used in the case of a one-to-one projec- 
tion of the surface S onto the zz-plane; in this case, 
S may be specified by the equation y =» (z, 2) and 
then 
ap ay 5 


grad {y—w(z, 2)) _ Page an 


EE SO —_ ;SxWu 
| grad [y—‘p(z, z)] | V 1+ (2) +(2) 3! 


cos Bf is equal to the coefficient of the unit vector j in 
this formula, that is, 


cos Bp = + 


n? = + 


1 
v4 (28)'4 (2B) 
V (2) +(2) 
If the angle 6 between the y-axis and the normal n° is 
acute, we take the plus sign, and if the angle 6 is obtuse, 
we take the minus sign. 


Remark. When the surface S is specified implicitly by 
the equation D(z, y, z) = 0, the unit vector of the 
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normal 
n° =icosa +jcosB + kcosy 


is found from the formula 


o_ grad © (z, y, 2) 
net | grad M (z, y, 2) | 


ao om 
a ita 


aD \2, / @ \2, { 2’ 
V (2) + (a) (ae) 
where the sign on the right is determined by the choice 
of the normal to the surface S. 

To compute the flux II of a vector field a through a 
surface S, it is necessary to project the surface one-to-one 
onto one of the zy-, zz-, yz-planes; this is possible if the 
equation ® (z, y, z) = 0 is uniquely solvable with respect 
to z(2 = f(z, y), yy¥=V, 2)) or r(c= ly, 2)) 
respectively. Then take advantage of one of the formulas 


(1), (4), (). 
Example 4. Find the flux of the vector field 


a= (¢ — 22) i+ @ + 3y+2j+ Obc+y)k 


through the upper side of the triangle ABC with vertices 
at the points A (1, 0, 0), B (0, 1, 0), C (0, 0, 1). 
Solution. The equation of the plane in which the triangle 
ABC lies is of the form z+ y-+2=1, whence z = 
= 1— 2 -— y. The triangle ABC is projected one-to-one 
onto the zy-plane into the region D,,, which is the triangle 
OAB (Fig. 18). 
It is given that the normal n° to the plane in which th 
triangle ABC lies forms an acute angle y with the z-axis 
and so we take the plus sign in (2) and obtain 
0 sred(etute—t) of gy ty 
e=Tged tyre) ys tysit ys © 
We find the scalar product 


(a, n°) =(2—22) Fe + (2+ By +4) Fe + (2 +H) oe 


_ Tz+4y—2z 
V3 — 


ad 
j+On 
aoe dz 
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From formula (6) we find that cosy = 1/) 3 >0 and, 
hence, 


dS = aa =V3drdy. 


Using formula (1), we aaipate the desired flux: 


= \ j (a, n°) dS = \ (7z-+ 4y — 2) |2at-2-y dz dy 
s 


xu 
{-x 


= | | @r+5y—1)ardy= dz | (82+ 5y—1)dy=3. 
Dey 0 


Example 5. Find the flux of the vector a = y*j + zk 
through the portion of the surface z = 2? + y* cut off by 


Fig. 18 


the plane z = 2. The outer normal is taken with respect 
to the region bounded by the paraboloid. 

Solution. The given surface (a paraboloid of revolution) 
is projected one-to-one onto the zy-plane into the circle 
D,, (Fig. 19). We find the unit vector of the normal n° 
to the surface S: 


—2zi—2yj+k 


grad (z— 2? — y?) 
V 429+4y?+1 | 


[grad (z—z?—y3)] 


n? = + 
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It is given that the normal n° forms an obtuse angle y 
with the z-axis and therefore the minus sign is taken in 
front of the fraction. Thus, 

no 22+ ike 
V 428+ 4y?+9 ” 


whence 
1 
ON Vara ~ 
and so 
dzd To TLD LA 
=TcoeT = V 422+ 4y? +1 dz dy. 
We find the scalar product 
(a, n°) ie 2y3—z 


V 428+ 4yt Hi * 
The desired flux is, by (4), equal to 
= j j (a, n°) dS = iF) (2y3—z) |amaey.je dx dy 


s se 


= j j (2y8 — y2—~ x?) dz dy, 


Dry 
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The domain of integration D,, is a circle of radius R = 
= ) 2 with centre at the coordinate origin. Introducing 
polar coordinates zr = pcos gq, y = psin g, we have 


I] = j j (2p3 sin’ p— p?) p dp dp 
Dey 
an V2 
=| dp j (2p sins @ — p*) dp = —2n& i —2n. 
0 0 


Example 6. Find the flux of the vector field a = 
= i — j + zyzk through the circle S obtained by cutting 


AY 
\ 


\ 


RW 
\\ 


A\ 


Z 


the sphere z? + y* + z*< R? with a plane y = z. Take 
the side of the circle facing the positive portion of the 
z-axis. 

Solution. Since the plane y = z is perpendicular to the 
zy-plane, the circle S lying in that plane is projected onto 
the zy-plane into the line segment A,A, and so the one- 
to-oneness of the projection is disrupted. The circle S 
is projected one-to-one onto the other coordinate planes. 
Projecting the circle onto the zz-plane, we obtain a region 
D,, bounded by an ellipse (Fig. 20). The equation of the 
ellipse can be found by eliminating y from the system 
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of equations 
ety? += Ri, 


y=a, 
whence 
Qat 4 ot = RP or 24 a1, 
3 


It is given that the normal to the circle S forms au 
obtuse angle f with the y-axis (see Fig. 20) and so we 
take 

n= —grad (y—z) =i—j, 


From the latter equation we have cos Bp = —1//2 < 0. 
The area element dS of the circle is equal to 
dS =—2% _V2drdz. 


™~ TeosB | 
We find the scalar product: (a, n°) = V2. 
The desired flux is, using formula (5), 
= = ( = . aR? = 9 2 
ll 1 2dr dz=2 \ drde= 2.7 = V2 Rim 
since the area Q of the region D,, bounded by an ellipse 
with semi-axes a = R//2 and b = R is equal to 


Q= \ dz dz=nab= 7 . 


Example 7. Compute the flux of the vector a = 
= zi+ yj + 2k through the outer side of the lateral 
surface of the circular cylinder 2? + y? = R*® bounded 
by the planes z = 0 and z = H (H > 0). 

Solution. The given cylinder is projected onto the 
zty-plane into a line, namely, into the circle (Fig. 21) 


re+y= R2, 
z=0. 


£24 
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We will therefore project the cylinder onto the other 
coordinate planes, for instance, the yz-plane. Since the 
cylinder does not project one-to-one onto the yz-plane, we 
take advantage of the additivity property of the flux of 
the vector and represent the de- 
sired flux IJ as a sum of 
fluxes: IT=II,-+I],, where I], is 
the flux ofthe field through the 
portion S, of the cylinder located 
in the region where y > 0, and 
Il, is the flux of that field 
through the portion S, of the 
cylinder located in the region 
where y < 0. On S, we have 


ne = Zul , (a, n°) = z ae =R 


and so 
Il, = RdS=R dS = RS, 
[fawn 
Fig. 21 where § is the area of the por- 


tion S, of the cylinder. Since 
S = 1RH, it follows that TI, = 1 R?H. 
On S, we again have 
_ tity) _ zty? 
ee eS 


and so 


ll, = \j RdS = RS =n. 


5S, 


The desired flux is I] = 2nR°H. 

Remark. The solution is made simpler if we introduce 
curvilinear coordinates z = Reosg, y= Rsing, z= 2 
on the cylinder (see item 3° below). 

To find the flux of the vector field a = P (z, y, z)i + 
+ Q(z, y, z)i + R (z, y, z)k through the surface S 
specified by the equation z = f (x, y) by the method of 
projection onto a coordinate plane, it is not necessary to 
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find the unit vector of the normal n°, and we can take 
the vector 
7) of . 
n= + grad[z—f(z, yJ=+ (—H1—-Fi+k). 


Formula (1) for finding the flux IIT becomes 
N= j j (a, n°) dS = §j (a, n) dz dy. (7) 
Ss Dry 


In similar fashion we obtain formulas for computing 
the fluxes through surfaces given by the equation z = 
= 9(y, z) or y= (z, 2). 

Formula (7) is written thus in coordinate form: 


n= +f {—Piz.y, f@ mi —Olz, y, fe, Wiz 


Dry 


z=J(x, y) 


+ Riz, y, f(z, y))} dx dy. 
Example 8. Compute the flux of the vector field 


a=zityj+V2+y—1k 
through the outer side of the hyperboloid of one sheet 
z= Vz? + y*—1 bounded by the planes z = 0,z = V3. 
Solution. The given surface is projected one-to-one onto 
the zy-plane into the region D,, bounded by the circles 
P+ y=, and 2 ts } 
z=0, z=0. 
We find the outer normal n: 
= —~Vriyp—i)= oat had | 
n= + grad (2—V2?+y2?—1)=+ ( Vato +k). 
Since n forms an obtuse angle y with the z-axis (Fig. 22), 
we take the minus sign and, hence, 
mens | 
Vite 
We find the scalar product 


eee el Lee 24 i_s. 
MOMsivempyaq ce eo >= yet 
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Using formula (7), we obtain 


n= {f (@nyas—{f dedy 


7at,3_4° 
: Bee | 


Passing to the polar coordinates z = pcosq@ and y = 
=p sing, we have 


2 ¢ 
= pdpdp _ p dp 
n= {J Vp?—1 = | de | V p§—1 
Dry 0 1 
=2n V p?—1|, =2V3n. 
Example 9. Compute the flux of the vector field 
a= yi+2j+ 2k 
through a closed surface bounded by the cylinder z?+ y* 
= R? and the planes z = z, z = 0 (2 > 0). 


Fig. 22 


Solution. The surface S is piecewise smooth and so we 
take advantage of the additivity property of a flux re- 
presenting the desired flux I as a sum of fluxes I], JI,, Hs 
through the smooth portions S, (semicircle z? + y*? < R?, 
0<z< R, z= 0), S, (part of the plane z = z), and 
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Ss (part of the cylinder z? + y? = R*): I = I, + I, + 
+ II,. Since S is closed, we take the outer normal to it 
(Fig. 23). 
(4) On S,, where z = 0, we have n® = —k and so 
(a, n°) = —z. 
This means that the flux 


ese jj za =—J J zazdy. 


s, S, 


Passing to the polar coordinates tr = pcos mg and y = 


=p sin g, we find 


= — j j pcos Gp dp dp 
" n/2 R 
= = eee | e2dp= — Ra, 


(2) On S,, where z = z, we have 
n = +erad (2 — z) = + (-i +k) 
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and since the normal n to S, forms an acute angle with 

the z-axis, we take the plus sign in the right-hand member. 

Thus, n = —i+ k and, hence, (a, n) = z — y. 
Projecting S, onto the zy-plane, we get the semicircle 


Day: OS T<KV R2—y?. 
Then by (6) we have 


T= | f @, n)|_aray, 


Dry 
and again passing to polar coordinates, we find 


1/2 R 
Nl, = j (cos p— sin @) dp j p2dp= - R, 
7/2 0 


(3) On S;, where z? + y* = R?, that is, on the lateral 
surface of the cylinder, we have 


o_ grad (z?+y8—R*) zityj — , zityj 
w= tT ed@ty—R| + Yap + R 
In this case, the normal n forms a right angle with the 
z-axis and therefore cos y = 0 and so the choice of sign 
in the right-hand member is arbitrary. Take the plus 

sign and then 


n= StH (a, n°) = ery 


and so 


It is impossible to project the surface S, (right cylinder) 
onto the zy-plane since it projects into a line, a semicircle 

“(the one-to-one nature of the projection will be upset). 
The same occurs when projecting onto the zz-plane. We 
therefore project the surface S, onto the yz-plane, onto 
which it is projected one-to-one into the region D,, 
bounded by the line 


a+ yes RP, 
Z=T.- 


Eliminating z from this system, we obtain the equation 
for the projection of this line onto the zy-plane: z? +- y? = 
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= R? is a circle. Since 


=|== (x50), 


reone=tecnme. | i)| =| (nj, i) |= | R 


we will have 


Naar | § feeetl Peat hag 
ws yz 
= {| SPayaeme | j y dy dz. 
Dy: yz 


Using the polar coordinates y = p cos » and z=¢ sin q, 
we get 


R 
N= 2 2 d do = 0. 
3 SJ pcos pp dp dp= feos oo 0 


Thus, 
2 2 3 
T= —3R+35aR +0=0. 


108. Compute the flux of the vector field a = yi + 2j + 
+ zk through the upper side of a triangle bounded by the 
planes 


ttytz=a, r=0, y=0, z=0. 


109. Compute the flux of the vector field a = xzi through 
the outer side of the paraboloid z = 1 — zx? — y? bounded 
by the plane z= 0 (z > 0). 

110. Compute the flux a the vector field a = zi + zk 
through the lateral surface of the circular cylinder y = 
= V R? — 2? bounded by the planes z = 0, z = h (h > 0)o 
111. Compute the flux of the vector field a = zi + yj + 
+ zk through the upper side of a circle cut out of the 
plane z = h (kh > 0) by the cone z = Vz" + y?. 

112. Compute the flux of the vector field a = 3zi — yj — 
— zk through the outer side of the paraboloid 2? + y? = 
= 9 — z located in the first octant. 

113. Compute the flux of the vector field 


a=@+yYityYt+a”)j+@r+xryk 
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through the portion of the plane z = 0 bounded by the 
circle z? + y? = 1 in the direction of the unit vector k. 
114. Compute the flux of the vector field a = yzi — zj — 
— yk through the total surface of the cone z* + y? = 2? 
bounded by the plane z=1 (0< z < 1). t 
115. Compute the flux of the vector field a = 2zi + 
+ (1 — 2y) j + 2zk through the closed surface bounded 
by the paraboloid 2? + 2? = 1— 2y (y > 0) and the 
plane z = 0 (z > 0). 
116. Compute the flux of the vector field a = zi + 
+ y*4j + 2k through the total surface of a pyramid 
bounded by the planes z+ y+2z2=1, r=0, y=0, 
Z2— UV. 
117. Compute the flux of the vector field a = zi + yj + 
+ zk through the sphere 2? + y* + 2? = R?. 

2°. The method of projection onto three coordinate planes. 
Suppose a surface S is projected one-to-one onto all three 
coordinate planes. Denote by D;,, D:,, Dy, the projec- 
tions of S onto the zy-, zz-, yz-planes respectively. 

In that case the equation F (z, y, z) = 0 of the surface 
S is uniquely solvable for each of the arguments z, y, z 
so that 


z=cz(y,z) y=y(z,z) z=2(z, y). 
Then the flux of the vector 
a = P(g, y; z)i+ Q(z, y, z)j+ Ro, Yy, z)k 


through the surface S, the unit vector of the normal of 
which is 


n° = cos a@-i + cos B-j + cos y-k, 
can be written thus: 


= \5 (a, n°) dS 


= j j [P (x, y, z)cosa+ Q(z, y, 2) cosB 
s 


+ R(x, y, z)cosy]dS. (8) 
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We know that 
dS cosa = +dy dz, 


dS cos B = +dz dz, (9) 
dS cos y = +dr dy, 
the sign in each of the formulas of (9) being chosen to 


coincide with the sign of cosa, cos f, and cos y on the 
surface S. Substituting (9) into (8), we get 


N= +] jPiey, z),Y, 2) dy dz + { } tz, y (z, 2), 2] dzdz 


Dyz Dxz 


cf in) R{z, y,2(z, y)|\dedy. (10) 
Dry 
Example 10. Find the flux of the vector 
a = zyi + yzj + zzk 


through the portion of the external side of the sphere 
z?+ y* + 2? = 1 located in the first octant. 
Solution. We have 


—Tered @+yete—)l ~ yarn tut 2k, 
whence, taking into accqunt that the surface S lies in 
the first octant, we obtain 


cosa@a=2z2>0, csB=y>0, cosy=z> 0. 


We therefore take the plus sign in (10) in front of the 
integrals, and putting 


P=2y, Q=yz, R=2, 


we obtain 
l= j j zy dy a) j yedrds+ | j azdzdy. (11) 
Dyz ez xy 


From the equation of the sphere z? + y? + z* = 1 we get 
z=2(z,y)=Vi—2—y, y=y(z, 2)=V 1—22—2, 
z=2(y,2)=V1—ye— 2. 


Substituting these expressions for z, y, z respectively 
into the third, second, and first integrals on the right of 
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(11), we get 


Tl = {—x2?—ydrd {—2x—22drd 
je" r—y rat \ Sav z?— 22 dzrdz 


fe §j yVi—y—2tdydz. (12) 
Dyz 


Let us compute the first integral on the right and pass to 
the polar coordinates t=pcosq@ and y =opsinp, 
where O< gp <i x/2, OS p <1. This yields 


I,= j j rVi—2—ydzdy= j j p? V1— p2 cos pdedp 
Dry Dry 
n/2 t 1 
= j cos pdg | eV 1—ptdp= j 2? V/1— p2 dp. 
0 0 0 
Setting p = sint, dp = cost dt in the last integral, we 
have 
nj2 n/2 
I,= | sin? t cost dt = + | sin?2tdt= =. 
The second and third integrals in (12) are computed in 
similar fashion and we obtain 


I= j j 2 Vi-2—2drdz=5,, 


Dxz 


hm | \eVir Bos ay dem a. 
yz 
The desired flux is 


N=Nt+h+l=F- 


118. Use the method of projecting onto all three coordinate 
planes to compute the flux of the vector field through a 
surface S. 
(a) a = 2i — 2j + yk; 
S is the upper side of a bounded portion of the plane 
3z + 6y — 2z = 6 cut out by the coordinate planes. 
(b) a=(et+yt2z)it@t+ytz—1)j—%; 
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S is the upper side of part of the planez + y+z2—1=0 
lying in the first octant. 


(c) a=(2—V y—2*) i+ 5+ (V y—a?—2z) k; 


S is the outer side of the paraboloid of revolution y = 
= x? + 2? bounded by the plane y = 4 and lying in the 
first octant. 

3°. The method of introducing curvilinear coordinates 
on a surface. In certain cases, when calculating the flux 
of a vector field through a given surface S it is possible 
to choose a simple coordinate system on the surface 
itself to compute the flux instead of projecting onto 
coordinate planes. 

Let us consider some special cases. 

Case (1). Suppose a surface S is part of the circular 
cylinder z? + y? = R? bounded by the surfaces z = 
=f, (z, y) and z = f, (z, y), and we have f, (zt, y) < 
< fa (z, y)- 

Setting 

z=Reosg, y=Rsing, 2 =z, 


we have for the given surface 
0< p<2n, f, (Reosg, Rsing)<az< 
< f, (R cos 9, R sin 9), 
and for the element of area dS we obtain the following 
expression (Fig. 24): 
dS = R dq dz. 
Then the flux of the vector field a through the outer 
side of the surface S is computed from the formula 
Qn !,(Rcos®, R sin 9) 
Il=R \ dp (a, n°) dz, (13) 
0 f,(R cos @, R sing) 
where 
o_ _ grad (x*+y2—R*)  __zi+yi 
~ [grad (z?+y8—R%)| RO 
Example 11. Find the flux of the vector 
r= zi+ yit zk 
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through the outer side of the lateral surface of the cir- 
cular cylinder z? + y*? = R? bounded by the planes 
z=0andz=H (H>0). 
Solution. Here we have 
0< px 2u; f,(Reosg, Rsing) =0 
f,(Reosg, Rsin 9) = 7 


Fig. 24 


Introducing curvilinear coordinates on the cylinder, 
we get 
z=Reoeosq, y=Rsing, 2=2 

Then the desired flux of the vector r is 

on H 

N=R j dg | (r, n°) dz. 

0 0 

But since 
r=2i+ yit+ zk = Reosqg-i+ Rsing- j+ & 

and the normal n° on the cylinder is 


n= Tu) Ass ght Reet —cosg. -i+sing: j 


it follows that the scalar product on the cylinder will be 
(r, n°) = Roos? ga + Rsin*? gp = R. 
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Finally, we obtain 
2n H 


T= R? j dg | dz— 2nF2H. 
0 0 
Example 12. Compute the flux of the radius vector 
r= 2i+ yj + zk 

through the lateral surface of the circular cylinder 
xz? + y? =1 bounded from below by the plane x + y + 
+z = 1 and from above by 
the plane z+y+2z=2. 


Solution. Here (Fig. 25) 
we have 


R= 1,f,(z, y) =1—2- 
— Hh my=2—r—y. 


Passing to coordinates on 
the cylinder 
z=cosg, y=sing, 


2=2, 
we get 
f, (z, y) =1—cos@—sin 9, 
fa (@ y) = 2—cos p — 
— sing. 


According to (13), the flux 
of the vector r is 


on 2-cos @-sin ® 
N= [ dy \ (r, n°) dz. 


0 1-cos g-sing@ 
But since on the cylinder z* + y* = 1 we have 
n°? = zi + yj = cos g-i + sin Q'j, 
it follows that 
(r, n°) = 2? + y? = cos* pt sin?’ g = 1 


6 -910 
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and, hence, 


2n 2—cos g-sin@ 2x 
n= { ap dz = | d= 2a. 
0 i-cos g-sin @ 0 


119. Find the flux of the vector 
a= yi+ zj — evk 


through the outer side of the lateral surface of the cylinder 
zx? + y* = 4 bounded by the planes z = 0 and z+ y+ 
+2=4, 
120. Find the flux of the vector 

a = zi — zyj + 2k 
through the outer side of the cylindrical surface z? + 2? = 


= R? bounded by the planes y = 1 and 2+ y= 4. 
121. Find the flux of the vector 


a= i — yj + rk 


through the outer side of the cylindrical surface z? +- y* = 
= 9 bounded by the sphere x? + y*? + 2? = 25 
122. Find the flux of the vector field 


a = zi — yj — ry2k 


through the outer side of the lateral surface of the cylin- 
der z* + y* =1 bounded by the plane z = 0 and the 
hyperbolic paraboloid z = 2? — y’. 

123. Find the flux of the vector field 


a = (zy — y*) it (22 — z* + zy) j + zk 


through the outer side of the lateral surface of the cylin- 
der 2? + y?=1 bounded by the elliptic cone 2? = 
= 2/2 + y’. 

Case (2). Suppose the surface S is a part of the sphere 
xz? + y? + 2? = R*® bounded by conical surfaces whose 
equations in spherical coordinates have the form 0 = 
= 7, (9), 9 =f. (g) and by the half-planes 9 = q,, 
P = Po. 

For the points of the given sphere, set 


z= RceosgsinOo, y=Rsingsin§, z= Roos8, 
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where 9, < 9 XS Q, and 0, < 89 < B,. Then for the 
element of area dS we obtain (Fig. 26) 

dS = R?* sin 0 d6 dg. 


In this case, the flux of the vector field a through the 
outer part S of the sphere is found from the formula 


UPN 6, 
T= R? j dg | (a, n°) sin 0d9, (14) 
e: 0, 


where 
o_ _gtad (22+ y3+22—R%) zit yjtzk 
“| grad (z3++ y?+ 2?— R*) | R , 


Example 13. Find the flux of the vector 


Fig. 26 


a= (t# — 2y + 1)i+ Q24+ y — 32)3 + By + 2)k 


through a part of the surface of the sphere 2? + y*? + 2? = 
= 1 located in the first octant into a region where 
e+t+yt+2> 1. 


6* 
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Solution. Here we have 
R=1, 9, =0, gp = 7/2, 6, = 0, 0, = n/2, 
n°? = zi + yj + zk, (a, n) = 2+ y+ 24+ 2, 
On the sphere z? + y? + z* = 1 we introduce coordinates 
@ and 6 so that 
z=cosgsin®, y=singsinO, z= cos 0. 
We then have 
(a, n°) = 1 + cos m sin 0 
and, using (14), we obtain 


m/2 n/2 
fl = j dp \ (1+ cos 9 sin 8) sin 6d0 
0 
nm/2 1/2 m/2 m/2 
~ Perf snoan+ Posed] suroao=3 
=  \ sin + \ cosq@dg | sin 6dd=—> a2. 
0 0 0 0 
124. Find the flux of the vector field 
a= Fi — yj + zk 
through the outer side of that part of the sphere z* + y? + 
+ 2? =1 cut out by the conical surface 2? = z? + y?, 
z> Va + y*. 
125. Find the flux of the vector field 
a = ya + zzj + zyk 


through the outer side of that part of the sphere 2? + y? + 
+ z? = R® located in the first octant. 
126. Find the flux of the radius vector 
= zi + yj + zk 

through the outer side of that part of the sphere x? + y* + 
+ 2? = 2 that is bounded by the planes z = 0 andz = y. 
127. Find the flux of the vector 

a = 22i + yzj + 2k 
through the outer part of the sphere z* + y2? + 22 = 9 
cut off by the plane z = 2 (z > 2). 
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Sec. 12. The flux of a vector 
through a closed surface. 
The Gauss-Ostrogradsky theorem 


Theorem. If in some region G of space the coordinates 
of a vector 


a=P(z,y,z)i+ Q(z, y, 2j+ Rx, y, 2k 


are continuous and have continuous partial derivatives 
OP/dz, dQ/dy, AR/dz, then the flux of the vector a through 
any closed piecewise smooth surface = located in G is equal 
to the triple integral of OP/dx + dQ/dy + AR/dz over 
the region V bounded by the surface >: 


n= (, n°) do = ny) ($+2+2)H 
z= 


(the Gauss-Ostrogradsky formula). 

The normal n to the surface = is taken to be the outer 
normal. 

Example 1. Compute the flux of the vector 


a= ri+ yj + 2k 
through the closed surface 
P+yt+ a= Rt, z=0 (2 >0). 
Solution. By formula (1), 
I = (2x + 2y + 22) dv. (2) 
i 


The integral (2) is conveniently computed in the spherical 
coordinates r, 8, o. We have 


z=rsinOcosg, y=rsin@6@sing, z=rcos0 
and the element of volume is 


dv = r* sin 8 dr d8 dp 
so that 


Il=2 ee) (rsinOcosg+rsin6Osin@ 
Vv 


+rcos 8) r?sin 8 drd@dqp 
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on n/2 R 
= 2 dp j sin 6 (sin 6 cos g + sin 9 sin @ + cos 8) 40 {r°dr 
0 0 0 


2n 
= 2K { a o| cos@ sin 0.d0 = 78" ; 
0 0 
Example 2. Compute the flux of the vector 
a = 4zi — yj + zk 


through the surface of a torus. 
Solution. Taking advantage of the Gauss-Ostrogradsky 
theorem we find that the desired flux is 


n= Poo n°) do = ff (4+2+23)av—w, 


where V is the volume of the torus. To compute the 
volume V, let us take advantage of the Guldin theorem 


Fig. 27 


on the volume of a solid of revolution, by virtue of which 
the volume is equal to the product of the area of the 
rotating figure into the path covered by the centre of 
mass of the figure during the rotation. 

Let A, and R, be the inner and outer radii of the torus 
(Fig. 27). The area S of the circle, which during rotation 
forms the torus, is equal to 


San(Bhy’ 


Sec. 12] The Gauss-Ostrogradsky theorem 87 


The path length described by the centre of mass (the 
centre of that circle) is the circumference 1 of a circle 
of radius (R, + R,)/2, that is, 


L=2n ERs _ (K+ R,). 


Thus, the volume V of the torus is equal to 


= 2 
Vina (F85*8)'n (Re + Ry) = (Ro— Ry)? (Rot A). 
The desired flux is 
= 1? (Ry — Ry)? (R.+ Ry). 


Example 3. Using the Gauss-Ostrogradsky theorem, 
compute the flux of the vector field 


.  «t 1+y? 
a= (774 )i i+ 2zarctan y-j— CFP SITY 5 
through the outer side of that part of the surface z = 
= 1— 2? — y* located above the zy-plane. 
Solution. In order to be able to apply the Gauss- 
Ostrogradsky theorem, close the given surface from 


below with a portion of the zy-plane that is bounded by 
the circle 


P+y=1, } 
z=0. 


Let V be the volume of the resulting solid bounded by 
a closed piecewise smooth surface o consisting of a part 
o, of the paraboloid of revolution z = 1 — x? — y? and 
a part o, of the plane z = 0 (Fig. 28). 

The flux of the given vector through the surface o is, 
by the Gauss-Ostrogradsky theorem, equal to 


= Yea nao= | SS (FE 2 +2) av. 


We find a sum 


P , a , OR 2z (1+) 
at+Et fee Si: 


tee ity i+y3 
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Consequently, the flux is 
m=Q (a, n°) do = 0. 
fe] 


By virtue of the additivity of the flux we have 
Il= i) (a, n°) do+ j j (a, n°) do=0, 


Fig. 28 


From this the desired flux is 


IT, = (a, n°?)do = — (a, n°) do. 

J) JS 
The flux I], of the vector a through the circle 22 + y?>< 
<1, 2 =0 is equal to 


Tl, = (a, n°) do. 
i 
Since on the plane z = 0, we have 
as Ae 


ay ri 


i+ 2zarctany-j—k, n°=—k, 
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and hence (a, n°) = 1, it follows that the flux IT, through 
the circle o, is equal to the area of the circle o,: 


N= f do=n. 


o; 


The desired flux I], = —II, = —n. 

By appropriately closing the given unclosed surfaces 
and making use of the Gauss-Ostrogradsky theorem, com- 
pute the fluxes of the vector fields through the indicated 
surfaces (the normal to the closed surface is the outer 
norma)). 
128.a=(1—2z)it+yjt2k; 8S: 2+ y=2 
O<z< 4). 

1299. a=2i+aj+yk; S: 2 +y=4—2 D> 0). 
130. a= (y+ 2)i— yj+ yk; S: r+2=y' 
(O0<y€ 1). 


Sec. 13. The divergence of a vector field. 
Solenoidal fields 


The notion of the flux of a vector through a closed 
surface leads to the concept of the divergence of a field. 
This concept yields a certain quantitative characteristic 
of a field at each point in the field. 

Let M be a point of the field under study. Surround it 
by a surface 2 of arbitrary shape, for instance, a sphere 
of sufficiently small radius. Let the region bounded by 
the surface = be (V) and its volume V. We consider 


the ratio 
vy) (a, n°) do 
Lt 


+—. (1) 


Definition. If the ratio (1) has a finite limit when the 
region (V) contracts to the point M, then this limit is 
termed the divergence of the vector field (the divergence 
of the vector a) at the point M and is designated as 
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Ny (a, n°) do 


diva(M)= lin =. (2) 
(V)-M 


div a(M). We thus have 


Formula (2) yields an invariant definition of diver- 
gence. This definition means that the divergence of the 
field a at the point M constitutes the volume density of 
the flux of the vector a at that point. 

The points M of the vector field a (M) at which div 
a>0O are termed sources, while the points at which 
div a < 0 are termed sinks of the vector field. 

The divergence of a vector field is a scalar function of 
the points of the field. 

If the coordinates of the vector 


a (M) = P (z,y,z)i+ Q(z, y,2)) + FR (z,y,2)k 


have continuous partial derivatives 0P/dz, dQ/dy, AR/dz 
in the neighbourhood of the point M (z, y, z), then, using 
the invariant definition of divergence, we find from the 
Gauss-Ostrogradsky theorem that 


. P 
diva= T+ 24%. (3) 


All the quantities in (3) are considered at the same point 
M (a, y, 2). 

Using (3) for divergence, we can write the Gauss- 
Ostrogradsky theorem (see Sec. 12) in vector form: 


Wy (a, n°) do = SJ div adv (4) 


=z 


Example 1. Using the invariant definition, compute 
the divergence of the vector a = zi at the point O (0, 0, 0) 
taking for the surfaces o surrounding the point O the 
spheres o, of radius e centred at that point. 
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Solution. By the definition of divergence, we have at 


the given point 
) (a, n°) do 


diva(0)= lim — 


(0,)-+0 Ye ; 


where v, is the volume of a ball bounded by a sphere 


Oe, OF 
§ (a, n°) do 


div a(0)=lim — 
e-0 


But since the volume of the ball is v, = 4ne/3, it fol- 


lows that 
) (a, n°) do 


div a(0) =lim 2 
e-0 nes 


CS 


Let us compute the flux Ha, n°) do of the given vector 


¢, . 
through the sphere o,. The unit vector of the normal n?® 
to the sphere o, is directed along a radius of the sphere, 
and so we can put 
“—S r 
be 


where r° is the unit vector of the radius vector r = 
= zi+ yj + zk, or 


n°= r? 


iF 
=~ 


zit+yjt+zk 


€ 


n° — 
The desired flux is 


(a, n°) do = i) ca do. 


% 
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Passing to coordinates on the sphere o,, 
z=ecosqgsinO, y=esingsin®@, z= ecos@, 
we get 


3 2 n3 Qe; 
G ca, op aoa [ { steatgsinteeruneceas 
Ge fo 


e 
2x n 4 
= e3 j cos? p dp j sin’ @ dg=z ne’, 
0 0 
Consequently 


Bone 
a me’ 
div a(0)=lim 


2-0 4 a¢3 
3 


={. 


Example 2. Compute div r. 
Solution. We have r = zi + yj + 2k, so that P = x, 
Q@ =y, R =z and, hence, by formula (3), 


atte ae Ore es — 
divr= Bait Oy + = = 3. 
Example 3. Compute div (u- a), where u (M) is a scalar 
function and a(M) = P(z, y, z)it+Q(z, y, z)jt 
+ R (z, y, z)k is a vector function. 
Solution. Using formula (3), we get 
div (ua) 
_9(uP) , 9(uQ) , 9(uR)__ so ou 9Q ou 
ee ap oe oe Be Op ee oe 
oP 6Q OR du 
pian te 
Ou du 


+05 +R a =udiv a+ (a, grad wu). 


OR Ou 
tha, Ta a 


Thus 
div (ua) = u div a + (a, grad u). (5) 


Example 4. Find the divergence of the vector 


a= @(r) r= 2 Tr, 
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where r = |r | is the distance from the coordinate origin 
tu the variable point M (gz, y, 2). 
Solution. Using formula (5), we obtain 


diva = divr+ (r, grad 27) ) ; 
Furthermore, 
divr=3, grad at) = ( ete) ) eee A290 o. 
And so 
div ax 2 34 (z =e) r?, r) 


=3 20) ee =e) 2 ot) +9’ (r). 


131. For what function »p (r) will we have div (r)r = 
= 2p (r)? 
132. Find div (rr). 
133. Find the divergence of the vector field 
a = [e, rl, 

where ¢ is a constant vector. 
134. Find 

div (r [w, r]), 


where w is a constant vector. 
135. For what function + (z) will the divergence of the 
field 
a= z2i + yj + »p (z)k 
be equal to 2? 
136. Find the flux of the radius vector r through the 
surface of a sphere. 
137. The electrostatic field of a point charge q is 


227 8 
Compute div E. 
138. Show that 


+ (r, n°) do=V, 
z= 
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where V is the volume bounded by a closed surface 2. 
139. Prove that if 2 is a closed piecewise smooth surface 
and c is a nonzero constant vector, then 


i) cos (u, c)do =0, 

z 
where n is a vector normal to the surface 2. 
140. Prove the formula 


yp (pa, n°) do= j J j (p div a+ (a, grad g)) du, 


where ~ = 9 (2, y, 2) and 2 is a surface bounding the 
volume V. Establish the conditions under which this 
formula is applicable. 
141. Prove that if the function u (z, y, z) satisfies the 
Laplace equation 

03u Ou 04u 

oat + art oar = 


i) # do=0, 
= 


where du/dn is the derivative with respect to the direction 
of the outer normal to the piecewise smooth closed sur- 
face 2. 

142. Prove that if the function u (z, y, 2) is a second- 
degree polynomial and 2 is a piecewise smooth closed 
surface, then the integral 


then 


is proportional to the volume bounded by the surface 2. 

Find the flux of a vector field through the indicated 
closed surfaces: (1) directly, (2) via the Gauss-Ostrogradsky 
theorem. 


zer*+y’, 


143. a=zi+zk; S: | eer 
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144. 


145. 


146. 


147. 


148. 


149. 


150. 
151. 


152. 


- a=zi—2yj—zk; S: { 


22 = 22+ y?, 
z=4H, H>0. 
2=6—2?— y?, 
2=g+y’, z>0. 
z+ 22= y?, 
= z+ y?, 
z=2?+y?. 
a+ y= 4, 
a=2zri—(z—1)k; S: ee z=1. 
+ y+ = 4, 
32 = 224 y? (z>=T). 


a=2zi+ 2yj—zk; S: { 
a= zi— zj; s:{ 
a=yzi—zj—yk; S: { 


a=zi+2yj—zk; S: { 


a=2zri—yj+ zk; 8S: { 


a=yzi+ 2Qyj—zk; S: 22+y?+22=4, 
a= (z+ z)i+(y+z)j+(z+2)k; 
S: lace 
* [z=y, z>0. 
9—z=7?4y, 
a=3zi—yj—zk; S: 7 z=0, y=0, z=0 (first 
octant). 


r+y+z=1, 


. a=(y—z)i+(z—y)j+(z—z)k; 8S: {epee 


z=0, z=0. 
1—z=27+ 7’, 
z=0. 


Solenoidal fields 


Definition. If at all points M of a certain region G 
the divergence of a vector field (specified in G) is zero, 


div a (M) = 0, 


then we say that the field is solenoidal in that region. 
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Thus, a solenoidal field is, by definition, without 
sources and sinks. 
From the Gauss-Ostrogradsky theorem it follows that 
in a solenoidal field the flux of a vector a = a (M) 
through any closed surface o 
n lying in the field is zero: 


(a, n°) do =0, 
: 


In a solenoidal field G, the 
vector lines cannot begin or 
end. They can either be closed 


L curves or have ends on the 
boundary of the field. 
Fig. 29 The equation 
div a (M) = 0 


is encountered in hydrodynamics and is termed the con- 
tinuity equation of an incompressible fluid. 

In this case, the amount of fluid exiting through some 
closed surface o is always ‘equal to the amount entering, 
and the total flux is zero. 

Which of the following vector fields is a solenoidal 
field? 

155. as2e(2?—y/)ity@—z)jtzy’?—2)k. 
156. a= yi — (2? + y*)j tz (By? + 1)k. 

157. a = (1 + Qzy) i — y®ej + (2’y — Qzy + 1)k. 
158. Show that the field of the vector 


E=47P (r=VP+yP+2) 


is solenoidal throughout the region that does not contain 
the coordinate origin O (0, 0, 

159. Under what condition will the vector field a = 
= p(r)r be solenoidal? 

Suppose we have a field a (M), which is not necessarily 
solenoidal. In the field, consider a closed oriented con- 
tour L. The surface 2 containing the curve L as its edge 
will be called a surface spanned by the contour L. Let us 
agree to orient the normal n to the surface 2 so that the 
chosen circuit about the contour LZ will be seen from the 
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end of the normal as being counterclockwise (Fig. 29). 
160. Show that in a solenoidal field the flux of a vector 
a (M) does not depend on the type of surface = spanned 
by the given contour Z and depends solely on the contour 
itself. 


Sec. 14. A line integral in a vector field. 
The circulation of a vector field 


Suppose we have a continuous vector field a = a (M) 
and a piecewise smooth curve LZ on which a positive 
direction has been chosen (in other words, L is an oriented 
curve). 

Definition 1. The line integral from the vector a = 
= a(M) along the oriented curve L is termed a line 
integral of the first kind (an integral over the arc length 
of a curve) of a scalar product (a, 1°): 


j (a, t°) ds, 

L 
where t® = +° (M) is the unit vector of the vector tangent 
to the curve L whose orientation coincides with that 
of L; ds is the differential of arc length s of L. 

Ifr = r (M) is the radius vector of an arbitrary point M 

of the curve L, then the line integral in the field a (M) 
may be written thus: 


J (ae) as—= { (a, dr), (1) 
L L 


If a rectangular- coordinate zyz-system is introduced 
in the vector field, then r = zi + yj + zk, 
a=P(z,y,z)it+Q@y2z2i+ RG, y, zk 


and the line integral (1) can be expressed in terms of a 
line integral of the second kind: 


j (a, dr) = f P (z,y, 2) dz+Q(z, y, 2) dy +R (z, y, 2) dz. 
L L 
When a = a(M) is a force field, the line integral (1) 


yields the magnitude of the work of that field over the 
curve L. 


7-910 
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Properties of a line integral 
(a) Linearity: 


J (ay t pas dr) =A | (ay, dr) + \ (a, dr), 
L L 


where A and p are constants. 
(b) Additivity: 
(a, dr) = { (a,dr) + { (a, dr), 
Ly +L? Ly Le 
(c) The integral reverses sign with a change in the 
orientation of L: 


J (a, dr) = — j (a, dr), 


where A is the initial and B is the termina] point of the 
curve L, 


Calculating a line integral in a vector field 
Let the curve L be specified by the parametric equations 
z= 9 (t), y=» (t), z= x(t) to <t< ty). 
Here, at the initial point A of Z the parameter ¢ assumes 
the value ¢ = ¢, and at the terminal point B of L it 
assumes the value ¢t = ft, (the direction on L corresponds 
to increasing values of the parameter ¢ from ty to t,); 


the functions @ (t),  (¢), x (¢) have continuous deriva- 
tives on the interval [t,, ¢,). Then 


ty 
fardn= | @dr=f (Pie, 00. xOle' 
L AB to 
+Q [0 (t), p(t), x ()) Y (+ Rig (1), x (1) x ()} at. 
If the curve L is given by a system of equations 
y=*(z), 2=x(z%),¢@<2< 6, then 
b 
(a, dr) = J {P la, (2) (21+ Qe, (2), 4 (2)} W' (2) 
AB a 


+ Riz, p(z), x (2) x" (z)} az. 
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Similar formulas may also be written for cases where 
the line is specified by one of the following systems 
of equations: 


z=o(y), z=xty Woycn) 


zr=@(z), y=vp(z) (<2< 4). 


Example 1. Find the line integral of the vector a= 
=r/|r |, where r is a radius vector on the line segment 
from point A (ra) to point B (rz). 

Solution. The desired line integral is 


f (adr) = w (1) 
AB AB 


or 


From 
d (r, r) = (dr, r) + (r, dr) = 2 (r, dr) 
we find 
(t, dr) => a(n) = d (rl) =4-2 |r] d |e] = |r] |r], 
whence 


(r,dr) 
aa ie =d|rl. (2) 


Substituting (2) into the integral (1), we get 


"p 


J (ade) = fairl = J dir) = isl — Ira. 
AB AB ra 
Note that 
| dr | 4dr |. 


Find the line integral over the line segment bounded 
by the points A (r,) and B (r,) for the following vector 
fields: 

161. a =r. 


162. a=—_ 
r 


163. a= ia , f° is the unit vector. 


164. Compute the line integral over the straight line 
passing through the points O (0, 0, 0) and M, (1, 1, 1) 


7* 


100 Vector Fields (Ch. III 


from point O to point M, if a = [b, r] , where b is a con- 
stant vector. 
165. Prove that 


{ (grad u, dr) =u (B)—u(A). 
4B 


Example 2. Find the line integral of the vector 
a=2z+ aj + yk 
over an arc L of the heli- 
cal curve 
z=Reost, y=Rsint, 
t 
t= On 
from point A, the point of in- 
tersection of the curve with 
the plane z = 0, to point B, 
the point of intersection with 
the plane z = 1 (Fig. 30). 
Solution. Here, the line in- 
tegral is of the form 


| (ar) 
L 
= j zdx+xzdy+ydz. 
L 


The helical curve is located 
Fig. 30 on a circular cylinder 2z* + 
+ y? = R?. At point A we 

have t) = 0, at point B we have t, = 2x. Since 


dz= —Rsintdt, dy=Rcostdt, a=—, 


it follows that the integral is equal to 
on 
j (a, dr) = j ( = Rsint + Rcos*t-+ + sint) dt 
L 0 
on 2n 
= R2 j cos?t dt — j tsintdt=nR2+R 
0 


2n 
0 
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because 
2n 
j cos*tdt= 7; 
0 
an 
j tsintdt= —2n. 
0 


Example 3. Find the line integral of the vector (see 


example 2) feat lege 
a=z y 


over the straight line AB (see Fig. 30) in the direction 
from point A to point B. 

Solution. Since the straight line AB (the generatrix 
of the cylinder z? + y* = R?) is located in the zz-plane 
and passes through the point A (R, 0, 0), it follows that 
y=0, z = R, dz = 0, and for the radius vector r of 
the points of AB we will haver = Ri + zk, dr = dz-k. 
Therefore the scalar product 


(a, dr) =zdr+ardyt+yd 


on AB will be zero. 
Hence, the desired line integral 


| (a, dr) = 5 (a, dr) 


on AB will be zero. 

From examples 2 and 3 it follows that in the general 
case a line integral in a vector field depends not only on 
the initial and terminal points of the path of integration 
but also on the shape of the path. 

Example 4. Compute the work done by the force field 

F=yitat+@t+yt+z2)k 
along the segment AB of a straight line passing through 
the points M, (2, 3, 4) and M,(3, 4, 5). 

Solution. The work done by the force field will be 

equal to the line integral along M,M,: 


A= Lie j ydz+<2dy+(zr4-y+ 2) dz, 
MM, 


MiMg 
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Let us find the canonical equations of the straight line 
M,M,. We have 


1 eT 
whence 
y=2+1, 
aes 


Here, z varies from 2 to 3 (since the abscissa of M, is 
equal to 2 and the abscissa of M, is 3). The desired work is 
3 


A= LORere er trate 


3 
= \ (524+ 4)dr= 3. 

2 
166. In the plane vector field 
yi-zij 
Vite 
compute the line integral along the semicircle 

z=Reost, y=Rsint QOStgna). 
. In the plane vector field 
a=(?+ y)i+@—y)j 
compute the line integral over the curve y = | z | from 
the point (—1, 1) to the point (2, 2). 
168. In the plane vector field 
a = (2° — dry) i + (y? — 2zy) j 

compute the line integral: 

(a) along the parabola y = 2? from the point (—1, 1) 
to the point (1, 1); 

(b) along a segment of the straight line joining the points 
(—1, 1) and (4, 1). 
169. Compute the work of the force field F = 2ryi + 2°j 


along an arc of the circle z* + y* = 1 from the point 
(1, 0) to the point (0, 1) counterclockwise. 
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170. Compute the work of the force field 

F= (2 + dzy)it+t(@t+y’)j 
along the parabola y = 2? from the point (0, 0) to the 
point (4, 1). 
171. In the vector field 

_ zi+yj+7zk 

compute the line integral along a segment of the straight 
line from the point (1, 1, 1) to the point (4, 4, 4). 
172. In the vector field 

a = (y? — 2%) i + Qyzj — 2k 
compute the line integral over the line 


r=t 
inet (0<t<1) 
z= 


in the direction of increasing values of the parameter t. 
173. In the vector field 

a= yi+ a+ zk 
compute the line integral along a turn of the helical curve 


z=acost, 
y=asint, (0<t<2n) 
z=bdt 
in the direction of increasing values of the parameter ¢. 
174. In the vector field 

a=ri+ yt 2k 
compute the line integral in the direction from the point 
(0, 0, 0) to the point (1, 1, 1) along the line segment 
between the two points. 


Computing the circulation of a vector field 


Definition 2. The circulation C of a vector field a = 
= a (M) is a line integral taken around a closed oriented 
curve L Thus, by definition, 


C=§ (a, dr), 
L 
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where the symbol § denotes the integral around the 


> 


closed curve L. 
If the vector field a = a(M) is given in coordinate 
form, 


a = P (z, y, zji+ Q(z, y, z)j a R (z, y, z) k, 
then the circulation of the vector field is 


C=§ Pdr+Qdy+ Raz. 
L 


For the positive direction of traverse around the closed 
curve L we take the direction in which the region bounded 
by the curve is on the left. 

Example 5. Compute the circulation of the vector 


field 
a = —yi + 25 


around the ellipse L: 27/a? + y?/b? = 1. 
Solution. By the definition of circulation we have 


C=§ (a, dr) = —ydz+ 2 dy. (3) 
L L 
The parametric equations of this ellipse are of the form 
z=acost, 
vam (0<t<2a), (4) 
whence 
dz = —asint dt, dy = bcost dt. (5) 


Substituting (4) and (5) into (3), we get 
on 
C =ab j (62 sin‘ ¢ + a? cos* t)dt= = stab (a? + b?) 
0 


since 
on 2n 


) satel j (1 —cos 22)? dt 
0 


=f (1—2cos 2¢ + tts0e% ) ae 
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73 { 163 3 
J aes Daal Cee oe sa Tl. 
Similarly, we find that 

2n 

j cos tdt= > nt. 

0 

Example 6. Calculate the circulation of the vector 
field 
a = ye™i + ze™j + cyzk 


around the curve L obtained by cutting the cone z? + y? = 
= (2 — 1)? with the coordinate planes (Fig. 31). 


Fig. 31. 


Solution. The curve L consists of two parts BC and CA 

located in the coordinate yz- and zz-planes respectively, 
~ rv+y=1 

and the arc AB of the circle z=Of° 


the circulation of the given vector field will be 


Therefore 


C= § (adr) = j (@.dr) + { (a, dr) + | (a, dr). 
L Be CA “4B 
(1) On the line segment BC we have 


z=0,dr=0; z=1—y, dz=—dy; 1 >y>0. 
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Consequently, 
j (a, dr) = j ydr=0. 
BC BC 
(2) On the line segment CA we have 
y=0, dy=0; z=1—2,d2=>-—dz; 0g r<i1 
and so 


j (a, dr) = { zdy=0. 


CA CA 
-_~ z+y?=1, 
(3) On the arc AB of the circle 0 we have 
z=0, dz=0, which means that 
j (a, dr) = j ee" (ydx+a2dy) = { et¥d(2y) 
AB AB 4B 
= xv sy[BO Ogg 
= j ae oral RT Cr a 


AB 
The desired circulation of the vector field is zero. 
Example 7. Compute the circulation of the vector 
field 
a = ryi + y2zj + rzk 
if 
is 2 +y= 1, 


t+y+z=1. 
Solution. We have 


C=O (a, dr) = j zy dx + yzdy+a2zdz. 
L L 


The curve L is an ellipse obtained by cutting the cylinder 
z7+ y? =1 by the plane r+y+z= 41. Let us find 
the parametric equations of the curve. The projection 
of any point of the curve on the zy-plane lies on the 
circle 27+ y3=1. From this we obtain z = cost, 
y = sin t. But the ellipse lies in the planex + y +2 = 1, 
whence z = 1 —z—y or 2 = 1 — cost — sint. Thus 
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the parametric equations of the curve L are: 


zxz=cost, 
y=sint, (0<t<2n). 
2=1—cost—sint 


From this we find 
dx = —sint dt, dy = cost dt, dz = (sin t — cost)dt 
and so the circulation is 
on 
C= j [—costsin?é¢+ sint(1—cost —sin2Z) cost 
0 
+cost(1—cos¢—sin ¢) (sint —cos ¢)] dt 


2n 
= j (—3 sin?écost+ sin 2¢— cos? ¢ sin t—cos? ¢ 4+- cos? t) dt 
0 


on 
= — | cos?¢dt= — n. 
0 


Compute the circulation C of a vector a around the 
given curve L: 
175. a = (zz + y)it (yz—2z)j—@+y)k; 
’+y=t1, 
L: { +y 


2=3. 
176. a = yi + 275 + 27k; 


L z+ y?+4 22= R2, 
{ P+y= Rr 


177. a = (2x + 2)i + (2y — 2) j + zyzk. L is the line of 
intersection of the paraboloid of revolution z? + y? = 
= 1-—2z with the coordinate planes. 

178. Show that if in a vector field the circulation of a vec- 
tor around any closed circuit is equal to zero, then there 
can be no closed vector lines in such a field. 


(z>0). 
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Sec. 15. The curl (rotation) 
of a vector field 


Suppose we have the vector field 

a(M) = P(z, y, 21+ Q(z, y, 2jI+ RG, y, zk. 
We will assume that the coordinates P,Q, R of the vector 
a (M) are continuous and have continuous partial deriv- 
atives of the first order with respect to all its arguments. 

Definition 1. The curl (or rotation) of a vector a (M/) 
is a vector (symbolized: curl a (M) or rot a (M)) defined 
by the equation 


OR\. 

curl a= (5 —) i+ (4-4) 

30 

+ (35 or =) k 

or, in easy-to-remember symbolic form, 
ij k 
6.6068 

curl a= a oy a I: (2) 

P QR 


This determinant is ordinarily expanded in terms of 
elements of the first row, in which case the operations of 
multiplication of elements of the second row by elements 
of the third row are regarded as operations of differen- 
tiation; for example, 


Definition 2. If in some region G we have curl a = 0, 
then the field of the vector a in G is said to be irrotational. 
Example 1. Find the curl of the vector 


a=(r+2)i+ yt2j+ @+2)k. 
Solution. Using formula (2), we obtain 
i j k 
0 te] 0 
Oz Oy 
r+zy+2 22+2 


curl a= 
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Expanding the determinant in terms of elements of the 
first row and .regarding the operation of multiplication, 
say 0/dy by x? + z, as an operation of partial differen- 
tiation, we obtain 

curl a = —i— (2z — 1) j. 


Example 2. Find the curl of H, the intensity vector of 
a magnetic field under the conditions of example 3 of 
Sec. 10. 

Solution. The magnetic-field intensity vector H is 


2 
H=—-y Ihr] 
or 
ijk 
2 2h. gs. 2h 
H=—r OOo] = — or it oF th 
zyz 
where p? = z* + y?, whence, by (2), 
i j k 
0 G) C] 0 2lz 
calH=| 3 a wl=le (age) 
aly 2iz 0 


eye ay 

0 2ly _ z?4y2—2z3 x24 y?— 2y2 = 

og (sree) k= 2 [Seca + “perigee |k=0 
(z?+y? # 0). 


Thus, curl H = 0 everywhere except the z-axis, at the 
points of which the last formulas are meaningless (the 
denominator vanishes), that is, the field of the vector H 
is irrotational everywhere outside the points of the z-axis. 

Find the curl of the following vectors: 

179. a = ane ge Uy + 2) i+ (24 2’)k. 
180. ‘a = a + yj 
181. a = . (—y*i + a 
182. Show that if the coordinates of the vector a (M) 
have continuous partial derivatives of the second order, 
then 

div cur] a = 0, 
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that is, thefield of the vector curl a (M) is a solenoidal 
field. 
183. Show that 

(a) curl (a + b) = curl a + curl b, 

(b) curl (Aa) = A curl a. 
where A is a numerical constant. 
184. Show that if u = u(M) is a scalar function and 
a = a(M) is a vector, then 


curl (ua) = u curl a + [grad u, al. 


185. Show that if a and b are constant vectors, and r 
is the radius vector of the point M (z, y, z), then 


curl (r, a) b = [a, bl. 
186. Show that 


curl (ra) = + (r, al, 


where a is a_ constant vector and r= |r| = 
=V2+/+ 2. 
187. Show that curl (f(r) a) = a [r, al, where f (r) 


is an arbitrary differentiable function of its argument and 
a is a constant vector. 

188. Show that the vector field a = f (r) r is irrotational 
that is, curl a = 0. 

189. Show that 


div [a, b] = (b, curl a) — (a, curl b). 


190. Show that the curl of a field of linear velocities v of 
a rotating solid is a constant vector in the direction paral- 
lel to the axis of rotation, the modulus of which is equal 
to twice the angular velocity of rotation: curl v = 2. 
191. Determine the angular velocity w of rotation of 
a solid about a fixed axis passing through some point of 
the solid if its linear velocity is 


v = dai + yj + ok. 


192. Show that the field of the curl of the vector a (M) 
is free of sources and sinks. 
193. What must the function f (z, z) be so that the curl 
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of the vector field 
a = yzi + f(z, 2) j + zyk 


is coincident with the vector k — i? 


Sec. 16. Stokes’ theorem 


Suppose the coordinates of a vector 
a(M) = P(z,y, z)it+ Q(z, y, 2j+ Re, y, zk 


are continuous and have continuous partial derivatives. 

Theorem. The circulation of 
the vector a aroundaclosed con- 
tour L is equal to the flux of the 
curl of the vector through any 
surface X spanning the con- 
tour L: 


§ (a, dr) = I (cur) a, n°) do. 
(1) 


It is assumed that the orien- 
tation of the normal n° to the 
surface 2 is matched with the 
orientation of the contour L @ 
so that, when viewed from the Fig. 32 
end of the normal, thecontour 
is traversed in the chosen direction counterclockwise. 

Example 1. Compute the circulation of the vector 


a= yi + vi — zk 


+ 
around the contour L: . 
(2) via the Stokes theorem. 

Solution. (1) The contour L is a circle of radius R = 2 
lying in the plane z = 3 (Fig. 32). We choose the orienta- 
tion as shown in the drawing. The parametric equations of 


is (1) directly and 


442 Vector Fields (Ch. II! 


the curve L are 


z=2cost, 
y=2sint, (0<t< 2n), 
z=3 


so that 
dx = —2sintdt, dy=2costdt, dz=0. 
For the circulation of the vector a we have 
on 
C= j {2sin¢t(—2sin t) +4 cos?¢2 cost—3-0] dt = —4n. 
0 
(2) To compute the circulation via the Stokes theorem, 
choose some surface 2 spanning the contour L. For 2 
it is natural to take a circle having. L as its boundary. 


According to the chosen orientation of the contour the 
normal n° to the circle has to be taken equal tok: n? = k. 


Then 
ij k 
6.6C«@#é# 
cur] a = de. oy ae = (27—1)k 
y 2 —2 


Therefore, by the Stokes theorem, 


C= \I (curl a, n°) do = \J (22—1) do 
on 2 
2 


2 
= j ae | (2p cos ~—1) pdp= —2n — —An. 
0 

194. Show that the flux of the rotor through an open 
surface spanning a given contour does not depend on the 

shape of the surface. 
Find the circulations of the vectors around the indi- 
cated contours (1) directly and (2) via the Stokes theorem. 


24 y2— 4, 
195. a=zitajtyk; L: {* ie 
z=0. 
P+ y+ = 4, 


196. a=yi—zj+zk; L: { e+y2=z? (2>0). 
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197. a = 2x2i — yj + zk around a contour formed by 
intersection of the plane z + y + 2z = 2 with the coor- 
dinate planes. 


198. a=yi—zj+(c+y)k; L: aa 
Pry? +2 = 16, 
z=0,y=0,z=0. 
r=yr+ 2, 
z=9. 
e+y=9, 
3y+4z =5. 
e+ y+ z=, 
rT=Z. 


199. a=2zi;_ L: 
200. a=zy2i+ r22j+ 22yk; L: 
201. a=y*i+ 22j;  L: { 


202. a=yi—azj+z2k; L: | 


203. Given a vector field of velocities v of the points 
of a solid rotating with a constant angular velocity 
about the z-axis. Calculate the circulation of the field 
around the circle 


z=acost, 
L: Sy=asint, (0<t<2n), 
z=0 


directly and via the Stokes theorem. 

From the Stokes theorem we find that the projection 
of the vector curl a on any direction n is independent of 
the choice of the system of coordinates and is equal to 
the surface density of the circulation of the vector a 
around the contour of the area perpendicular to that 
direction: 


$ (a, dr) 


Pra cur! ala = (curla, n°)| = lim + 5 
M (2)-+M 


(2) 


Here, (2) is a plane area perpendicular to the vector n; 
S is the area of (2); D is the contour of the area and is 
oriented so that the traverse of the contour is counter- 
clockwise as seen from the end of.the vector n; the nota- 


8-910 
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tion (2) — M means that the area (2) contracts to the 
point M at which we consider the vector curl a, and the 
direction of the normal n to that area is always the same. 

Example 2. Compute the density of the circulation of 
the vector a = yi around the circle 


xr=acost, 
L: 4y=asint, (O<t<2n), 
z=0 
at the centre of the circle, M (0, 0, 0), in the positive 
direction of the z-axis. 
Solution. Here, (2) is a circle of radius a with centre 
at M, so that S = na?. 
The desired density of the circulation is 


; 4 3 1 
peng ti oe 


ox 


= lint sat | (— a?) sint?¢dt= —1, 
On the other hand, 
ij k 
Q 0 
curl a= ae oy = —k 
y 0 0 


and, hence, 
(curl a, n°) |x, = (—k, k) = —1, 


which, by virtue of (2), confirms the correctness of the 
result. 

204. Compute the density of the circulation of the vector 
a=a-+a2j+ yk around the circle L: {y = acost, 
z=asint, z=0, (0 <t< 2n)} at the centre of the 
circle, C (0, 0, 0), in the positive direction of the z-axis. 
205. Compute the density of the circulation of the vector 
a = 2yi + 5z2j around the ellipse L: {r = acost, y = 
= bsint, z= 1(0 < t < 2n)} at the centre C (0, 0, 1) 
of the ellipse in the positive direction of the z-axis. 
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Sec. 17. The independence 
of a line integral of the path 
of integration. Green's formula 


Definition. A region G of three-dimensional space is 
said to be simply connected (more precisely, it isa simply 
connected plane region) if any closed contour lying in the 
region can be spanned by a surface lying entirely in G. 
For example, the whole of .three-dimensional space and 
the interior of a sphere are simply connected regions; the 
interior of a torus and three-dimensional space with 
a straight line deleted are not simply connected regions. 

Theorem. In order that the line integral 


\ (a, dr) = j P (z,y, 2) dz+Q (z,y, 2) dy+R (z, y, 2) dz 
i L 


should not depend on the form of the path of integration L, 
it is necessary and sufficient that the vector field 


a= P(g, y, zji+Q(z, y, z)j+ Ro, Y, z)k 
be irrotational, that is, 
curl a (M) = 0. 
It is assumed here that the coordinates P (z, y, 2), 
Q (xz, y, 2), R(z, y, 2) of the vector a have continuous 


partial derivatives of the first order and the domain of 
definition of the vector a(M) is simply connected. -— 


In that case, the line integral J (a, dr) will depend 
solely on the position of the initial and terminal points 
of the path of integration L. 

If the theorem is complied with, “the circulation of 


the vector a(M) around any closed contour C located in 
the field of the vector a(M) is equal to zero: 


$a, dr) =0. 


c 
Example 1. Show that in the field of the vector 


a = zy?2i + xz*y2j + Satytk 
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the line integral j (a, dr) is independent of the shape 
L 
of the path of integration L. 

Solution. The coordinates of the vector a are everywhere 
continuous functions so that the domain of definition G of 
the vector a is the entire space (a simply connected region). 
In this region we have 


ij k 
0 0 7] 
curla=| @z dg @ |=0 


zy?z xyz ry? 
Consequently, the line integral 


See dr) = aye dz+ xtyz dy+ > xy? dz 

L 

is independent of the shape of the path of integration L. 
In particular, for the plane vector field 


a(M) = P(z, yi+Q, y)j (1) 
we have 
ij k 
curl a(M)= EER =(4-F)k. 
P@o0 


Therefore, for the plane vector field (1) that is defined 

in asimply connected region G, the condition curl a ( 

= 0 is written in coordinate form thus: 0P/dy = dQ oa. 
This means that in order for a line integral 


JPe, y)dz+Q(z, y) dy 


in a plane field defined in a simply connected region G to 
be independent of the shape of the path of integration, it 
is necessary and sufficient that the following relation hold: 


aP __ 2 
y Oz 

Remark. The requirement that the region G, where the 
vector a = a (M) is defined, be simply connected is essen- 
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tial. If the region G is nonsimply connected, then, provid- 
ed_that curl a (M@) = 0, the line integral may depend 
on the shape of the path of integration. 


Example 2. Let us consider the line integral 
—ydz zdy 
eet ee 


Solution. The integrand is meaningless at the point 
O (0, 0); and so we eliminate that point. In the remaining 
part of the plane (which is then a nonsimply connected 
region), the coefficients of dz and dy are continuous and 
have continuous partial derivatives, and the following 
identity holds true: 


6 z _ @ y 
Or (aaa) "ay ( —xta) ’ 
On the other hand, if we compute the integral around the 


circle L: 2? + y* = R?, then by parametrizing the equa- 
tion of the circle we get - 


on 2n 
—yd a R* sin? t+ R? cos? t 
§ eet ee — [tat | ar—2n. 


We find that the circulation is nonzero and, hence, the 

line integral depends on the path of integration. 
Determine in which of the vector fields indicated below 

the integral is independent of the shape of the path of 


integration: 
206. a = 24 + 27f + y’k. 
207. a= itu te 


Vitatepa 
_ yi zjtak 
208. a= “+ yi+a 


Green's formula 
Given in a region D with boundary L a plane vector 


field 
a= P(z, yit Qe, y)i, 


where the coordinates P (z, y), Q(z, y) are continuous 
and have continuous partial derivatives dP/dy, dQ/dz. 
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Then Green’s formula, 
Gp ast Qay= | I(a— ay) dray, (2) 
D 


holds true. Here, the boundary L is traversed so that the 
region D is on the left. 


—— 


a 


L, 
Fig. 33 


The region D may also be nonsimply connected so that 
the boundary may consist of several components (Fig. 33). 
In that case, the integral sy 


JPas+ Ody 


is understood to mean the sum of the integrals over all 
components of the boundary of D. 

Green's formula (2) is a special case of the Stokes theo- 
rem (see Sec. 16). 

In certain cases, Green’s formula permits simplifying 
the computation of the circulation of a vector field. 

Example 3. Compute the circulation of the vector 


a=Vitxe+y ity (zy +1n (c+ V1+2?4 y)I)j 


around the circle 27 + y? = R?. 
Solution. The circulation of the given vector is 


C= §(a, dr) = b Vi+z2-+ y2 dz 
L L 


+y [zy +1n (2+ V1+ 27+ y? )] dy. 
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Here 


P=Vitety, Q=aytyln(c+Vite+y*). 
We find the partial derivatives 
aE, 
wy Vitate or Yt 
Using Green's formula, we obtain 


— 
Vitaty 


C= \ (e+ yap yipeee) oe 


= j jy dz dy. 
D 
Passing to polar coordinates, 


=pcosp, y=psing, 
we have 


C= \ \p2sin? gp dpdp= |\ \ p* sin? pdpdg. 

i r 

Since O<S Q< a 0<oe = R, it follows that 
C= Jintao ot p= 20" 


Compute the circulation of the following vectors around 
the given contours using Green's formula: 
209. a=(yt+z)i+Yy—2zj 2: e+y=1, 2 =0, 
y =v. 
210. a= (x — y)i+ Qzyj; L: y=2, y= 2. 
zy 
211. a=aln (1 +y¥)i+7 pp ii L: 22 + y*? = 22. 
212, a= yi — 2}; L: ee yee —1,z=0,y=0. 
(32—y* VIF 284 4y*) 1+ (18y8+ 2° Viet tay) | 
213. a= = 
3V 1+ 29+4y3 
L: 22+y?=1. 
214. Use Green’s formula to compute the difference be- 
tween the integrals 
I= J @+ydz—(2—y)ay 


AmB 
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and 
I, | (e+ydz—(2—y)Pdy, 
AnB 


where AmB is a line segment joining the points A (0, 0) 
and B (1, 1), and AnB is an arc of the parabola y = 2’. 
215. Prove that the integral 


(22+ y)dz+ 2x dy, 


L 


where L is a closed contour, yields the area of the region 
bounded by that contour. 
216. Using Green’s formula, compute the line integral 


(a, dr) in the vector field 


a = (et sin y — y) i + (e* cos y — 1) j, 


where the curve L is the upper semicircle 2? + y? = 2r 
traversed in the direction from the point A (2, 0) to the 
point O (0, 0). 


CHAPTER IV 
POTENTIAL FIELDS 


Sec. 18. The criterion for the potentiality 
of a vector field 


Definition. A vector field 
a (M) = P (z, y, z)i + Q(z, y, i+ RC, y, zk 
specified in a spatial region V is said to be potential if 
there exists a scalar function @(M) such that at all 
points of V the following equality holds: 
a (M) = grad @ (M). (1) 

The function g (M) = @ (z, y, z) that satisfies (4) in 
V is termed the potential (or the potential function) of 
the vector field a. 

The relation (1) is equivalent to the following three 
scalar equations: 

OP: es J oY: 
P(z.y.2= Se, Q(zy2=se, Rizy, = 2. 
(2) 

The potential of a field is not defined uniquely but 
only up to an additive constant. 

Remark. For force fields, the function g (M) is ordi- 
narily called a force function, and the potential is the func- 
tion —q (@M). 

Example 1. (The electrostatic field of a point charge.) 
Show that the field of electric intensity E set up by a point 
charge g located at the coordinate origin, 


7 E=4r, r=V 22+ y2 + 22, 


is a potential field. 
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Solution. The problem is posed thus: show that there 
exists a function 9 (z, y, z) such that relations (2) hold. 
In our case, we have 


P(z.y.2=S, Oy ad=4, Rey, 2=%. 


Since 


Cn me | 4 Or z 
(=)= noe Pr 
and, analogously, 
$(t)=-4.4(4)=-5. 


it follows that the function 
a 
@ (z, y, 2) — r V + y3+23 
is the potential of the given field: 
grad ( —+) =E. 


In this case, the coordinate origin (this is where the 
charge q is located) is a singular point of the field E. 
Theorem. For a vector field a (M) specified in a simply 
connected region V to be potential, it is necessary and suf- 
ficient that the following condition hold at every point of V: 


curl a = 0. (3) 


In other words, for a vector field specified in a simply 
connected region to be potential, it is necessary and suf- 
ficient that it be irrotational. 

The potential  (z, y, z) of the vector field 

a= P(z,y, zi +Q(z, y, z2j+RG, y, 2k 
is defined by the formula 
(x, y, z) 
(2,y2)= | Pde+Qdy+Raz, (4) 
(Xo, Vo. 20) 
where (Zp, Yo, Zo) is some fixed point of the field and 
(rz, y, 2) is an arbitrary current point. 

Example 2. Show that the field of the vector 

a=rit+ yj + 2k 


is a potential field. 
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Solution. The coordinates P = 27, Q = y*, R = 2° of 
the vector a are infinitely differentiable functions through- 
out the space so that a is an infinitely differentiable 
vector defined throughout three-dimensional space. We 
have 


ij sk 
a) a a 
curl a= Ve Oy os =(42-#¥)I 
a2 y2 22 
6 a a 
+(¢2-Z2)1+(Zy—-Ze)k=0 


By virtue of the theorem given on p. 122 the field of the 
vector a is a potential field. It is readily seen that the 
function 


3 
9(z,y, 2) ==tEt* 4 ¢, 


where C is an arbitrary constant, is the potential of the 
given field. 

Check to see whether the following vector fields are 
potential fields: 
217. a = zzi + 2yj + zyk. 
218. a = (2ry + 2*) i + (Qyz + 2°) j + (222 + y")k. 
219. a = + (ci + yj + 22k). 


220. a = yz cos zy-i + 22 cos zy-j + sin zy-k. 
221. a = In (4 + 2") i + In (1 + 2?) j + ack. 
a 


a a= (Set) t4(S44) i+ (H+4)e 


Zz 
223. H= "0 (—yi+2j), Paot+y, rX0. 
224. Prove that the field a =f (r)-r, where f(r) is a 
differentiable function, is a potential field. 
225. Show that the vector lines in the potential field 


a = grad @ are perpendicular to the level surface of 
the function g. 
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Sec. 19. Computing a line integral 
in a potential field 


Theorem. A line integral in a potential field a (M) is 
equal to the difference between the values of the potential 
@ (M) of the field at the terminal and initial points of 
the path of integration: 

M, 

f (a, dr) =9 (M2) —9(M)). (t) 

M, 

Example 1. Compute the line integral in the field 
of the vector 
r= zi + yj + zk 

along the line segment bounded by the points M, (—1, 0, 
3) and M, (2, —1, 0). 

Solution. We will show that the field of the given vector 
is a potential field. Indeed, 


i j k 


It is easy to see that the potential of this field is 


Q(z, y, y= Stet? 1. 

Using formula (1), we obtain 

Ms 

{ (a, dr) =9(2, —1, 0)—9(—1,0,3)=2-5=—8. 
M, 


Note that it is immaterial what line joins the points M, 
and M,; for fixed M, and M, the integral 


Ms M, 
| (a, dr) = j zdz+ydy+2dz 
M, M, 


always has the same value. 
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Computing the potential of a field in Cartesian coordinates 


The formula 
Q(z, Y, 2) 
(x,y, z) 


= j P(z, y, 2)dz+ Q(z, y, 2) dy+R(z, y, 2) dz 
(Xe, Yor Zo) 
(2) 


may be used to find the potential function g (M) = 
= (z, y, 2) of a specified potential field 
a(z, y, 2) = P (x, y, 21+ Q(, y, 2)j 
+ R(z, y, 2)k. 


To do this, fix theinitial point M, (zo, Yo, 29) and join 
it with the current point M (z, y, z) with a broken line 


M(2,Y,2) 


Fig. 34 


M,ABM whose segments are parallel to the coordinate 
axes; namely, M,A || Oz, AB || Oy, BM || Oz (Fig. 34). 
Then formula (2) takes the form 

= y 
@ (2, ys 2) = | P(x, Yo, 20) d2+ [ Q(z, ys 20) dy 


Xe Vo 
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+ ja (z,y,2)d2z, (3) 
Ze 


where z, y, z are the coordinates of the current point on 
the segments of the broken line along which the integra- 
tion is performed. 

Example 2. Prove that the vector field 


a=(ytzjit+t@+2j+et+yk 
is a potential field and find its potential. 
Solution. Ist method. A necessary and sufficient condi- 


tion for the potentiality of a field a(M) is that curl 
a (M) be zero. In our case, 


i j k 
curl a= Oy Bs, 205 
= Oz oy Oz 
yt2z2+4+2a2+y 


=(1-—1)i+ (1-1) j+ (1-1) k =0. 
That is, the field is a potential field. The potential of 
the field can be found with the aid of formula (3). For 
the initial fixed point we take the coordinate origin 
O (0, 0, 0). We thus have 


x , y 2 
@ (2, ys 2)= (040) drt | (240) dy+ [(e+y) dz 
0 0 0 


= ry +22 + y2. 
To summarize, 
Q(z, y, 2) =zytaz+ y2t+C, 
where C is an arbitrary constant. 
2nd method. By definition, the potential @ (z, y, 2) 
is a scalar function for which grad g = a. This vector 
equality is equivalent to three scalar equalities: 


Bayts, (4) 
Bartz, (5) 
arty. 6) 
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Integrating (4) with respect to z, we obtain 
9 (ty. 2)=\(y+a)de=zy+22+fly,2), (7) 
0 


where f (y, z) is an arbitrary differentiable function of 
y and z. Differentiating both sides of (7) with respect to 
y and taking into account (5), we obtain a relation for 
finding the as yet undetermined function f (y, z). We 
have 


ap 
ent e 
or 
s+e=atZ, 
whence 
of 
a ay” (8) 


Integrating (8) with respect to y, we have 


y 
fy. 2) = |2dy=zy+F (2), (9) 
iy 
where F (z) is an as yet undetermined function of z. 
Substituting (9) into (7), we get 


@ (z, Y, 2) = zy + 22 + zy + F (2). 


Differentiating this equation with respect to z and taking 
into account (6), we obtain an equation for finding F (z): 


rty=2tyt, 


whence dF/dz = 0 so that F (z) =C = constant. 
Thus we have 


gp, yz) =szytyztarte. 


3rd method. By the definition of the complete differen- 
tial of the function 9 (z, y, z) we have 


dp= Fe drt 2 dy+ Baz. 


128 Potential Pields (Ch. IV 


Substituting in place of the partial derivatives d/0z, 
dq/dy, Op/dz their expressions taken from (4), (5), (6), 
we obtain 


dp = (y + 2) dz + (x + z) dy + (x + y) dz 
or, after some simple algebra, 
dp = (y dx + x dy) + (dx + xdz) + (y dz + 2 dy) 


= d (zy) + d (zz) + d (yz) = d (zy + zz + yz). 
Thus 
dp = d (ry + yz + 22), 
whence it follows that 


p(z, y, 2) =szytyztatc. 


In the following problems, establish the potentiality of 
the given vector fields a (M) and find their potentials 


p (M): 
226. a = Qryzi + 222j + ryk. 
227. a = (yz + 1) i + zzj + cyk. 
228. a = (Qzy + z) i + (2? — Qy) j + zk. 
_ i+j+k 
229. a= Sys: 
i k 
230, a = #4 silk 


231. a 


r 
; 
232. a=. 


233. a=r-r. 
When the region 2 is a star with centre at the coordi- 
nate origin O (0, 0, 0)*, the potential @ (M)of some 


* The region Q is called a star-shaped fepion relative to the 
point O belonging to Q if any ray emanating from that point cuts 
the boundary of the region at one poiiit at most. For example, star- 
shaped regions in the plane are the plane itself, a parallelogram, 
a circle; in three-dimensional space, the space itself, a parallele- 
piped, a sphere. 
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vector field a = a(M) at the point M (z, y, 2) may 
be found from the formula 


1 
9(M) = fea (M’), r(M))dt+-C, C=constant, (40) 
u 


where r (M) = zi + yj + zk is the radius vector of the 
point M (z, y, z), and the point M’ (tz, ty, tz)forO < t 
< 1runs over the segment OM of the straight line passing 
through the points O and M. 

Example 3. Find the potential of the vector field 


a= y2i + zzj + zyk. 


Solution. It is readily seen that curl a = 0, which 
means the given vector field is a potential field. This 
field is defined throughout three-dimensional space and 
is star-shaped with centre at the coordinate origin 
O (0, 0, 0), and therefore to find the potential we take 
advantage of formula (10). Since in this case 


a (M') = a (tz, ty, tz) = t?yzi + t?x2j + t2zyk, 
it follows that the scalar product of the vectors a(M’) 
and r (M) is equal to 
(a (M'), vr (M)) = t? (cyz + xyz + xyz) = 3i2zyz. 
The desired potential is 
1 1 
@(M) = ft (M’), r(M)) dt = zyz | 32° dt + C=ayz+C. 
0 0 


@ (M) = zyz + C. 


Using formula (10), find the potentials of the following 
vector fields: 

234. a = ai + Bj + yk, where a, B, y are constants. 
235. a= (y+ 2z)i+ (ec +2z)j+yta2)k. 

236. a = yi + 2j + e*k. 

237. a= e sin y-i + e*cosy-jf +k. 


Thus, 
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CHAPTER V 


THE HAMILTONIAN OPERATOR. 


SECOND-ORDER DIFFERENTIAL 
OPERATIONS. 


THE LAPLACE OPERATOR 


Sec. 20. The Hamiltonian operator del 


Many operations of vector analysis may be written in 
abbreviated form and in a form convenient for calcula- 
tions; this is done through the use of a symbolic operator 
called the Hamiltonian operator del: 


+i eke (1) 


This operator combines both differential and vectorial 
properties. We will regard the formal multiplication of 
d/dx by the function f (z, y, z) as the partial differen- 
tiation Of /dz. 

Within the framework of vector algebra, we will per- 
form the formal operations involving del as if it were 
a vector. Using this formalism, we obtain the following. 

4. Ifu =u (z, y, z) is a scalar differentiable function, 
then, by the rule of multiplying a vector by a scalar, we 
have 


. O . @ 0 
w=istigtks)u 


Ou 


. Ou , . du 
=ig - tig; +k z-=eradu. (2) 


2. fa=P(z,y, 2z)i+ QO, y, z)j + R(z, y, 2)k, 
where P, Q, A are differentiable functions, then, by 
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the familiar formula for a scalar product, we have 
a,:;4@ a F ‘ 
(Via)= (I +iZ +ks, Pit Qi+ Rk) 
_ OP , 00, aR 4. 
=a tay tae =diva (3) 
in particular, (V, c) = 0, where c is a constant vector. 


If a = P(z, y; z)i+ Q(z, y, z)j+ R(z, y, z)k, 
then 


{V, a]=|=— = = |=curl a; (4) 


in particular (V, c) = 0, where ¢c is a constant vector. 

Continuing the formalism of operations with V as 
a vector, we obtain the following from the distributive 
property of scalar and vector products: 


(V, a + b) = (V; a) + (V, b), (5) 
that is, div (a + b) = div a + div b 
(v, a+ bl =Iv, al+[v, bl] (6) 


or curl (a + b) = curl a + curl b. 

Formulas (5) and (6) may also be interpreted as an 
exhibition of differential properties of the del operator 
(vV is a linear differential operator). 

When using the formalism of operations involving the 
del operator regarded as a vector, one must bear in mind 
that del is not a vector, for it has neither magnitude 
nor direction, so that, for example, the vector [V, al 
will not, generally speaking, be perpendicular to the 
vector a (however, for the plane field a= P (z, y)i+ 
+ Q(z, y) j the vector 


[V, a}=curl a= (2-+) k 
will be perpendicular to the zy-plane and, hence, to the 
vector a). In the same way, the concept of collinearity is 


meaningless with respect to the symbolic vector V. For 
example, the expression [Vp, Vip], where @ and p are 


Qe 
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scalar functions, formally resembles a vector product of 
two collinear vectors, which product is always equal to 
zero. But this is not true in the general case. Indeed, the 
vector Vo = grad @ is directed along the normal to the 
level surface @ = constant, while the vector Vw = grad p 
defines the normal to the level surface p = constant, 


Fig. 35 


and in the general case these normals are not necessarily 
collinear (Fig. 35). On the other hand, in any differentiable 
scalar field p we have [Vg, Vp] = 0. These examples 
show that the del operator must be handled with care. 

Besides its vectorial nature, the Hamiltonian operator 
de! has a differential aspect to it. Taking into account 
the differential aspect of V, we will agree that the operator 
V acts on all quantities that follow it. In this sense, 
(V, a) (a, V). Indeed, 

(V, a) = diva, 

whereas 


C] 7] 7] 


is a scalar differential operator. 
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When applying the del operator to a product of any 
quantities, one must bear in mind the rule for differenti- 
ating a product: 


0 Ou Ou. 
Gq WSU tus 


From this it follows that the del operator must be applied 
in succession to each factor, leaving the other factors 
unchanged, and then the sum of the resulting expressions 
taken. In this procedure we are guided by the following 
rules. 

1°. If the del operator acts on some product, first 
take into account its differential nature and only then 
its vectorial properties. 

In order to note the fact that the del operator 
does not act on some quantity involved in a complicated 
formula, that quantity is labelled with the subscript c 
(constant), which may be removed in the final result. 

3°. All quantities not acted upon by the del operator 
are placed in front of the operator (that is, to the left) 
in the final result. 

Example 1. Show that 


div (ua) = u div a + (a, grad uw). 
Here, wu is a scalar function and a is a vector function. 
Solution. In symbolic notation, 
div (ua) = (V, wa). 
Taking into account first the differential nature of V, 
we have to write 
(Vv, wa) = (V, u,a) + (Vv, wa,). 
Regarding the expression (V, u,a), we can take the con- 


stant factor u, outside the del sign and, as a scalar, outside 
the sign of the scalar product; this yields 


(V, u,a) = (u.V, a) = ue, (Vv, a) =u (Vv, a). 


(The last step is to drop the subscript c). 
In the expression (V, wa,), the del operator acts only 
on the scalar function uv, and so we can write 


(V, ua,) = (Vu, a.) = (a,, Vu) = (a, Vu) 
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to get the formula 
(Vv, wa) = u(V, a) + (a, Vu) 
or 
div (ua) = u div a + (a, grad u). 
Example 2. Show that 
curl (ua) = u curl a — [a, grad uJ. 
Solution. In symbolic notation, 
curl (ua) = (V, wa). 
Taking into account the differential properties of V, we 
first write 
[v, ua] = [V, u.a] + [V, va,]. (7) 
Then in the first term on the right we take the scalar 
factor u, outside the del sign and outside the sign of the 
vector product, which yields 
Iv, ua] = u, [V, a] = ulv, al. 
In the second term in (7) we refer u to the operator V and 
change the order of the factors so that the vector a,, which 
del does not act on, is in front of V. This yields 
[v, wa,] = [vu, a] = —[a, Vu). 
Thus 
[v, ua] = ulv, a] — la, Vu) 
or 
curl (ua] = u curl a — [a, grad ul. 
Example 3. Use the symbolic method to find div [a, b] 
Solution. We have 
div [a, b] = (V, [a, bl) = (V, [a, b,]) 
+ (Vv, [a., b)). (8) 
Using the property of cyclic permutation of factors in 
a mixed product, we transform the expression on the 
right of (8) so that all constant quantities are in front 
of the del operator and all variable quantities follow it. 
This yields 
div(a, b]=([V, a], b.)—(V, [b, ac]) 
=({[V, a], b.)—([V. b}. a) =(b. [V, a]) —(a, [V, b]) 
or div{a, b] =(b, curla)—(a, curl b). 
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Remark. The use of the symbolic method enables us to 
avoid cumbersome analytical transformations and obtain 
the final results very quickly. On the other hand, the 
various formal transformations involving the del opera- 
tor must be performed with extreme caution in order to 
avoid serious mistakes. For this reason, if there is any 
doubt about the final result, it is wise to verify it by 
the analytical method. 

238. Show that 

(a) v(+)=4 

(b) Vf (u) = f’ (u) Vu. 

239. Prove that the vector [Vu, Vv] is solenoidal if 
u and v are differentiable scalar functions. 

Use the Hamiltonian operator del (V) to prove the 
following equations: 

240. (a) grad (uv) = v grad u + wu grad v; 


(b) curl (a, b) = (b, V) a — (a, V) b + a div b — 


— b diva. 

241. curl [a, r] = 2a, where a is a constant vector. 
242. Prove that the vector a = u grad v is orthogonal to 
curl a. 


Sec. 21. Second-order differential 
operations. The Laplace operator 


Second-order differential operations are obtained as 
a result of a twofold application of the operator V to 
fields. 

Suppose we have a scalar field u = u(M). In this 
field, the operator V generates a vector field Vu = grad u. 

In the vector field Vu, the operator V, when applied 
a second time to Vu, yields the scalar field 


(V, Vu) = div grad u (1) 
and the vector field 
(V, Vu) = curl grad u. (2) 


If a vector field a = a(M) is given, the operator V 
generates in it a scalar field (V, a) = div a. In the scalar 
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ield div a the operator V generates a vector field 
Vv (Vv, a) = grad div a. (3) 


In the vector field a = a (M), the operator V also gener- 
ates a vector field [v, al = curl a. If the operator V is 
again applied to this field, we obtain the scalar field 


(v, [Vv, al) = div curl a (4) 
and the vector field 
[v, [v, al] = curl curl a. (5) 


The formulas (1) to (5) define what are called differen- 
tial operations of the second order. 

Example 1. Suppose a function u = u (z, y, z) has 
continuous partial derivatives up to second order inclu- 
sive. Prove that 


curl] grad u = 0. 
Solution. Ist method. Operating formally, we obtain 


curl grad u = {[v, Yul = Iv, vVlu=0 


since [V, V] = 0 being the vector product of two identical 
“vectors”. 

2nd method. Using the expressions for the gradient and 
curl in Cartesian coordinates and taking into account the 
given conditions, we have 


ijk 
7) 06.620 
eam es: 0 { du 0 { du 
curlgradu= | a %|=|3,(32)—ze (47) | 
au du au 
Or Oy a2 


+[()-2 (2) @)-$ (Es 
= ( dtu au yi atu Au )i ( a3u d%u )k=0, 


ay dz az dy a20n Oz az azdy dyor 
since the mixed derivatives are equal in this case, 
au atu au au 8u tu 


dydz d2dy’ d20x  dzdz* Odzdy dyad’ 
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In similar fashion proof is given that for the vector 
field 


a=P(z,y,z)i+Q (zy 2j+ Re, y, 2k, 


the coordinates P, Q, R of which have continuous partial 
derivatives of the second order, we obtain div curl a = 0. 

Note particularly the second-order differential opera- 
tion div grad u = (V, Vu). Assuming that the function 
u (z, y, z) has second partial derivatives with respect to 
Zz, y, 2, we obtain 


on : a .. 8 


= E(B) (B) +E (B) Bi emo 


Thus, (V, Vu) = Au, where the symbol 
_ fu , Au Hu 

Au=sat+ att aa 
is termed the Laplace operator (or Laplacian). It may be 
represented as a scalar product of the Hamiltonian oper- 
ator V into itself, that is, 

a a, a 
A=(V, V)=V?= 528+ at oa: 

This operator plays an important role in mathematical 
physics. 

Let us examine another second-order operation curl 
curl a. We have curl curl a=[v, [v, al]. Let us take 
advantage of the formula for a vector triple product writ- 
ten as [A, [B, C]] = B(A, C) — (A, B) C. Replacing 
A by V, B by V, and C by a, we obtain 


[v,[v, all =vV(v, a)—(V, V) a=V(V, a) — Aa (6) 


or 
curl (curl a) = grad div a — Aa, where Aa = 


= AP-t + AQ-j + AR-k. 


The following table is a pictorial display of second- 
order differential operations: 
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Sealar field u Vector field a 


grad div curl 


+ grad div a 

7 

z div grad u=Au |__| ___ div curl a=0 

a duet curl curl a 

3 eur prea =grad div a— Aa 


Example 2. The laws of the classical theorem of electro- 
magnetism are postulated in the form of a system of 
Maxwell equation. 

In the most elementary case of a nonconducting homo- 
geneous and isotropic medium and in the absence of 
charges and currents, this system is of the form 


+ & =Iv, Hi). (7) 
—~+ FB iiy, E), (8) 
(V, E)=0, (9) 
(Vv, H)=0. (10) 


Here, E and H are vectors of the electric-field and mag- 
netic-field intensity; e and p are coefficients of the permit- 
tivity and permeability (our assumptions are that e 
and p are constants); c is the velocity of light in empty 
space. 

Since the spatial and time derivatives commute, that is, 


lv, H=[v, 3], 
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it follows that by differentiating (7) with respect to ¢, 
we obtain 
e FE = oH | 


c of | %* at 
Replacing ae taken from (8), we get < a = 
= - £19, Iv, Ell or 
oo SF —(v, IV. Ell. (11) 
By virtue of formula (6), [V, |v, E]] =V(V, E) — AE. 
Since (V, E) = 0, it follows from (41) that 4 o= = AE. 


To summarize, for the vector field E we obtain the 
equation 
OE 3 
“Of ~ en 
This is one of the basic equations of mathematical physics 
and is called the wave equation. 
It is easy to see (check this!) that the vector field H 
satisfies the same kind of wave equation 


Cds | e 
ot ~~ en AH. 


Thus, under our conditions, each of the coordinates 
E,, E,, E, and H,, H,, H, of the vectors E and H 
satisfies the equation 


Au Bu tu 
— = q2 f —— pant a: 
ot? e (aa+art ar) 


Here, a = c/Y ep is the velocity of propagation of the 
process. In a vacuum, where e = p = 1, we have a = ¢, 
that is to say, in a vacuum, electromagnetic processes are 
propagated with the velocity of light. 

243. Show that any solution of the equation [v, [v, All — 
— HA = 0 that satisfies the solenoidal condition satisfies 
the Helmholtz vector equation 


VA + A = 0. 
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Definition. The scalar field u = u (z, y, 2) that satis- 
fies the condition Au =0 is termed a Laplace (or harmon- 
ic) field. 

Example 3. An important instance of a harmonic field 
is the scalar field u=k/r, k=constant, r= V 2? + y? + 2?. 
This function is the potential of a gravitational 
field generated by a point mass placed at the coordinate 
origin. It is‘easy to verify that the function u = k/r is 
harmonic everywhere except at the coordinate origin, 
where it is not defined. Indeed, 


= 3, FAV, 1) = 
for all r 40 since 
1 1 
(V, M=(V,)=(V—, r)+—=(V, 9) 


S(-ae i ets 


Example 4. Prove that in the potential field of a vec- 
tor a its potential function u (z,°y, z) satisties the Pois- 
son equation 


au dtu, at 
Aus Fat Fat sane (t Ys 2), (12) 
where p (z, y, z) is the divergence of the vector a. 
Solution. It is given that 
diva=o. (43) 


Since the field of the vector a is a potential field, 
it follows that a = grad u, where u is the potential of 
the field. Substituting a = grad u into 13), we obtain 
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div gradu =o or, since div grad u = Au, we have 
Au = p. 

In the special case of points of the field where the diver- 
gence is zero, equation (12) turns into the Laplace equa- 
tion Au = 0. The Laplace-Poisson equation permits find- 
ing the potential function u by integrating the partial 


BW 


++t+4+4+] [+++ 4+4+ 444 


TTT TITS TLTSILISIIVA YT OOO OO La aaa ahhh hh chchinhenkn 


A 0 Ay 
Fig. 36 


Wy 


differential equation. This turns out to be more convenient 
in some problems. 

In electrostatics, preference is often given to finding 
the function v = —u instead of u. Then a = —grad v. 
Accordingly, in the theory of the electrostatic field, the 
Poisson equation is of the form 


ay aw , Fv 4np 
att oat ae me a) 

Let us consider an elementary case where the Poisson 
equation is employed. 

Example 5. Suppose two infinite parallel plates AA, 
and BB, with opposite charges have potentials v, and v,; 
for the sake of definiteness, v, > v,. Find the field E 
between them (Fig. 36). 

Solution. Send the z-axis at right angles to the plates 
in the direction of decreasing potential, and bring the 
yz-plane to coincidence with the positively charged plate 
AA,. We now seek the potential function from the Pois- 
son equation. By virtue of the symmetry of the problem 
about the z-axis and due to the infinity of the plates, 
we can conclude that the equipotential surfaces are planes 
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parallel to the plates, and the function v depends solely 
on the variable z. Equation (14) takes the form 

a? 

= 0 (45) 
since space charges are absent throughout the space be- 
tween the plates. Integrating (15), we find 

v=C,r4+C,, (16) 

where C, and C, are arbitrary constants. 

We require that for z = 0 the function v take the value 
v, and for z = d, where d is the distance between’ the 
plates, that it take the value v,. This yields C,=v,, vs= 
= Cyd + C2, whence C, — Vi, C, = (v, — v,)/d. Substi- 
tuting these values of C, and C, into (16), we obtain 

Ug—v, 


v= z. 


vy U. 
I+U,;=VYy— 7 2 


The vector E is found from the formula E = —grad v, 
which yields 


E=4>*i 


so that the field is homogeneous and in the direction of 
the z-axis. The magnitude of E is equal at every point to 
| E | = (v, — v,)/d, that is, it is equal to the potential 
drop per unit length of the shortest distance between the 
plates. 

244. Suppose a scalar function @ (M) satisfies the Laplace 
equation. Show that the vector V@q is solenoidal and irrota- 
tional. 

245. Show that A (wv) = u Av + v Au + 2 (Vu, Vv). 
246. Prove that if r is the radius vector, then 


2 in space, 
Ar= r 
= in the plane. 


247. Check to see whether the following scalar fields 
are barmonic or not: 

(a) us z+ Ary — y*, 

(b) u = zy + y*2 + 22, 

(c) u=2— y’*. 
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248. Show that the scalar field 
u=in+, where r=/Vz?+y2 (r0), 


is harmonic. 
249. Find all harmonic fields that depend solely on z. 
250. Find the general form of a homogeneous harmonic 
polynomial of second degree in z and y. 
251. Find all solutions of the Poisson equation Au = 
= 2"? that depend solely on z. 

Example 6. Green's formulas. Let 9, » be two scalar 
functions of a point. Set up the vector a = 9 grad yp. 
Then 


div a — div (9 grad p) = @ div grad p + (grad q, grad »p) 
= o- Ap + (grad 9, grad»). 
Now apply the Gauss-Ostrogradsky formula 


9 (a, n°) do = J) j div adv. 


Note that in our case 


(a, n°) =(pgrad p, n°) = (grad p, n°) =p, 
We thus obtain Green’s first formula: 


JJ Jieav+(erade, grad wiav= Q) 95s oft do, (17) 
which for @ = tp turns into 
JS) (9 be +lerad odo = Q) 9 5 a) 
If we put p = 1 in (17) we get 
J) j Aydv =) 2 as. 


In (17), interchange p and » and subtract the resulting 
formula 


V5 [p Ag+ (grad p, grad @)) dv = Ove 9% ao 
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from (17). This then yields Green’s second formula: 


Jf J ay van d= Q | pt _ yp 2) as. 
f z 


It is assumed here that all functions that we have to deal 

with and also all their derivatives that occur in the for- 

mulas are continuous in the region under consideration. 
Example 7. Find the surface integral 


I= (0-5 a0 


taken over the sphere = ’?+y+2=1 for 9 = 
= 27+ y* and p = y? + 27. 

Solution. By Green’: s first formula, the desired integral 
is 


=|) | (pay+ (grade, grad y)} dv, 
v 
where the region of integration V is a sphere: 2? + y? + 
+2 < 1. 
We have Ap = 4, grad p = 2zi + 2yj, grad p = 
= 2yj + 22k, (grad g, grad p) = 4y? and therefore 


I= j \j (422 + 4y2 + 4y2) dv =4 hf 5 (22 + 2y2) dv. 


Passing to the spherical coordinates z = rcos 9 sin 0, 
y =rsin pmsin 6, z = rcos 8, we obtain 


IT=4 ( j j (r2 cos? p sin? 6 + 2r? sin? p sin? 8) r2 sin @ dr dO dp 
v 
on x i 
=4 j (cos? p + 2 sin? p) dp j sin? 0 d0 | rédr 
0 0 0 
on 7 
= j (1-4 sin? @) dp j (1 —cos?6) sin 6 d0 
0 0 


ea oe 1 39) |*_ 16 
=> x ( cos 8+ — cos 6) on Be 
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Example 8. Find the surface integral 


I= Y ( (o—» 2) do 


taken over the Stee Z:2+y= R?,2=0,2=4H, 
H > 0, provided g = 2? + y? Peay e+ y + 
+ 22+ 2. 

Solution. By Green's second formula, the desired inte- 


gral is 
I= { \J (p Ap— Ag) dv. 


For the given functions @ and we have Ag = 4, Ay = 4 


and therefore 
I=—4 J \J 2 dv. 


Passing to the cylindrical coordinates z = p cos g, y = 
=p sin Q, 2 = 2, we obtain 


= -4§{f napdode= 49 to | pap | as 


. = —2nR?H?2, 
Example 9. Find the surface integral 


I= Des 2 ao 


over a closed surface 2 pounded by the planes z + y + z= 
=1,z2=0, y=0, z2=0, provided p =e siny + 1. 

Solution. The given function is harmonic since Ag = 
= eé* sin y — e* sin y = 0. Therefore, by (18) we have 


ja ny) | grad @ |? dv. 
We find the modulus of the gradient of the function 9: 
grad p = ée sin y-i + e* cos y-j, | grad p | = & 
The desired integral is — to 
1 -x- 
I= e2* dy= | e2* dz ax fa [ dz= — 9). 
jieen| 


10-910 
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252, Compute the surface integral J = § og a do over 


a closed surface 2: ap piest y=0, y>0}, 
provided p=z?, p=2?+y?—2?2. 


253. Compute the surface integral J/= { (o x —» a ) do 


z 
taken over the entire surface of a closed cylinder &: 
{x24 y2=1, z=0, z=1}, provided p=22?, p=2?+22. 


254. Compute the surface integral J = em do, pro- 


vided p=(z+y-+2)/V 3and & isa sphere: pense R. 
255. Find the surface integral J = Ys do, provided 


p=e*siny+e"sinz+z and 2 is ; triaxial ellipsoid 
(2?/a2) + (y?/b2) + (22/c2) = 1. 


Sec. 22. Vector potential 

Let a vector field 
a=a(M)=P(z,y,2)i+Q(z,y,2)) + A(z, y,2)k 
be solenoidal in the region G, that is, div a (M) = OinG. 

Definition. The vector potential of a vector field a = 
= a (M) isa vector b (M)=P, (z, y, 2) i+ Qy (2, y, 2) + 
+ R, (z, y, 2) k that satisfies in G the condition 

curl b (M) = a (M) (1) 

or, in coordinate form, 
aR, 9Q, _ oP, OR, _ aQ oP 
Dg ee ee be ae ay 

For a solenoidal vector field a (J) the an potential 
b (M) is not defined uniquely: the condition (1) is also 
satisfied by the vector B(M) = b(M) + grad f (M), 
where f (//) is an arbitrary differentiable scalar function, 
since curl (grad f (M)) = 0. 

Thus, two vector potentials of the solenoidal field a (M) 
differ by the gradient of the scalar field. 
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Finding the vector potential b (7) of the solenoidal 
field a (M) reduces to finding some particular solution of 
system (2) of three partial differential equations for the 
three unknown functions P, (z, y, z), Q, (2, y, 2), 
R, (z, y, 2). 

The vector potential b (M) may be constructed in the 
following manner: Taking advantage of the arbitrariness 
of choice of the vector b (M),. we will simplify matters by 
setting P, (z, ys z) = 0, that is, the vector b (M) will 
be sought in the form b (M) = Q, (a, y, 2) j + 
+ R, (z, y, z)k. Then the system of differential equa- 
tions (2) for finding the unknown functions Q, (z, y, 2) 
and R, (z, y, z) takes the form 


OR =O: (3) 


From the second and third equations of this system we 
find 


Ry (2, ¥, 2)=—| Q(z, y, 2) dz +C, (y, 2), 
Orla, y, 2)= J R(x, y, 2)dz+Cr(y, 2) 


where C, (y, z) and C, (y, 2) are any differentiable func- 
tions of y and z. For the sake of simplicity, set C, (y, z) = 
= 0 and choose the function C, (y, z) so that the first 
equation of system (3) is satisfied. To do this, we sub- 
stitute into the first equation the expressions that were 
ia for Q, and Ry: 

ac, @ 


Hy J Ole, y, 2)dz+4—F | A(z, y, 2) dz 
=P (z, y, 2). 
From this we get 


muZ j Q(z, y, 2)dz+5— 5 | R (z, y, z)dx+P(z, y, 2). 


fo* 
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It is easy to verify that the right-hand side of this 
equation does not depend on z; this is because div a (M) = 
= 0 in G. 

Integrating the last equation with respect to y, we find 


Ci =S [J Oley 2) de 
+3 J ley, 2)de+ Plz, y, 2) ]dy4 Cals). (A) 


Setting C, (z) =0 in (4) and substituting (4) into the 
expression for R, (z, y, z), we get a particular solution 
of system (3): 


P,=0, (5) 
Qi= J Rl, y, 2) ae, (6) 
Rr=V[Z JQ y ade+Z | Riley, 2)de 


+P(z, y, 2)]ay—J Q(z. y, 2)dz. (7) 

The vector b(M), whose coordinates P, (z, y, 2), 
Q, (x, y, 2), Ry (z, y, 2) are defined by formulas (5), 
(6), (7), is the vector potential since it satisfies the condi- 
tion curl b = a. -- 

Example 1. Find the vector potential b = b (z, y, 2) 
for a solenoidal field given by the vector 

a = 2yi — aj + Qzk. 
Solution. We seek the potential b in the form 
b = bg, y, 2) a Qi (t, y, z) j + R,(z, y, 2) k, 

where Q, (z, y, z) and A, (z, y, z) are found from (7) 
and (8). Since in the given case P = 2y, Q = —z, R= 
= 2z, we obtain 


Q(z, ¥, 2= j 22 dz=2z?, 


Ry(z, y, a= j zdz+ j 2y dy =2z+ y?. 


Thus 
b (zt, y, 2) = 2j + (zz + y*)k 
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It is clear, by direct verification, that curl b = a, 
and, hence, this vector is the vector potential of the given 
field. 

Remark. Due to the arbitrariness in choosing the vector 
b, we could require Q, (z, y, z) =0 or R, (z, y, 2) =O 
instead of P, (z, y, 2) = 0. The system of equations (2) 
and formulas (5), (6), (7) would naturally change. 

Find the vector potentials of the following solenoidal 
fields: 

26. a=i+j+k. 

257. a = Qyi + 277. 

258. a = (e* — ev) k. 

259. a = 6y7i + 62) + 62k. 

260. a = 3y*%i — 3z7j — (y? + 2z)k. 

261. a = ye™i + Qyzj — (Qryze™* + 2%) k. 

If the vector field a = a (M) is solenoidal in the re- 
gion G, which is star-shaped (see Sec. 19, Chapter IV) with 
centre at the coordinate origin O (0, 0, 0) [the field a (M) 
may not be defined at the point O], then one of the vector 
potentials b = b(M) may be found from the formula 


1 
b(M) = { (a(M"), (Myler, (8) 
0 


where r (M@) = zi + yj + zk is the radius vector of the 
point M (z, y, 2) and the point M’ (tz, ty, tz) runs 
over the line segment OM as the parameter ¢ varies 
from 0 to 1. 

Example 2. Use formula (8) to find the vector poten- 
tial of the solenoidal field 


a = 2yi — a + Qzk. 


Solution. The given vector field is defined throughout 
three-dimensional space, which is a star-shaped region 
with centre at the coordinate origin, and so we can use (8) 
to find the vector potential. At the point M’ (tz, ty, tz) 
we have 


a (M’) = 2tyi — tz + 2tzk. 
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We find the vector product 


i j k 
[a(M'), r(M))=| 2ty —tz 2tz 
yy 


= [— (2ry + 22) i+ (22% — yz) j + (2y? + 22) k} et. 
Using (8), we obtain 
1 
b(M) = | | —(2ay + 22)i + (222 2y2)§ + (2y?+ 22) k| dt 
0 
= — 4 (Qay + 2142 (2? ~y2)j++ (2y?+ zz) k. 


It is easy to establish that curl b (M) = a (M). 

Remark. In examples 1 and 2 we obtain different vector 
potentials for one and the same solenoidal field a = 
= Qyi — zj + 2Qck: 


by (M) = 2%) + (22 +y?)k, 
be (M) = —-3 (ay + 2%) +5 (22—y2)j +4 (2y?+ 22) k. 


They differ by a term equal to the gradient of some 
scalar field f (M). This term plays the role of an arbitrary 
constant (when acted on by the curl). It may be repre- 
sented as the gradient of some scalar function_f (M). 
Let us find this function in our example. We have 


grad f (M) =b, (M) —b, (M) 
= 5 (20y + 22) i+ (e+ Qys) j++ (222+ y)k. 
To find the scalar field f (M), use formula (3) of Sec. 19, 


in which we take the coordinate origin O (0, 0, 0) for 
the point (Zo, Yo, 29). This yields - 


x v z 
f(M)= j O-de+ { Fedr+ j 4 (202+ y?) dz+C 
0 0 0 


= (xy + yz +22) +C, 


where C is an arbitrary constant. 
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Example 3. Find the vector potential b of a magnetic 
field H set up by an electric charge e that is moving with 
a constant velocity v. 

Solution. By the Biot-Savart law, the intensity of the 
magnetic field is 


H(M) = 31 (9) 


where r is the distance of point M from the charge e. 

Since H is a solenoidal vector, that is, div H = 0, 
there exists for it a vector potential b such that H = 
= curl b or, taking into account formula (9), 


(ev, r)_ e [v, r] 
curl b= 4ne 4x PrP ° 


Let o rewrite this formula as 


curl b= 7— + {[ev. St] +[ev, 4 = se) 
=e {li -]+[5. a a as a 
ar {lt &(F)] +a (F)] +e eG 


Employing the readily verifiable equation 


curla= [i, = ]+[i, 3 ]+[k 


we obtain 


x 48 
—— 


)}}. 


1 Vv 
curl b=—— curl —_, 
4n r 
whence 


1 ev 
b = Gn = 7 7? 

Using formula (8), find the vector potentials of the 
following solenoidal fields defined in star-shaped regions: 
262. a = i. 
263. a = 62i — 15yj + 9zk. 
264. a = Sa*yi — 10zyzk. 
265. a = 2 cos zz-j. 


—yitzi 
266. a= a 2+y?>0. 


CHAPTER VI 


CURVILINEAR COORDINATES. 
BASIC OPERATIONS 
OF VECTOR ANALYSIS 
IN CURVILINEAR COORDINATES 


Sec. 23. Curvilinear coordinates 


In many problems it is more convenient to define the 
position of a point M in space by three numbers (4q;, 73, 9s) 
instead of the three Cartesian coordinates (z, y, z). These 
numbers often prove to be more suitable to the problem at 
hand. 

Let every point M be associated with a definite number 
triple (q,, 92, 93) and, conversely, let every number triple 
correspond to a unique point M. Then the quantities 
91: Ya, 73 ave termed the curvilinear coordinates of the 
point 

The coordinate surfaces in a system of curvilinear 
coordinates 4;, Yo, 73 aré the surfaces 


mW =C,, (4) 
qe >= C., (2) 
I= Cs, (3) 


on which one of the coordinates remains constant. 

The line of intersection of two coordinate surfaces is 
called a coordinate line (azis). 

The coordinates g, and q, maintain constant values 
along the line of intersection of the coordinate surfaces 
(2) and (3); it is only the coordinate g, that varies. Simi- 
larly, on the lines of intersection of the surfaces (1) and (3) 
and (1) and (2), it is g, and gs that vary respectively. 

Let us introduce the unit vectors e,, e,, es directed 
along tangents to the coordinate axeg (q;), (G2), (gs) at 
the point © in the direction of increasing variables q,, q., 
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qs Tespectively (Fig. 37). Let us agree to take the unit 
vectors @, @,, es always in that order so that, taken 
together, they constitute a right-handed trihedral. 
The basic difference between curvilinear coordinates and 
Cartesian coordinates is this. In the Cartesian system, the 
vectors €,, @€,, €3 are constant for all points of space 
and are equal, respectively, to i, j, k. In any other 


Fig. 37 


system, they will, generally speaking, change their direc- 
tions when passing from one point, ™M, to other points. 
’As examples of curvilinear coordinates we consider 
cylindrical and spherical coordinates. 
1°. Cylindrical coordinates. The position of a point M 

in space is defined in cylindrical coordinates by three 
coordinates: 

1 = ?P» 0 < p < +0, 

@=9, DO Q9<2n, (4) 

qs = 2, —oo << z< +00. 


The coordinate surfaces are: 
p = constant: circular cylinders with the z-axis; 
@ = constant: half-planes adjoining the z-axis; 
z = constant: planes perpendicular to the z-axis. 
The coordinate lines (or axes) are: 
(p): rays perpendicular to the z-axis and having their 
origin on that axis; 
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(gp): circles with centre on the z-axis and lying in 
planes perpendicular to that axis; 
(z): straight lines parallel to the z-axis (Fig. 38). 
Cartesian coordinates are 
related to cylindrical coor- 
dinates by the following 
formulas: 


Z = Cos Q, 
y=psing, (5) 
2=2. 


2°. Spherical coordinates. 
In spherical coordinates, 
the position of a point M 
in space is defined by the 
following coordinates: 


aT, 0<r<+o, 
q. = 9, 0<6<xa, (6) 
Z@=97, DOK Q<2n. 


The coordinate surfaces 
are (Fig. 39); 
Fig. 38 r = constant: spheres 
centred at O: 
® = constant: circular half-angle cones with the 
2-axis; 
~ = constant: half-planes adjoining the z-axis. 
The coordinate lines are: 
(r): rays emanating from the point O; 
(9): meridians on a sphere; 
(p): parallels on a sphere. 
Cartesian coordinates are related to spherical coordi- 
nates via the following formulas: 


Zz =rcos @sin @, 
y =rsin g sin 9, (7) 
2=rcos 80. 
A system of curvilinear coordinates is said to be or- 
thogonal if at every point M the unit vectors C1, 2, 05 
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are pairwise orthogonal. In such a system, the coordinate 
lines and the coordinate surfaces are also orthogonal. 
Systems of cylindrical and spherical coordinates are instan- 
ces of orthogonal curvilinear systems of coordinates. 


Fig. 39 


Henceforth we consider only orthogonal systems of coor- 
dinates. 
Suppose r =r (q,, Y2, Ys) is the radius vector of 
a point M. Then 
dr = H, dq, &, + Hy dq, e, + Hy gz eg. (8) 


Here 
m=V (x) +(AZ)+(E). tt. 2,3 


are the Lamé coefficients of the given curvilinear system 
of coordinates. 
In cylindrical coordinates, 


a=? I=, Ia = 2. 
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By virtue of (5) we ae 


nome Ue ey Te yt =p, 


nan, (ECR CE) =! 


In stile coordinates, 
n=!, q2 = 8, Is = 9- 
By virtue of (7) we have 


eC) 


68 
_p _ az \2 dy\2 , (oz \2 
H,=H,= V (=) + (32) + (ss) =rsin@, 
The quantities 
dl, = H, dq, i=1, 2,3 

that appear in formula (8) are differentials of the arc 
lengths of the coordinates lines, In a number of cases, this 
consideration permits of a more simple computation of 
the Lamé coefficients. For instance, in the case of cylin- 
drical coordinates (4) (see Fig. 38), the differentials of 
the arc lengths of the coordinates lines (p), (@). (z) are 

d(p) =1-dp, whence H, = 1; 

d(~) = p-dg, whence H, = 9; 

d(z) =1-dz, whence H, = 1. 
It is just as easy to obtain expressions for the Lamé 
coefficients in the case of spherical coordinates (6). 


Sec. 24. Basic operations 
of vector analysis in curvilinear 
coordinates 


1°. Differential equations of vector lines. Suppose we 
have a vector field 


& = 4G; (91, 2, Fs) C1 + Fe (915 Yas Ys) Ce + 49 (91s V2» Ys) Cs- 
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The vector-line equations in curvilinear coordinates 
91» V2. 4g are of the form 
Hy dq = Hy dq, = Hs dq 
91 (911 92» 9s) 22 (91> Gar F9) 3 (91 Gar 95) 
In particular, in cylindrical coordinates (¢, = p, 9. = 9, 
Ya = 2): 
ee, Onn 5, eee ee (4) 
a, (P, @, 2) a, (P, 9, 2) a3 (P, 9, z) ’ 
in spherical coordinates (¢, = r, gg = 9, gy = 9): 
dr = r d0 _ _rsin6d@ 
a,(r, 8, @) ~— ag(r, 8, @) ~— ag (r, 8, @) * 
_ Example 1. A vector field is given in cylindrical coor- 


dinates 
a (M) = e, + Gey. 
Find the vector lines of the field. 
Solution. It is given that a, = 1, a, = 9, ag = 0. By 
virtue of formula (1) we have 


whence 


which are Archimedean spirals lying in planes parallel 
to the zy-plane. 

2°. The gradient in orthogonal coordinates. Suppose we 
have a scalar field 


U =U (%, Ys, Js). 
Then 


1 @ 0 Ou 
gradu = Fr 5q, 1+ Hy age tt Hy O05 © 
In particular, in cylindrical coordinates (q, = p, 9. = 
= 9 Ys = 2): 
1 6u 


gradu= Se e+ ete ez; (2) 


in spherical coordinates (q, = r, g. = 8, qs = 9): 


1 du Ou 
gradu=S"e,+— 00 00+ abe (3) 
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Example 2. Compute the gradient of the following sca- 
lar field specified in cylindrical coordinates (p, g, 2): 


u=p + 2zCos@. 
Solution. Using formula (2), we obtain 
1 , 
gradu=1 "€p— — 2SiN D-€y + COS H-e,. 

Example 3. Find the gradient of the following scalar 
field given in spherical coordinates (r, 6, g): 
sin8 

r 


u=r-+ —sin@cos g. 


Solution. Using formula (3), we have 
in 6 6/1 i 
gradu= (1-—) e, + —— (——cos@) eo+ 28 ey. 
3°. The curl in orthogonal coordinates. Suppose 
B=, (91, Go» 9s) €1 + 42 (Mis Yar Ys) Co + Os (Ms Jes Ys) Css 
Then i ; 


1 4 4 
TA, °' HH, ° HA, ® 
curla=| _@ _@° ee : 
041 43 09 


a,H, a,H, a;H3 


In particular, in cylindrical coordinates (q¢, = p, q. = 
=P 9s = 2): 


1 1 
> lo > & 

curla=|@ a @ |; (4) 
ap Op az 


a; paz, ay 
in spherical coordinates (q, = Tr, g. = 8, gs = Q): 


Eas oh a es 
sing" rsino °r 9% 

cula=| 2 a af, (4”) 
ar 0p 


a, ra, rsin®@-a, 
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Example 4. Compute the curl of the following vector 
field specified in cylindrical coordinates: 


€y — P2e,. 


Solution. Taking advantage of formula (4), we get 


1 1 
. 7 eo P e, 
curla=|@  @ @ |=+e,(0—0)—e,(—2—0) 
ap 0 aptCz p 


sing cos@ —pz 


+e, (0 —cos @) = zey— = 


e,. 
4°. The divergence in orthogonal coordinates. Suppose we 
have a vector field 


A=, (Gi, Ya» 9s) €r + 22 (Gi, Fer 9s) C2 + 2s (91, Yes 9s) es 
Then 


div a= 


1 O(a,HsHs) | 9(¢sHiHs) | 0 (agi H) 
HHH, [ 04, + 093 pees rae 093 ]. 
In particular, in cylindrical coordinates (q, =p, 92 = 9, 
qa = 2): 
1 4 (pa,) 1 da, | Gag. 
p to ap tas 
in spherical coordinates (g, = Tr, g3 = 9, qs = 9): 
eae | fa) 1 0(sin 8-a,) 1 das 
ava= YY or * + rain 00 + Tain 6 aq ° (5) 
Example 5. ae that the vector field 
= ieee sinve 


e+ 


is solenoidal. 
Solution. Using formula (5), we have 


1 94 2cos 8 1 ant 
diva=—4 > (7 r )trint & 38 (sin 0 = )+0 
zi (Aen) Nea g esta 0 


wherever r =£0, which means that the field of the vector 
a is solenoidal at all points with the exception of r = 0. 
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267. Find the equations of the following vector fields: 
(a) a=en+teyteys; 


(b) a= pep + pe, + 2e,; 

—— e.+ —_ 
Find the gradients of the following scalar fields. 
(a) In cylindrical coordinates: 

268. u = p*? + 2p cos p — é’ sin g. 

269. u = pcos g + zsin* g — 3°, 
(b) In spherical coordinates: 

270. u = r* cos 8. 

271. u = 3r* sin 8 + e’ cos @ — r. 

272. u=y cos’ , p=constant. 
Compute the divergence of the following vectors. 
(a) In cylindrical coordinates: 

273. a = pep + 2 Sin M-e, + e€® Cos z-e,. 

274. a = arctan p-e, + 2e, — 27e’e,. 
(b) In spherical coordinates: 


(c) a= €o, @=constant. 


275. a= re, —2cos?p-e9-++ aes ee. 


Compute the curl of the following vector fields: 
276. a = (2r + @ cos p) e, — @ sin 8-e9 + rcos 0-e,, 
@ = constant. 

277. a = re, + 2 cos 0-e9 — Gey. 


278. a=cos M-e, — at €o+ pre,. 
279. Show that the vector field 
2 0 in 0 
a= = e, +e, 


is a potential field. 
280. Show that the vector field 


a =f (re, 


where f is any differentiable function, is a potential field. 
5°. Computing the flux in curvilinear coordinates. Let 
S be a part of the coordinate surface gq, =C, where C = 
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= constant, bounded by the coordinate lines 


M1 =, Jo= A, (a, < Gq); 
93 = Bi, 93 = Be (Bi < By). 
Then the flux of the vector 


&@ = 4 (%» ar Fs) ©1 +e (915 Ye) Ys) C2 + 5 (41s Yes Fs) es 
through the surface S in the direction of the vector e, is 
computed from the formula 

@, Bs 


n= f{ J 05 (C, gar 99) Ha(C, 2s 93) Ha(C, Qo» 93) das an. 


” (6) 


The calculation is similar for the flux of a vector 
through a part of the surface g, = C or through a part of 
the surface g, = C, where C = constant. 

Example 6. Compute the flux of the vector field, speci- 
fied in cylindrical coordinates, 

a = pe, + 2, 
through the outer part of the lateral surface of the cylin- 
der p = 1, bounded by the planes z = 0 and z = 1. 

Solution. The cylinder is the coordinate surface p = 

= C = constant and so the desired flux 


2n 1 
= j j C2 dz dg = 2nC?, 
0.0 
whence for the surface p = 1 we obtain 
Tl = 2x. 


Example 7. Find the flux of the vector field, specified 

in spherical coordinates, 
a= re, + re Beg 

through the outer side of an upper hemisphere S of ra- 
dius R with centre at the coordinate origin. 

Solution. The hemisphere S is part of the coordinate 
surface r = constant, namely, r = R. On the surface S 
we have 


M=r=R; g.=8, 0<0< + ; 
g=97, 0SP<2x. 
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Taking into account that in spherical coordinates 
AL, =A, =1, H,=He=r, #3 =H, =rsin8, 
we find, via (6), that 
nm/2 an m/2 
il = j ao [ R*0 sin 0 dg = 2nRS j sin 0d0=2nR', 
c 0 0 
Compute the flux of the vector field, specified in 
cylindrical coordinates, through the given surface S. 
281. a = pe, — cos p-e, + ze,; S is a closed surface 
formed by the cylinder p = 2 and by the planes z = 0 
and z = 2. 
282. a = pe, + pye, — 2ze,; S is a closed surface 
formed by the cylinder p =1, the half-planes gp =0 and 
@ = n/2, and by the planes z = —1 and z = 1. 
283. Find the flux of the vector field a = (1/r*) e, through 
a sphere of radius R with centre at the coordinate origin. 
284. Find the flux of the vector, specified in spherical 
coordinates, 
a =re, + rsin 0-e9 — 3rq sin 0-e, 


through an upper. hemisphere of radius R. 
285. Find the flux of the vector, specified in spherical 
coordinates, 

a =r’e, + R* cos Q-eg 
through the sphere r = RA. 
286. Find the flux of the vector, specified in spherical 
coordinates, 

a = re, —rsin 0-e, 


through a semicircle of radius R located in the half-plane 
gp = n/4 (the flux is taken in the direction of the vector 


€)- 
287, Find the flux of the vector, specified in spherical 
coordinates, 


a=rsin Feotrsin 8 cos M-e, 
through the outer side of part of the half-angle cone 


V 32 = z’ + y’, bounded from above by the plane z = 
=V30<z< V8). 
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6°. Finding the potential in curvilinear coordinates. Giv- 
en, in curvilinear coordinates g,, gs, qs, the vector field 


a (M) = a (9), 92) Ys) C1 + 2 (91s G2, Ys) &s + 
+ a3 (G15 G2» Js) &a: 
this is a potential field in some region 9 over which 
the variables q,, 92, ¢3 range, that is, curl a = 0 in Q. 

To find the potential u = u (q, ge, gs) of this field, 
write the equation a al = oe u (M) as 

1 du 1 Ou 

G,@, + A2@, + a3e3 = Hi, On on a+ Zz 7m Be ey Ap ons, 3. 

From this it follows that 

du du Ou 

en Ga = Ha, gg = Cals. (7) 
This is a system of partial differential equations whose 
integration yields the desired potential u=w (q,, Ga: Ya) -+ 
+ C, where C is an arbitrary constant. 

The system (7) of differential equations is solved in 
the same way as in finding the potential in Cartesian 
coordinates 

The system (7) of differential equations is of the fol- 
lowing form: 

(1) In cylindrical coordinates (g, = p, gg = 9, 3 =2), 

7 u ? 
Fp 90, Fe = Pg, Sr = Gy. (7’) 
(2) In spherical coordinates (q, =r, ¢3 = 8, gs = 9), 
8 a, Bara, S=raind- (7) 
or or 96 3 oes 

Example 8. Find the potential of the following vector 

field specified in cylindrical coordinates: 


= a,H,, 


a= (= *+cos 9) @p— Sin 9-e9-+ 773 rere 
Solution. By formula (4) we find 
4 4 
> a) > & 
a a a 
curla= Op ap os =0 (p>0). 
= F4 . Inp 


——+cos@ —psin@ +s 
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This is a potential field. The desired potential u = 
= u(p, 9, 2) is the solution of the following system of 
differential equations: 


ou — arctan z 

ao > + cos @, | 
wee = sin 

ap a p gq, 

Ou imp 

az 142? 


From the first equation, integrating with respect to p 
we obtain 


u = In p-arctan z + pcos @ + C (Q, 2). (8) 
Differentiating (8) with respect to BA we get 
du 
3p —psing +2. 


and since du/dp = —p sin Q, it follows that dC/ag = 0, 
that is, C = C, (z). Thus 
u = In p-arctan z + p cos p + C, (2), 
whence 
& ] ’ 
= Tor tC (2). 
By virtue of the and equation of the system we have 
l l , 
Tess peat Gi) 
or C; (z) =0, whence C, (z) = C = constant. 
To. summarize, the potential of the given field is 
u (p, 9, z) = Inp-arctanz + pcosg+C. 
In the following problems, verify that the vector fields 
given in cylindrical coordinates are potential fields and 
find their potentials. 


288. a= e+ e+e, 
289. a= pep+—ey+ ze,. 
290. a= pze,+ ze, + pge,. 


291. a= e?sing-e,+ . e COS ~ -€g + 22e,. 
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292. a= @cosz-e, + cos z-€, + pp sin z-e,. 
Example 9. Find the potential of the following vector 
field given in spherical coordinates: 


a= + e%e, + sine ee, + mr peeg. 
Solution. Using (4’), we find that 
e, ree r sin 8-e, 
ts) (7) ts) 
curla=y—4| or rr) op «=| = 0. 


= e olnr-e 6 Ilnr-e% 


This is a potential field in the region where r > 0,8 s«nx 
(n = 0, +1,...). 

The system (7) of differential equations for finding the 
potential u = u(r, 6, g) is of the form 


% = pe Inr, (9) 


Integrating the first equation of system (9), we obtain 
u = elnr+C (Qg, 8). (10) 

Differentiating (10) with respect to 6 and taking into 

account the second equation of the system, we have 
pe In r = pee Inr-+-o, 


That is, 0C/00=0, whence C (9, 6) =C, (g) and 


therefore 
u =e Inr+ C, (Q). (11) 


Differentiating (11) with respect to @ and taking into 
account the third equation of system (9), we obtain 
Ge8 In r = Be Inr+ Ci (q) 
or C,(p)=0, whence C, (p)=C = constant. The 

desired potential is 
u(r, 6, g) = &Inr+C. 
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Establish the potentiality of the following vector 
fields, specified in spherical coordinates, and find their 
potentials. 


293. a= 6e, + eg. 


1 1 
294. a= 2re, + 4 ea t+ — eo. 
ag. ) 0 
295. a= 9'e,+ sno Cot eo. 
296. a=cos@sin 8-e, + cos pcos 8 -e9 — Sin P- eg. 


, 2 
297. a=e'sin O-e, ++ 6" cos 8-€9-+ Ga ae ey: 


7°. Computing the line integral and the circulation of 
a vector field in curvilinear coordinates. Suppose a vector 
field 
a (M) = a, (1, Ge. 9s) €1 + 22 (G1 92 9s) Ce 
+ ay (91, Ye, 9s) &s 


is defined and is continuous in a region 2 over which the 
orthogonal curvilinear coordinates q,, 92, 93 range. 

As we know [see Sec. 23, (8)], the differential dr of 
the radius vector r of any point M (q, 92, 9s) E 2 is 
equal to 

dr = H, dq, e, + Hy dq, e, + Hy dq; ey. 


Therefore the line integral of the vector a (M) over an 
oriented smooth or piecewise smooth curve Lc 2 is 


f(a, de) = | ay dqy-+agHedgg+ 05H dgs. (12) 
L L 


In particular, for the cylindrical coordinates g, = 0, 
q2 = 9, 73 = 2 we have 


a = dy (p, @, 2) ep + ay (Pp, , Zz) ee +a. (P, Y, 2) ez» 
dr = dp-e, + pdp-e, + dz-e,, 
and therefore 


J (, dr) = J a,dp + asp do +a, dz. (13) 
L L 
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For the spherical coordinates gq, = r, g. = 8, gy = @ we 
have 
a= a, (r, 8, @) e, + a (r, 6, P) €e + dy (r, 6, P) eo, 
dr = dr-e, + r dO-eg + rsin 8 dp -e,, 
and consequently 


| (a, dr) = | a, dr + rag d8+ ra,sin 9 dg. (14) 


The circulation C of the vector field a (M) is computed 
in the curvilinear coordinates g,, go, 73 via formula (12) 
in the general case; for cylindrical or spherical coordi- 
nates it is computed from (43) and (14) respectively. 
Example 10. Compute the line integral in the vector 
field, given in cylindrical coordinates, 
a = 4p sin p-e, + zee, + (p + @)e, 


along the straight line 


— a . 
L: ee ’ 
z=0, 


from the point O (0, 1/4, 0) to the point A (1, 1/4, 0). 
Solution. In the given example, 
a, = 49 sing, dy = 28°, a4, =p +o. 
By formula (13), the desired line integral is 


j (a, dr) = j 4p sin p dp + pze? dp + (p+ @) dz. 
L L 


On the straight line L we have 
=z, dp=0; z=0, dz=0; 0<p<i. 


Therefore 
1 


j (a, dr) = j 2V2pdp=V2 j 2odp=V2. 
L L 0 
Example 11. Compute the line integral in the vector 
field, given in spherical coordinates, 


a = e sin 0-e, + 30? sin p-eg + robe, 
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along the line 


in the direction from the point M, (1, 0, x/2) to the 
point M, (1, n/2, x/2) (Fig. 40). 


Fig. 40 


Solution. The curve L is an arc of a circle with centre 
at the coordinate origin and radius r = 1 located in the 
yz-plane. The coordinates of this vector are 


a, =e’ sin 8, ap = 30* sing, ap = rq@d. 
By virtue of (14) the line integral is of the form 


f (a, dr) = j e’ sin 0dr + 362r sin p dO + r298 sin 6 dg, 
L L 


Taking into account that the following conditions hold 
on L, 
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we obtain 
n/2 5 
a 
f @, dr) = | 30? dd = j 362d0 = 
L L 0 


Example 12. Compute the circulation of the vector 
field, given in cylindrical coordinates, 
a =p sin p-e, + p zey + pe; 
over the curve 
=sin g, 
L: {Pro o<e<e, 
directly and via the Stokes theorem. 
Solution. The coordinates of the vector are 
a4 =psing, ag =—p2z, a,=*. 
The contour L is a closed curve located in the plane z = 0 
(Fig. 41). 
(1) Direct calculation of the circulation. 


Fig. 41 


Substituting the coordinates of the vector into (13), 
we obtain 


C= § p sin @ dp + p2z dg + p° dz. 
L 


On the curve LZ we have 
z=0, dz=0; p=sing, dp = cos @ dg, 
0< gan. 
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Therefore the desired circulation is equal to 
a 


C= §psingdp= | sin? pcos gdp =0. 
L 0 


(2) Computing the circulation via the Stokes theorem. 
By the Stokes theorem, the desired circulation is equal 
to 


C= § (a, dr) = j j (curla, n°) dS, 
L s 


where S is the surface spanning the contour L. 
We find the curl of the given field: 
Cp pep &, 
1 a nr) 
curla=—- op op ra = — pep — Spey 
psing p2z 
+ (2z—cos @) e,. 
At points where p = 0 we redefine the value of curl a 
with respect to continuity, setting 
curl a (0, g, z) = (22 — cos 9) e;. 
Thus, curl a is defined throughout three-dimensional 
space. Since the curve L lies in the plane z = 0, for 
the surface S spanning this curve we take that portion 
of the plane z = 0 that is bounded by the curve L. Then 
we can take the unit vector e, for the unit vector of the 
normal n° to the surface S, that is, n° = e,. We find 
the scalar product: 
(curl a, n°) = (—pe, — 3p’e, + (22 — cos 9) e,, 2) 
= 22 — cos @ 
because by virtue of the orthonormality of the basis 
€p, €g, €, we have 
(ep, e .) = (ey, e. = 0, (e., e,) = 1. 
The desired circulation is 


C= j (2z — cos g) dS. 
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Taking into account that z=0 on S and the element of 
area dS of the coordinate surface z = 0 is equal to 


dS = p dp dg, 
we finally get 


C=— \5 cos pdS = — \J cos gp dp dp 
nx sing nm 


=— j cos pd j pdp= -+ j sin? pcos @ dp =0. 
0 0 0 
Example 13. Compute the circulation of the vector, 
specified in spherical coordinates, 
a =re, + (R +r) sin 0-e, 
around the circle 


in the direction of increasing values of the angle q, 
directly and via the Stokes theorem. 
Solution. In this example, 


a,=r, a =0, ag=(R+7r)sin6. 


(1) Direct computation of the circulation. 
By formula (14) the desired circulation is equal to 


C=Ordr+(R+r)sinOr sin Odp= 
E 
=rdr+r(R+r) sin? 6 dq. 
L 


On the given circle L, the centre of which lies at the 
coordinate origin, we have 
r=, dr=0; O=+; 0<g<2a, 


and, consequently, 


C =2R° dpa on? | dg Anh 
L 0 
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(2) Computation of the circulation via the Stokes 


theorem. 
The desired circulation is, by the Stokes theorem, equal 
to 


C= § (a, dr) = j (curla, n°) dS, 
L 8 


where S is the surface spanning the circle L. 
We find the curl of the given vector: 


e; reg rsin@-e, 
cirla Se oO oO 50. 
~ risinO | dr 3 ap 
r 0 (Ar+r?)sin?6 


=2 (R +r) cos @-e,—+ (R+ 2r) sin 0-e9, 


For the surface S spanning the circle L we take, for 
example, the upper hemisphere of radius R: r= R, 
0< 8 < xn/2,0<q@ < 2n. The unit vector of the normal 
n° to the outer side of the hemisphere S is directed along 
the vector e,, and so we take n° =e,. We find the scalar 
product 
R+2r 


r 


2(R+r) cos 8-e. — 
r r 


(curl a, n°) = ( sin 8-e9, er} 


= AMT) cos 6 
since (e,,e,) = 1, (,, ee) = 0. 

Taking into account that r = R on the surface S, we 
obtain the following expression for the desired flux: 


n= {| 2E=" cosdds = 4) | cosOds. 
“s s 
In spherical coordinates, the element of area dS of the 


coordinate surface r = R, that is, the hemisphere S, is 
equal to 


dS = R* sin 0 d8 do 
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and, consequently, 


n=4{ j cos 8 R? sin8 d6 dp 


§ 
1/2 2n 
= 4Re{ cos @ sin @ 0 | dy = 4nR?, 
0 0 


By the Stokes theorem we get C = 4nR?. 

For the surface S spanning the circle L we could take 
the lower hemisphere, the unit vector of the normal of 
which, n° = —e,, and the result would be the same: 
C = 4nR?. 

Note that it is not desirable to take the circle bounded 
by Z for the surface S spanning the circle L because the 
circle has a point r = 0 (the centre of the circle) in which 
the curl of the given vector has a discontinuity. 

Compute the line integral along the given curves L in 
the following vector fields specified in cylindrical coordi- 
nates. 

298. a = ze, + ppe, + cos g-e,; L is a segment of the 
straight line: {9p =a, p=0,0<2< 4}. 

299. a = pe, + 2pqe, + ze,; L is the semicircle: {p = 1, 
z=0,0< qn}. 

300. a =e? cos g-e, + p sin g-e, + pe,; L is a turn 
of the helical curve: {op = R, z= 9, O< @< 2a}. 

Compute the line integral over the given curve L 
in the following vector fields given in spherical coordi- 
nates. 

301. a = e’ cos 6-e, + 26 cos g-eg + ge,; L is the semi- 
circle: {fr = 1, p=0,0< 0 < a}. 


3uZ. a=4r° tan + e, + Opes + cos? p-ey; Lis a segment 
of the straight line: {o=+; d=, 0<r<i}. 
303. a=sin? 0-e,+sin 6-e9-+rq8e,; L is a segment of 


the straight line: {e=F. r=, + <0<$}. 


Compute the circulation of the following vector fields, 
specified in cylindrical coordinates, about the given con- 
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tours, directly and with the aid of the Stokes theorem. 
304. a = ze, + pzeg + pe,; L is the circle: {p = 1, 
z = 0}. 
305. a = p sin g-e, — p’ze, + p’e,; L is the circle: 
{o = R, z = R}. 
306. a = z cos g-e, + peg + g’e,; L is the loop: {p = 
= sin g, 2 = 1}. 

Compute the circulation of the following vectors, given 
in spherical coordinates, along the given contours L, 
directly and with the aid of the Stokes theorem. 


307. a = re, +rsin 6-e,; L is the circle: {r =1,0= 
= 7} 

=F}. 

308. a = rsin 8-e, + @e%,; L is the loop: {r = sing, 
6 = +: 0<qog x} : 

309. a=rqe,; L is a contour bounded by the semi- 
circle: {r=R, e=>) 0<e<z} and its vertical 


diameter {o = > 6 = o} : 


Sec. 25. The Laplace operator 
in orthogonal coordinates 


If u = u (q, Ye, Qs) is a scalar function, then 


1 du du Ou 
rad u=—- —— 3. 1 
4 H, 8g; eC et+sz- an 095. €2 +z- nr O93. 3 ( ) 


If 
a=, (9), Ger Ia) 1 + Fe (9x5 Ger Ja) C2 + Gs (Qi, Ya» 9a) Cas 
then 


SZ (a\HaHs) + 2 (a,HsH,) 
1 2 


div a= sor 
H,HyH, | 09 


+ F(a H,)]. (2) 
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Using formulas (1) and (2), we obtain the following expres- 
sion for the Laplace operator Au: 


i 1 HH, du 
Au = div grad U=aaa- lee 2fts 
e HHH, ae | Hy, dq, 
4 2 (Aas Ou 2 (2M = 
9s Hy 993 995 Hy 445 


In cylindrical coordinates, 


1 t) du J] 1 ou 
du=— [ap (0 a) tay (5 ag) tae (0) 
4 @u , uy 
) p? aq? 028 * 
In spherical coordinates, 
4 a ee, 
Au=5— (+ (resine 2) +3 (sin 05) 
a 1 du ar) du 
+35 ay (aE a) Je ae (" oa 
1 ; du\ . 1 G7u 
+ 75nd x (sin 0) T Fisin?0 ag? * 

Example. Find all the solutions of the Laplace equation 
Au = 0 that depend solely on the distance r. 

Solution. Writing the Laplace equation in spherical 
coordinates and taking into account the spherical sym- 
metry of the solution (it must not depend on 9 or 9), 
we have 


1 @ ts) 
dues (rt)=0 (=u), 
whence 
ou 
r= Cy 
so that 
u=2+¢,, 


where C, and C, are constants. 
310. Given: a scalar field u = u (M) in cylindrical coor- 
dinates 
u (p, 9, 2) = pp + 2q* — pqz. 
Find Au. 
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311. Given: a scalar field uv = wu (M) in spherical coor- 


dinates 
u(r, 8, 9) = rg + Pq’? + op + 8. 
Find Au. 
312. Are the following functions harmonic? 
(4) u = p? cos 29. 
(2) u = rcos 20. 
313. Find all possible harmonic functions that 
(a) depend on 6 alone, 
(b) depend on @ alone 
(in the spherical system of coordinates). 
314. Find all solutions of the Poisson equation 


Au = r®-1 


in the spherical system of coordinates, provided u = u (r). 


ANSWERS 


; z=2, 
1. (a) The half-line y=~—z,y>0, z<0 is traversed 


twice when—co<t< + oo. (b) wet ass — 


U(—1, +00) the point r(t) = ia: 


war TF i 
twice traverses the half-line x+y =1, sot, y<y. 
) z 
al ae (d) y= 2=>. 
(©) z=1;)° 
(e) z?+-y2+22=1, r—y=0. 
7. i+k. 8& i+k. 9 —j+— k 
10. —i+k. 11. ei—j+2k. 12. No. 14. No. 
dr . 2 dr \. ar 
17. (a) 2(Fr) +(r, Sx)i (c) [r. ar |: 


21. Circles lying in planes perpendicular to the vector a. 

22. The hodograph of velocity is a helical curve: z = 

=acost, y=asint, z = 2bt; the hodograph of accele- 

ration is a circle: zr = —asint, y = acost, 2 = 2b. 
da da du , dia d'a / du da dtu 

26. =e ae ae = aur (ae) tae a 


28. (¢—A)eti+ + (tJ sin 2) j—arctant-k+c. 
29. Sin (1-12) -i+ + ef.j+sin t-k+ec. 


30. esint.j_— 4 sin 2.jtitk+e. 
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3 t 
31. £ i+ (tcost—sin tj + —pkte. 
2 


32. zs itak. 33. (1 —e-'/2) i+ (e!/?—1) j++ (e—1) k. 


; 2 
34. —In2-j+k. 35. 2n%+ nj+n%k. 36. R= tS ; 


37. R=— |t|(14+92)92, 38. R=6. 39. R=+an. 


40. R=2acosh?t. 44. z+y=0. 


5 1 1 
42. zc—y—V 2z=0. 43. T=T° 44. 


i eee: eee 
T 2acosh?t ° 


45. $4 H4 526, a family of triaxial ellipsoids. 


46. x? ++ y? —z=C, a family of paraboloids. 

47. x? + y* = Cz, a family of paraboloids. 

48. 2y? + 922 = C, a family of elliptical cylinders. 
49.x2+2y—2z=C, a family of parallel planes. 

90. A family of planes resulting from the sheaf of planes 
az + a,y + asz = C (bz + bey + 552) passing through 
the straight line 


a,r + ay +a,z=0, 
bz + boy + bz = 0. 

via the elimination of the straight line itself. Here, a,, a,, 

ag are the coordinates of the vector a; b,, by, b3 are the 


coordinates of the vector b. 
51. 22 + y? + 2? = C*, a family of concentric spheres. 


zy 2 

52. (a, b, r)=C or |a, a, a3}=C, a family of parallel 
by by by 

planes. 

33. 22 — y =C, a family of parallel straight lines. 

94. y = Cr, C>0, z #0, a family of rays. 

55. y? = Cz is a family of parabolas with the vertex 

O (0, 0) deleted. 

56. 2? — y? = C, a family of hyperbolas. 


57. 
58. 
62. 
67. 


70. 
73. 


75. 


81. 


87. 
88. 


90. 
91. 
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y = —zInC —C, C>0, a family of straight lines. 


V¥% 59 —V2 go, V3 os, sa. 


23, 63.4. 64. 0. 65. 2¥3 (y343). 66. 0. 


~2, 68, —S_., 69. S(i+j—h). 


k. 71. p=a. 72. p=0. 
og = 0.74. y = —z + Qnn, n= 0, +1, +2,.... 


at+y?+22=1. 78>. 79. a. 80. a(b, r) +b (a, 1). 


2|a|2r—2 (a, r)a. 86. ee 


or r° 
a dp. a 
lu cos (r, 1) . “i 
z= zZ=9 for r_Ll. 
Ou 4 du 
aa 8 a= 


du _ (gradu, grad vy), du_o. 
a Tgred vf iF =0 if grad ut grad v. 


(a) 1 in the direction of the y-axis; (b) 3 in the 


direction of the vector a= —i—2j+ 2k. 


92. 


94. 


y=e,r; z=e,2. 93. y=rtey; z= sh atCe. 
1 


4 4 1 1 
22+ y2+ 22=Ci, } 95. soya ztar=4 
zt+y+z=C,. 


96. v=Cy, 2=C,. 97. 2=Cyz, y=Cy. 98. zry=C,, 2=C,. 


99. 


101. 
103. 


T=Cy, 2y2—22=C,. 100. 4-{=¢, 2=C. 
y2+22=C,, z=C,. 102. r=Cyy, r=C,2. 
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where bo, Bos, 593 are cooordinates of the vector by. 
104. TI = —3. 105. I] = xR*y. 106. Il = nReh. 


107. = 4nR5f(R). 108. 1=+. 109. =4. 
110. M=+nReh. 111. = nhs. 

81 1 
112, =n, 113. T=, 


144. M=0. 145. 1=2, 116. n=. 1417. = 40k, 


2 

118. (a) M=>—2; (b) N=—F; (c) N=n. 119. =O. 
120. W=6nR. 121. 1=0. 

v2 
122, =m. 123. M=0. 124. 1=tn(1—-4*). 
125. 1=2 Rt. 
126. W=/V2n. 127. W=45n. 128. n= Sq. 
129. 1=0. 


130. =n. 134. p(r) =<. 132. Ir. 133, 0. 134. 0. 


135. p (z) = C — z, C = constant. 136. Il = 4nR°. 
137. divE = 0 (r #0). 


143, 16x. 144. H3. 145. — 146. 0. 147, 2. 
148. 4m. 149. on. 150. on. 151. 2R3. 


152. an. 153. —1. 154. — 


155. Solenoidal field. 
156. Nonsolenoidal field. 
157, aicuuaes field. 


159. g(r)=—, r>0, C=constant. 
161. ion 162. In, 163, ++, 164. 0. 
1 
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4 4 14 2 
166. —2R% 167, . 168. (a)—42; (b). 169. 0. 


170. 5. 171. 3/3. 172. 
173, —na?. 174. 1. 

175. —2n. 176. —2@ | 177. 4, 

179. —2(zi+-zj + yk). 

180. 3 (2222) j. 181. (2+y)k. 191. = curly = —4j, 


193. f(z, z)=22+2+2+C, C=constant. 
195. 4n. 196. —4n. 197. 4. 198. —2n. 199. x. 
200. 729. 201. 0. 202. —V2n. 


203. 2wna?. Hint: v=[o, r]. 
204. p.=1. 205. p,=3. 206. Dependent. 


207. Independent. 208. Dependent. 209. —1. 210. <<. 
2 n 1 
2it. 0. 212, 2. 213. 2. 244. 2. 


1 
35 ° 


246. +. Hint: Supplement the path of integration L 


with the line segment OA of the z-axis. 

217. No. 218. Yes. 219. No. 220. Yes. 221. No. 222. No. 
223. Yes. 226. p = z’yz. 227. p = 2x + zyz. 228. 9p = 
=2y—y+az. 229.9 =Injr+yt+2z|. 20. 9 = 
= arctan (zyz). 231. 9 = r. 232. @ = Inr. 


233. gaan. 234. p=axr+ fBy-+-yz+C, C =constant. 


235. p= sy tyztartCc. 236. p=aytet+C. 
237. p=e*siny+2z2+C. 247. (a) Yes, (b) No, 
(c) Yes. 

249, u= Cr+ C,. 

250. u = Az? + Bry + Ay’, where A and B are arbi- 
trary. 


wapr tts if n=~1, 


251. ve=4 
zlin [z]/+C,z+C, if n=1(r¢ 30). 
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252. T= a. 253. [= —*. 254. [= on’, 


3 
255. I = 0. 
256. b = zj + (y — z)k. 257. b = (y? — 2zz)k. 
258. b = (e* — ze”) j. 259. b = 32°j + (2y° — 622) k. 
260. b= —z(z+yY)j+@+y)k. 


261. b = — (xz? + yze**) j — Qryzk. 
262. b= (—2j+-yk). 

263. b = —8yzi + 2z2j + 7zyk. 

264. b = Qry*zi — 3z%yzj + zy*k. 
265. =+sin za-i—+ sin zz-k. 


266. b= atu z—k. 267. (a) p=Q+Cy, p=z+Cy; 


(b) p= ae, P= Cuts (0) P= Cy r=Cz sin. 


268. grad u=2(p+cos @) ey — (2sin 9+ —e' cos g) Cy 
—e’ sin -e,. 
269. grad u=(cos@—3* In 3)e,+ (4 sin 2p—sin 9) ey 
+ sin? p-e,. 
270. grad u=2rcos@-e,—r sin 8-eg. 
271. grad u=(6rsin® +e" cos gp —1) e, + 3r cos 8 -e9 


e’ sin @ 
rsin@ ~* 


272. grad u=—p (2S e rae sin’ ea). 


273. div Soo eg wains 
ivaat a Soe t 
274. div a= . arctan p+ T+? (22+ 22) e’. 
- 2 
275. div a= 4r— cost @cot8+ arg. 


__ cos 20 sin @ asin 6 
276. curla=—T>-e r— (20080457 % sine <a y) ee Cy. 


277. 
278. 
281. 


285. 
287. 
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curl a= — cot 0- e+e + 20088 Cy. 


curl a= —2pe,-+ —— sin oe 

24n. 282. $n. 283. 4x. 284. = n°. 

4nR‘. 286. —< R, 

48. Hint: Write the equations of the surfaces in 


spherical coordinates. 
28.u=p+qo+24+C. 


289. 


U=T(P+—+2)4C. 290. u=ppz+C. 


-u=ePsngt2+C. 292. u=ppceosz+C. 


293. u=r04+C. 2944. u=r+94+64C. 
eae gp + 6?) ++ C. 296. u = rcosqsin® + C. 
297. u = Coens # irs 

298. 1, 299. n®, 300. 2nR. 301. n®. 302. 4. 

303. a4 v8 

304. 0. 305. —2nR*. 306. x. 


314. 


. (4) u(0)=C, In 


. mt. 308. 0. 309. 0. 
: Au =49—+ om +298, 


. Au= = 698 + 12rg?+ = 4+ cot0+% cote+ ——— 
. (1) Yes, (2) No. 


eae (2) u(9) =C,p+Cr, 
pnt 

arnee te Cp an 

u(r) = Inr + —14+C¢,, n=—1, (r0) 


( ~L 4546, n=—2. 


APPENDIX | 
BASIC OPERATIONS OF VECTOR 
ANALYSIS IN ORTHOGONAL 


CURVILINEAR COORDINATES 


1. The scalar field is given in orthogonal curvilinear 
coordinates, u = wu (q;, gz, Ys). Then we have 
1 du 1 Ou 1 du 
ad u=-—- ——e,+-—-~—e,+ —-—— 
gradu Ay 09, 1TH, 092 ee Hs 095 
The Laplace operator is 
i A (oO ais Bey | Bas oe 
cae H,H3H, 3a ( A, a) ton ( Hy 095 
6 { A,H, du 
+ q5 ( Hy ale 


Special cases: (a) The scalar field is given in cylindrical 
coordinates, u = u (p, g, z). Then is nai 
Ou 1 du 

grad u=-—5— en +> . oD Pas, are 
The Laplace operator is 
ey 1 au, Hu, 
du=> 5 (0 e) tor age oa 
(b) The scalar field is ae in spherical coordinates, 


u= u(r, 9, 9). Then we have 
gradu=e,++ : = eto ee re 

The Laplace operator is 
=e ("3) + oramve ae (Ge 8" 8) 


Au= aa 
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2. The vector field is given in orthogonal curvilinear 
coordinates: 


@ = 4, (9, Ja, Ys) 1 + 2 (G1, Yer Fs) Ce 


ei as (a1, qe) qs) e3. 
Then we have 


A 1 4 (a,H,H,) 9 (as yHs) 9 (G3 Hs) H;) 
div a= H,HAH, [ 09; + 043 2493 a) 
H3H, ' H,H, * H,H, ° 
curl a= ee. a 9 
04, 093 043 


a,H, a,H, asfl; 


Special cases: (a) The vector field is given in cylindrical 
coordinates, 


a=Qy (p, q, 2) @p + a, (e, q, 2) €y + ay (p, q, z)e;. 
Then we have 


1 @ (pa) 1 das Bs 
V3.5 Tt ’ 
ar aay aT ot ee 
1 
p e@ —e&, 
curl a=|—& oy 
"| op ap Oz 
a, Pa, a 


(b) The vector field is given in spherical coordinates: 


a = a, (r, 8, p)e, + a, (r, 8, @) eo + ay (7, 9, @) ep. 
Then we have 


_ 1 @(a,r?) 4 8 (az sin 8) 4 Gay 
div a=—3 a + yam a0 + Tsind a9” 
1 1 é é 
resin rsino ° r@ 
eta | ee a. ee 
4 or 00 CT) 


a, ra, agr sin® 


APPENDIX II 


AREA ELEMENTS 
OF COORDINATE SURFACES 


Coordl- 
nates 


Coordinate surfaces 


Area elements 


General | 7! C=constant |dS,=H, (C, 99.93) Hs (C, 99593) 493 agg 
qg=C=constant |dS,=F, (q,,C. 93) Hs (41: C. 9s) 49, day 
G19 Yar I3 qg=C =constant |dSs=Ay (9), 92) C) He (415 92s C) dq; Bug 


Cylind- 
rical 
9,=pP | P=C =conslant 
I2= P| P=C=constant 
93 =2 | z=C=constant 
Spherical 


q: =r | r=C=constant 
q2=9 | 8=C=constant 
I=) p=C=constant 


dS =C dg dz 
dS = dp dz 
dS =e dp dp 


dS = C* sin 6 d6 dp 
dS =r sin C dr dp 
dS =r dr d8 


~-& WS w= 
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dents. The book contains 967 problems a con- 
siderable proportion of which have been com- 
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problems are provided with answers. There are 
also some exercises of a theoretical nature. Each 
section begins with a summary of basic facts 
and comprises worked-out examples of typical 
problems. The book concludes with two useful 
appendices. 
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of theory (pr 


with a brief review of the Sean 


os 
ila: . This is followed b oe 
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| are 


The book contains 314 problems to be 
solved by the student. They are all 
provided with hints and answers. The 


41 drawings that accompany the text 
serve to help the reader in analysing 
the theoretical material eee problems. 


